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ABSTRACT 


The purpose of the present study is to develop a most general structural dynamic 
model for a flexible rotor blade undergoing coupled axial, flap, lag and torsional 
deformations. The dynamic model includes the complex geometric parameters, such 
as pretwist, predroop, precone, presweep and torque offset, as well as hub motion. 
The ecpiations of motions are derived using the Hamilton’s principle. Using a beam 
type finite element model for the blade, the mass and stiffness matrices are obtained 
respectively, from the kinetic energy and the strain energy expressions. 

The structural dynamic model derived in this study was then used to calculate the 
undamped natural frequencies of a rotating blade. These results are compared with 
the data available in the literature. 

This dynamic model will be useful for the study of active control of vibration in 
helicopters as well as for rotor blade acroelastic analysis. 
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Chapter 1 


INTRODUCTION 


Tlic reduction of vibration levels in helico])t(us is one of the key problems in the 
design of helicopter. Increasing demand towards expanding the mission range of 
helicopter, such as high speed, high ‘g’ manoeuvres , coupled with improving the 
system reliability and reduced maintenence cost, has resulted in more stringent 
vibration specifications. Hence, the reduction of vibratory levels below an acceptable 
limit has become a critical design consideration in modern helicopters. Over the 
years, vibratory levels in the helicopter cabin have been reduced by suitable design 
modifications and also by incorporating vibration absorbers and isolators. However, 
in the next generation helicopters, with the stringent vibration control, it will become 
necessary to reduce the vibration levels below 0.05 g or even 0.02 g. 

It is well known that the time varying loads on the main rotor system contributes 
significantly to the vibration in helicopter. Therefore, the structural dynamic char- 
acteristics of the rotor blade and also the dynamic characteristics of the fuselage 
have a very strong influence on the vibratory levels in helicopters. Any analytical 
study of vibration in the helicopters requires the developement of a dynamic model 
of the helicopter system. The major components of this model are: 

• Rotor blade model 

• Fuselage model 
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• ll()tor-t'\isclage interface model 

The formulation of the rotor blade model requires the development of structural, 
inertial and aerodynamic operators associated with the rotor blade motion. The 
helicoptor fuselage model is represented by an idealized structural dynamic model 
of a three-dimensional structure. The rotor- fuselage interface model must represent 
both the geometry of the interface as well as the aerodynamic interaction in an 
appropriate manner. 

This study deals with the formulation and solution of a structural dynamic (struc- 
tural and inertia operators) model of a rotor blade incorporating all complex geo- 
metric parameters of a realistic rotor system. 


1.1 STRUCTURAL MODELLING OF ROTOR 
BLADE 

Most of the published work on the structural modelling of the helicopter rotor blades 
assume a beam type model. Nonlinear beam kinematics, which incorporate small 
strains and finite (moderate or large) rotations, are being used to account for the 
coupling effect between axial, bending and torsional deformations. In the derivation 
of the strain-displacement relationship, a small strain assumption (i.e. strains are 
small as compared to unity) is made. Such an assumption is consistent with the 
design requirement bfised on fatigue life consideration which states that the rotor 
blades must be designed to have an operating strain level well below the elastic limit 
of the blade material. 

The first analytical model for the flap-lag-torsion of pretwisted nonuniform rotor 
blades was developed by Houbolt and Brooks[l]. This model is based on a linear 
theory which neglects all the nonlinear displacement terms in the derivation. As 
a result, the bending-torsion coupling effects due to the geometric nonlinearities. 
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which arc important for the rotor blade analysis, are absent in the model. 

In order to incorporate the geometric nonlinearities due to the assumption of small 
strains and finite rotations, one should distinguish between the deformed and unde- 
formcd configurations of the blade, and derive the transformation relation between 
the triad of unit vectors associated with the undeformed configuration of the blade 
and the triad of unit vectors associated with the deformed blade configuration. In 
the moderate deflection beam theories [2-7], the transformation relations between 
deformed and the undeformed triad of unit vector were derived. These beam theories 
were later used to formulate the structural, inertial and aerodynamic operators for 
tlie aeroelastic stability analysis of rotor blades [8-12]. Structural modelling for 
composite rotor blade has also been developed recently by various investigators[13- 
22]. The earlier models presented in Refs. 2-7, were restricted to treatment of 
isotropic blades. In general, these models did not include the effects of cross- 
sectional transverse shear and warping. On the other hand, the composite blade 
models include these elfocts. 

During the derivation of the moderate deflection beam theory, a number of nonlinear 
terms are generated. Many of them are relatively small due to the assumption of the 
small deflections and moderate rotation. Therefore, an ordering scheme is needed 
to identify and neglect higher order nonlinear terms in a consistent manner. The 
commonly used ordering scheme [23,24] is based on assigning orders of magnitude to 
the nondimensional physical parameters governing the aeroelastic problem in terms 
of blade bending slopes, which are assumed to be of the order of e. A second order 
approximation assumes that the terms of order are neglected compared to the 
terms of order 1. In this study a third order approximation has been used, i.e, terms 
of order are neglected compared to the terms of order 1. i.e. 

1 + O(e^) 1 

In all the studies reported in Refs. [2-25], only certain geometric parameters such 
as precone, hinge offset and tip sweep have been included in modelling the rotor 
blade dynamics. In addition, most of these models, except [24,25], do not consider 
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hub motion. In the present formulation, an attempt has been made to include all 
the complex geometric parameters such as torque offset, blade root offset, precone 
angle, predrooi) angle, i)resweep angle aiul pretwist, as well as hub motion. 


1.2 OBJECTIVES 

The objectives of the present study are: 

• To develop a most general model for the rotor blade including various geomet- 
ric parameters like root offset, torciue offset, precone, predroop, presweep, and 
pretwist and hub motions. 

• To validate the model by comparing the results of this study with those 
available in the litrature. 



Chapter 2 


MODEL DESCRIPTION AND 
ORDERING SCHEME 


Helicopter rotor blades are long slender beams attached to the hub through a 
mechanical arrangement. There are several geometrical parameters describing the 
configuration of the rotor blade attachment to the hub. A schematic of the rotor 
configuration is shown in Fig.l. Figure 1 shows the torque offset ‘a’ which is the 
distance from the centre of rotation (hub centre) to the pitch axis of the blade, 
‘ei’ and ‘ 63 ’ refer to the blade root offset distance from the centre of the hub. (3p 
represents the precone angle defining the orientation of the blade pitch control axis 
with respect to the hub plane of rotation, pg is the presweep angle defining the 
orientation of the blade in the plane of rotation. j3d represents the predroop angle 
representing the inclination of the blade with respect to the pitch control axis of the 
blade. 


2.1 BASIC ASSUMPTIONS 

In the formulation of the dynamic model of the rotor blade, several assumptions 
have been made which are given below: 
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1. Tlu^ blade is t.real.(Hl a.s aa elastic Ixiaai. 

2. The blade is modeled by a series of straight beam finite elements located along 
the elastic axis of the blade. 

3. The speed of rotation (f2) of rotor is constant. 

4. The rotor shaft is rigid. 

5. The blade undergoes moderate deformation. 

6. The blade cross section can have a general shape with distinct shear centre, 
torque centre and the centre of mass. 

2.2 ORDERING SCHEME 

In the formulation of the equations of motion of a rotor blade undergoing moderate 
deformations, a large number of higher order terms are generated. In order to 
identify and eliminate liigher order terms in a consistent manner, an ordering scheme 
is employed. This ordering scheme is based on the assumption that the slopes of the 
deformed elastic blade are moderate, and of order e (0.10 < e < 0.20). Orders 
of magnitude are then assigned to various non-dimensional physical parameters 
governing the rotor blade dynamic problem, in terms of e. In the derivation of 
the governing equations, it is assumed that terms of the order are neglected with 
respect to terms of order 1, i.e., 

0(l) + 0(e^) ?«0(1) 

The order of magnitude of various non-dimensional parameters governing this prob- 
lem are given below ; 

ORDER 1 


cos.^i.sini/'t, y = 0(1) 

il,( ) = *()= 0(1) 
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l—i ) = — 

^ dxk ^ ' dxk 


0 ( 1 ) 


ORDER 

Ogk = 

ORDER e 

— n(f) 

I' I ' I ' n n I ■> i ~ 

Vk,w'k^^,Pp,Pd,Ps = 0{e) 

ORDER 

Im,,,,, = 0((:'^\^) 

ORD ER. 

^,u'k,rnr)m,mCm-0{e^) 

It is important to note that, ordering schemes are based on physical understanding 
of the behaviour of actual blade configurations. Hence care must be excercised in 
deleting the higher order terms, based on this ordering scheme. 



Chapter 3 


COORDINATE SYSTEMS 


The description of the complex deformation of a rotor blade requires several coor- 
dinate systems. The transformation relation between quantities referred in various 
inertial, non-inortial coordinate systems is to be established before deriving the 
equations of motion of the rotor blade. The relation between two orthogonal systems 
with axes Xi,Yi^Zi and Xj,Yj,Zj with exi,eyi,ezi and as unit vectors 

along the respective axes can be written as 






^xi 



£xj 


li 

ij] < 

* 

^zi 

✓ 





(3.1) 


where the transformation matrix [Tij] can be obtained using the Euler angles re- 
quired to rotate the j-system so as to make it parallel to i-system. The coordinate 
systems used in deriving the equation of motion for the rotor model are described 
below: 


3.1 INERTIAL SYSTEM -R 

The coordinate system ‘R’, shown in Fig.2, has its origin at the centre On of the 
unperturbed hub. The xr axis is pointing towards the helicoptor tail and zr is 
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pointing upwards. The unit vectors along the three axes are e^/^, tyR, e^R . 


3.2 BODY FIXED HUB SYSTEM -H 

The coordinate system H shown in Fig.3, is a body fixed system with its origin fixed 
at the centre of rotor hub Or. Prior to perturbational motion of the hub, the H 
-sj^stem coincides with R -system. 

If Oy and Oz represent the roll-pitch-yaw rotations, then the transformation matrix 
[Thr] can be written as - 



1 

0 

0 


cos 9 y 

0 

— sin 9 y 


cos 9z 

sin 9z 0 

[Trr] = 

0 

cos 9x 

sin 9x 


0 

1 

0 


— sin 9z 

cos 9z 0 


0 

— sin 9x 

cos 9x 


sin 9 y 

0 

cos 9 y 


0 

0 1 


I- J J L. J 

(3.2) 

Since dx, 9y and Qz are assumed to be of order the sine and cosine functions can 
be approximated as 

sin ^ « (9 


and 


cos 9 TH I 


Hence - 


[P77i^] = 


1 9 z — 9 y 

9y9z - 9z 1 9x 

9x9z + 9y 9y9z — 9x 1 


(3.3) 


3.3 ROTATING HUB FIXED SYSTEM -1 


The coordinate system 1 shown in Fig.4, rotates about zr axis with constant 
angular speed Q of the rotor. Its origin is fixed at hub centre Or. This system 
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can be obtained by rotating H system by an azimuthal angle ijjk of the blade 
about zji axis. The transformation matrix is given as - 

cos ‘il>k sin V-'fc 0 

[^17/]= -sin^fc cos-i/jfc 0 (3.4) 

0 0 1 

where, ipk is azimuth angle of blade 

Stt 

ipk = ^t + {k- 1) — 

■i^b 

3.4 ROTATING SYSTEM -2K 

The coordinate system 2k shown in Fig. 5 is a blade fixed coordinate system which 
rotates with blade. The origin of the 2K system is at the location of the blade 
root, which is at a distance aCyi+CiCx] from the hub centre. 1-system and 2k-system 
are parallel. 

3.5 PRECONED, ROTATING SYSTEM -3K 

The system 3k shown in Fig.6a and Fig.6b also rotates with blade. This system is 
obtained by rotating 2k -system by an angle — /?p (precone angle) about y^k axis. 
The transformation matrix between 2k and 3k systems is given as - 

1 0 /?p' 

[T32]= 0 1 0 (3.5) 

_ -idp 0 1 _ 

3.6 PREDROOPED, PRESWEPT, PITCHED, 
BLADE-FIXED ROTATING SYSTEM -4K 

The 4k system shown in Fig.fia and Fig.6b, is blade fixed system with its origin at 
pitch bearing of the blade . 
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It may be noted that the pitch axis of the blade is along direction and the blade 
reference elastic axis is along the direction. During the control pitch input, the 
clastic axis will move on surface of a cone whose vortex angle is described by the 
angles /?,, and as shown in Fig.6c. 

The 4k system is obtained by the following steps - 

Step-1 Translating the origin of 3k system by a distance ‘62’ along Cxzk- 

Step-2 Then rotating the system by an angle —Ps (presweep angle) about zzk axis. 

Step-3 Then rotating the system by an angle —pd (predroop angle) about yzk axis. 

Step-4 Then rotating the system by an angle 9i (pitch input) about Xzk axis. 

The transformation matrix is given as - 


1 




(/I, cos!?/ 4- /?,/sin Oj) 
{-Pd cos 0] -f /), sin 0i) 


—{Ps cos 01 -H Pd sin 0i) {Pd cos 0f — Pg sin 0[) 
cos 0j sin 01 

— sin 0] cos 01 

(3.6) 


3.7 ROTATING, BLADE-FIXED SYSTEM -5K 

The 5k system, shown in Fig. 7 is a cross-sectional coordinate system of the 
blade. In the undeformed state of the blade, both 4k and 5k systems are parallel. 
But the origin of the 5k system is at a distance Xk from the origin of the 4k system. 
During elastic deformation of the blade, the 5k system translates and rotates with 
the cross-section. After the deformation, the origin of the 5k system, from the origin 
of 4k system, is at the location given by 


{"^ r?. 1 ^ Zg -|- Xk T '(^^fc)Cx4fc d” '^k^yAk T '^k^xAk 


The transformation matrix between 4k and 5k system is obtained following a flap 
- lag sequence of rotation. The Euler angles are respectively -Pk and ^k which 
correspond to the local slope of the deformed blade. The transformation matrix is 
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given as: 



cos ^k sin 0 


1 

a 

0 

0 

1 

11 

- sin ^k cos ^k 0 


0 1 0 


1 

r— 1 

0 

1 


- sin Pk 0 cos pk 


(3.7) 


Since the angle and / 3 k are of the order 0(e), the transformation matrix can be 
simplified by assuming - 

sin 0 Ri 0 

and 

cos 0^1 

Defining the Euler angles in the terms of elastic deformation of the blade 

dwk 


Pk = fi’k 
— 't^k 


dxk 

dvk 

dXk 


Substituting the above relations in the transformation matrix [T54] yields 

1 

-w[ 


[Ts 4 ] = 


Vk W'k 


•' 1 




A 0 


(3.8) 


3.8 COORDINATE SYSTEM -6K 


The 6k system shown in Fig.8 represents the cross-sectional coordinate system in 
the deformed configuration of the blade. e,j,e^ represent the directions of the cross- 
sectional principal axes. 6k-system is obtained by rotating 5k system about e^^^k 
through the angle {(j)k + Og)i where 9 g represents the geometric twist angle of the 
cross-section and <j)k represents the elastic twist. 

The transformation relation is given as - 




1 0 0 


Cx 

1 ^1/ 

>. = 

0 cos{6a + (f>k) - sm{Oa + Pk) 

< 



5k 

0 sin(6'G + Pk) cos{9g + pk) 




6 jfc 
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Summary of Coordinate Systems 


Symbol 

Origin 

Coordinate System 

Unit vectors 

R 

Hub centre 

Oh 

Inertial. Fixed at undeformed 

hub centre ‘0^’ 

^yRi ^zR 

H 

Hub centre 

Oh 

Noninertial , body fixed. 

^zH 

1 

Hub centre 

Oh 

Rotates with blade with 

angular speed Q. 

^rrlj ^zl 

2k 

At point ‘A’ 

shown in Fig.6. 

Rotates with blade. Systems 

1 and 2k arc parallel. 

^x2k^ ^y2k} ^z2k 

3k 

At point ‘A’ 

shown in Fig.6. 

Preconed, rotates with k^^ blade. 

^xZki ^yZk^ ^zZk 

4k 

out board of 

pitch bearing 

Predrooped, preswept, pitched blade 
fixed system. Rotates with blade. 

^xAk} ^y4k^ ^z4k 

5k 

At elastic axis 

Rotates with blade. Deformed 

blade coordinate system. 

^x5kj ^y5ki ^z5k 

6k 

At elastic axis 

Rotates with k^^ blade. 

Obtained by rotating 5k system 
through an angle {Oq + (j>k) 

about 

^x6kj 


Dept, of aerospace engineering ' 

INDIAN INSTITUTE OF TECHNOLOGY 
KANPUR-208016 



Chapter 4 


KINEMATICS 


During operation, the rotor blade undergoes deformation in in-plane bending (lag), 
out-of-plane bending (flap), torsion and axial modes. In addition, the hub centre 
has both translational (i?.,,, Ry, R^) and rotational {Ox, Oy, 0^) perturbational motion. 
The formulation of inertia operator and aerodynamic operator requires a proper 
discription of kinematics of the blade motion. In this section, an expression for the 
absolute velocity vector at any arbitrary point ‘p’ on the blade is derived. 


4.1 POSITION VECTOR OF A POINT 

The position vector of any arbitrary point ‘p’ in the finite element of the blade 
(Fig. 9) with respect to the hub centre Ojj, is given by 


7^, = aCyX -f- CiCx-Zk + 626x3* + {< Tl — 1 > Ig + Xk + Uk)eTAk + ^kCyAk 

+Wkez^k + + Ce^ (4-1) 

Transforming all the unit vectors to the 4k - system and neglecting the higher order 
terms, the position vector can be written as: 
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rp = CAj[-a{ft,, cos 0t - ftd sin 0/) + + C 2 + (< n - 1 > /„ + Xk + Uk) 

+^{~•4cos(^?G + (l)k) - 'u4sin((?G + <i>k)} + C{'yfcSin(^G + (t>k) - w'k cos{6g + <?^fc)}] 
+e 4 yl[a cos 61 + ei (/?., cos 0/ - Pd sin 9j - PpSi) + e 2 {Ps cos 0/ - Pd sin 9i) 4- Vk 
+77cos(0g + (i>k) - Csin(^G + (pk)] 

H-eJ^^f-asin 9i - ei{Pd cos 9] + Ps sin 0/ + Pp cos 9i) - e 2 {Pd cos 9i + ps sin Or) + Wk 
+ri{—w'kVk cos{9g + 4>k) + sin(^G + 4>k)} + C{''^'k'^'k + cos(9g + (pk)}] (4.2) 


In Equation 4.2 all the length quantities are nondimensionalized with respect to 
the length of the blade ‘1’. For the sake of convenience, the nondimensional length 
quantities arc referred without an overbar ( ). E(iuation 4.2 can also be written in 
symbolic form as 


fp — ^['^x^xdk “b ‘^y^ydk "b T^C^/ik] 

Differentiating the position vector fp with respect to time 

fp = eiJillaOilPdCOsOj + PsSmOi) + Uk + r]{-VkCos(9G + ^k) 

+v'k<pk sin(0G + (Pk) - ui sin(^G + (pk) - w'k^k cos(0g + <Pk)} 

+C{Ufc sin(f?G + (pk) + v'k^k cos( 0 G + 0 jfc) 

-w'kCosiOa + 4>k) + Wk^kSm{9G + (pk)}] 

+edylQ[—a9r sin Oj — 9r{ei + e 2 ){Ps sin Or + Pd cos Or) — e^PpOr cos Or 

-^k{lsi\\{0G + + C cosiPo + h)) + yfc] 

+edzl^1[—a0r cos ^7 + Or (ei + e 2 ){Pd sin Or — Ps cos Or) + ciPpOr sin Or + tbk 
+{(^fc?7Cos(6'G + <Afc) - Csin(<?G + '5^*)}] 

(4.3) 

where (') indicates differentiation with respect to the nondimensional time 'tp, 


Ip — Q,t 
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Note that 


= = ) = «(■) 


4.2 ANGULAR VELOCITY VECTOR 

The angular velocity vector a) of blade consists of three components. They are 

(i) the constant rotational speed (fi) of the rotor; 

(ii) the rigid body angular velocity (cOrigid) of the hub due to perturbational rotation 
in roll-pitch-yaw (dx, Qy and 6^) 

(iii) the angular velocity contribution due to the rate of change of control pitch input 
0,0] to the blade 

Co = He zu + Corigid + HOje-xzk (4-4) 

where 

^rigid — H0x^x\k ~b HOyCyik "h HOz^zXk (4*b) 

Transforming all the unit vectors of Equation 4.4 to the 4k -system and neglecting 
the higher order terms, the angular velocity of blade can be written in symbolic 
form as: 

CO — H^COx^xAk ~h ^y^ydk T ^z^zAkl (4-6) 

where 

= [Oxoosipk + OySin'tpk + Pp + PdCosOi + PsSinOi + {1 + Pd < sinOi + cosOi > 
+Ps < sin 9i — cos 9/ >)9d] (4.7) 

ujy = [(— 9x sin t/;k + 9y cos ‘ipk) cos 9r + (l + 9z) sin 9/ 

-1-(1 - Pd sin 9/ + Ps cos 9/)9f] (4.8) 

= [(1 + 9z) cos 9/ + (1- Pd cos 9i - Ps sin 9i)9i] (4.9) 
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4.3 VELOCITY AT POINT ‘p’ 

— ♦ 

The absolute velocity vector V , at point ‘p’ on the deformed beam can be written 


as; 

V = Vii fp {oj X f) (4. 10) 

where Vjj is the rigid body perturbational translation of the hub centre Oh which 
is given as: 

Vh = 1[Rx^xR + iiy^yR + Rz^zl^ (4- 11) 

Transforming all the unit vectors in 4k -system 

yjj — iyu)x&x‘\k + {yu)y^y^k + (Kfl) (4-12) 

where 


{Vn)x = f3l//?,,.^{cosVifc - -+/5,jsin 0/)} 

+0,lRy{sinipk + cos-ipki—Ps cos9j-h p^sinOi)} 

+0,lRz{Pp + PdCOsOj + PsS'mSi} (4-13) 

{VH)y = ^ lRx{cos ipkiPs cos (9/ - /^rfsin 0/) - sin^/^jtcos 0j} 

-ffi IRy {sin ipk {Ps cos 9 1 — Pa sin 9j) -+ cos iik cos 0/ } 
+f^ii?zSin0/ (4.14) 

{Vh)z = -^lRx{cos'tpkiPp + Pd cos9i + p^sin9i) + sin iJkSi} 

— Q IRy {sin Ip kiPp + Pd cos9j -H PsSinBj) - cos ■^fcsi} 

-i-n in, cos 0, (4.15) 


and 


^xAk ^zAk 

{lO X f) = UJy CJ, 

Tx Ty r. 


(4.16) 


Where Vy, and Ux, ujy, are the x, y and z components of fp and uj respectively. 


Substituting various quantities in Equation 4.10 from Equations 4.2, 4.3, 4. 7-4.9 and 
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4.12-4.16, the velocity at point p is obtained. This can be written in symbolic form 
as: 

V = n + VyGy4k + Vz^zik] (4-17) 

The components of the velocity vector are given below. 
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Vx — "1“ (^7)):r H“ 

= Rx{cos ipk - sin i>k{-l3s cos Oj + f3d sin 9t)} 

+Ry {sin ipk -\- cos 15s cos 6 1 + PdSin 6 j)} 

+Rz{l3p + Pd cos 0 ! + ps sin 9 ] } 

+a9j{Pd cos 9t + Ps sin 9]) +Uk + v{-K cos{9g + pk) 

+?>fc 0 fcsin (^?(7 + (j^k) - lii'k sin (9g + <Afc) - -x’kPk cos{9a + pk)} 

+C{nfc sin(0(j + (j)k) + v'kpk cos(9g + ^k) 

-rb{ cos(9g + pk) + w{.pk sin((?o + pk)} 

+9x cos9/lsm'ipk{-Wk + rjsin(9G + pk) - Ccos(6'g + <?!'*)}] 
-^yCos(9;[cosT/)jt{-n;fc + 77sin(6>(7 + pk) -(cos(9g + (f>k)} 

-9;;Sm9!{-Wk + r7sin(i9G + (f>k) - Ccos(9g + ij^fc)}] 

-02 cos 0j [a cos + Uk + T] cos{9g + (l>k) + ( sin(0G + pk)] 

-9i[{-Wk + r)sin{9G + pk) -Ccos{9G + pk) + + 62 } [Ps cos 9 1 - pdsm9i} 

+eiPp sin 9i + {cj + 62 - a cos 9/ - Vk - rjcos{9G + pk) - Csin( 0 G + Pk)} 

[Pd cos 9 1 + ^^sin 6 >/} + 77 {cos( 6 Ig + pk) - sin{9G + pk)} 

+(^{sin((?G + pk) + cos(9g + pk) } + 4- Wk + ci(cos 0 / — sin 9i)\ 

- cos 9i[a cos 9j + (ei + e 2 ){Ps cos 9i - Pd sin 0/) + Vk 
+ijcos{9(; + Pk) + C sin(0f? + ^Afc)] 

-sin 0 /[(ei + 62 ) {Pd cos 9 1 - Ps sin 9 1 ) + Wk 
-r) sin( 0 G + ?!>*)+ C cos( 0 g + pk)] 


( 4 . 18 ) 
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Vy — (^H)y ”1“ iXp)y ^z^i) 

= Rx[cos -ipkiPs COS 6 j - Pd. sin Oj) - sin ■0A: cos 0/] 

+i?j/[sin pik {Ps cos 9i - Pd sin 9j) + cos xjjk cos 9i\ 

—Rz sin 9i — a9j sin 9i — ei9jPp cos 9j 
-{9i{e\ +e2)(/0sSin6>/ + /3rfCos6>/)} 

-<^ifc{77Sin(^G + ^^A:) + C cos(6»g + + ilk 

+4(cos’0fc{-‘UJfc + •nsm{9G + pk) - Ccos(6>g + pk)}] 

-0,j[sin 'pk{-Wk + r]sm{9G + Pk) - C cos{9g + ^fc)}] 

-9z cos 9i[ei + 62+ <n-l> le + Xk+Uk- v',^{r] cos{9g + pk) - Csin(6'G + Pk)} 
-tUk{vsm{9G + Pk) -C cos{9g + Pk)}] 

— 0/[{M(sin07 4- cos9j) + Ps{sm9j — cos9[) + 1}{ioa; — 77sin(^G + Pk) + C cos{9g + pk)} 
-l-(A/cos(?/ + /?., sin 0,){2{c] -I- 02)4- < ?/. - 1 > I,. 4- 4- Uk - v4(7/c.os < Oa 4- Pk > 

-C sin <9G + Pk >)} - w'kiv sin < > 4-C cos < 9g + Pk >)} 

-74(77 cos(^G + Pk) - Csin(0G + Pk)} - 7^(77 sin(0G + Pk) + Ccos(0g + Pk)} 

-'^^Wkiv cos{9g 4- Pk) - Csin(6>G 4- pk)} + a{PsCos9i - Pdsm9i) 

— (c] 4- 62 )— <n— l>le — Xk — Uk] 

-[(/?p 4- Pd cos 9 1 4- Ps sin 9i){wk-'n sin(6'G 4- 4- C cos(6’g + Pk)} 

4- cos 9i{-Vk{ri cos{9g + Pk) - C sin(6>G 4- pk)} 

-w'k.{vsin{9G + Pk) + ( cos{9g 4- <?ifc)} 

-a{Ps cos/?/ - /?rfsin 0,) + Ci 4- C24- < n ~ I > + Xk 4- 77,*;}] 


(4.19) 
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Vz — “I” (j" p) z “!■ (j^x'l’y ^y'^x) 

= -i?x[cos ipkiPp + Pd COS 6 1 + Ps sin 6/) + sin ipk sin dj] 

-Ry[sm i^kiPp + Pd cos d[ + Ps sin Op) - cos ipk sin 9i] 

+RzCosdi — di{acosdi — eiPpSinOj) 

+0i{ei + 62 ) {Pd sin Op - P^ cos Op) + {pkO cos{0g + Pk) - C sin(^G + ^it)} + Wk 
+4[cos '0fc{a cos Op-^-Vk-^r] cos(6'g + Pk) + <^ sin(6'G + (l>k)} + sin V-'fc cos Op {ci + 62 
+ <n-l>le + Xk + Uk- v'f.{T] cos < Og + Pk > -C sin < Og + pk >) 

-w[{r} sin < Og + Pk > +Ccos < Og + pk >)} 

+ 0 ,;[sin pk{a cos Op + Vk + ri cos( 6 »g + Pk) + C sin(0G + Pk)] - cos V-’jt cos Op {e, + 62 

+ < n - 1 > le + Xk + Uk - v[{r) cos < Og + pk > -C sin < Og + pk >) 

-w'k{r] sin < Og + Pk > +C cos < Og + Pk >)} 

+0 p[{bd{sin 0 p + cosOp) + Ps{sinOp — cosOp) + l}{acos07 + Vk -\-i]cos{0g + pk) 
+Csin(0G + Pk)} - eiPp sin Op 

-{Ps cos Op — PdSinOp + l){<n-l> Ig + Xk + Uk- v'k{r]cos < Og + pk > 

-( sin < Og + Pk >) - yj'k{v sin < Og + pk > +C cos < Og + pk >)} 

+n(/?,5 cos Op — Pfi sin Op) + (ci + C 2 )] 

+[{Pp + Pd cos Op + Ps sin Op) [a + 77 cos ( 6 >g + pk) + C sin(i9G + Pk)}] 

+(1 + Oz) sin Op [ei + C 2 + < n - 1 > + Xk + Uk 

-v'k{r} COS < 6'g + > -C sin <0G-\-pk>) 

-w'k{r]sin < Og + pk > +C cos < Og + pk >) - aOp{Ps cos Op — Pd sin (4.20) 



Chapter 5 


EQUATIONS OF MOTION FOR 
ROTOR BLADE 


The equations of motion for the rotor blade can be derived using Hamilton’s prin- 
ciple. The mathematical form of Hamilton’s principle is stated as follows: 

f\5U -6T-5We)dt = Q (5.1) 

where U is the strain energy; T is kinetic energy; Wg is the work done by the non- 
conservative external loads. Equation 4.1 is an integral equation which states that 
the total dynamic potential, {U — T — W^), is an extremum over the time interval ; 
ti< t < t2. 

In this chapter, the expressions for the variation of kinetic energy and strain energy 
of the rotor blade are derived. 


5.1 KINETIC ENERGY OF BLADE 

The total kinetic energy of the beam, T, is defined as : 

= IJ‘ j f^pV.Vdr,d(dx (6.2) 
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where V is the velocity of an arbitrary point ‘p’ on the blade cross section with 
respect to the inertial reference frame. The variation of kinetic energy, dT can be 
written as: 

5T = m l^pV.SVdrjdCdx (5.3) 

Substituting for the velocity V from Equation 4.17 and integrating dT by parts with 
respect to time, yields 

S T = J j j p\Zy, 6u + Zy Sv + Zyj Sw + Z'^ Sv' + Z'^ 6w' + Z^f, 5(j)]dT] d(^ dx (5.4) 

where the terms Zy, Zy,, Z', Z^ are the coefficients of 6u, 5v, Sw, Sv', Sw', 
S(j) in the variation of kinetic energy expression. 

Integration of the expression over the cross-section yields: 

ST = C [Zy Su + Zy Sv -k Sw + Z^ Sv' + Z' Sw' + Za, S(j)]dx (5.5) 

Jo 

After eliminating the higher order terms, the expressions for Zy, Zy, Zy,, Z'y, Z'y,, Zf^ 
are obtained. The detailed expressions of these coefficients are given in Appendix- A. 
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5.2 STRAIN ENERGY OF BLADE 

The formulation of strain energy outlined in this section essentially follows the 
procedure given in Refs. 23 and 28. 


5.2.1 Strain Energy 


The strain energy of a beam element is given by 


U 




T 

r N 


^xx 


Oxx 


^ ^XTJ 

> < 

^ XT] 

'> dT]d(dx 

^xC 

\ ^ > 


^xC, 

\ ^ J 



5.2.2 Explicit Strain-Displacement Relations 


(5.6) 


The expressions for non zero strain components written in terms of u,v,w and cj) are 
given as [23] 

-[ 77 COs(^p + <^) - Csin(^p + <^)]uix - [rjsm{$g + (p) + C cos( 6 ', + (f>)]wxx 
T^('Tit 7 ,i ’T'o^xd^) T CCTiC.i "h '^oIxT]) 


Cijj — “b ri C(*^i d" 4>o) 


^xC, — 7 xC + Vi^x + (po) 

where 

(j)g = [vxx cos 9g + Wxx sin 9g){—Vx sin 9g + Wx cos 9g) 

The underlined term in Cn represents the axial strain at the elastic axis. 
These strain expressions can be simplified using the following assumptions: 


• The transverse shear at the elastic axis is assumed to be zero. 
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• The warping amplitude a is assumed to be equal to —<px- 
The simplified strain components can be written as 

e.TX = + 2^'’f + C^)^x ~ 

-[i] CQs{9g + <p) - C sin((9p + (j))]vxx - [r]sin{6g + (^) + C cos{6g + (p)]wxx 
^xri = —{'^n~^0^x — C4>o 

exc = -('J'c + 
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-{•q sin 9g + C cos 9g){5wx:, - (j)6vxx - Vxx^<f>) 

5exri = r, 0^4>X - <,5(i)o 

Sexc = -{^^-q)5(px + vS(l>o 
where 


5(j>o = {—Svx sin 6g + 5wx cos dg) (vxx cos 6g + Wxx sin 9g) 

+{-Vx sin 9g + Wx cos 9g){6vxx cos 9g + 6wxx sin 9g) 

It is assumed that the variations of strain components are of the same order as 
the corresponding strain components. Substituting the above expressions in the 
strain energy variation and integrating over the cross-section, the variation in strain 
energy for a beam element is obtained. The expression for 5U, in terms of stress 
and moment resultants, is given as: 

rli _ 

5U = Eol^ {Vx{5ux + Vx5vx + WxSwx) 

Jo 

TqPx Rx ”b Sx d” '^x^x^^^x “b Rx Rx^'i^xx 

+[vxx{Mrj' COS 9g — sin 9g) -f WxxiMg sin 9g + M(^ cos 9g)\5(f) 
sin 9g -i- Mq cos 9g) + <i>{Mri cos 9g — M(^ sin 9g)]5vxx 
+[{—Mg cos 9g -f sin 9g) -f (f>{Mg sin 9g + Mc^ cos 9g)]6wxx}dx 

(5.9) 


Stress and Moment Resultants are defined as: 

K, = J j{Ecxx)di](K 

= EA[ux + ^vl + ^wl] 

+^EACo[(j>l] - EADo[(j>xx] - ToEAD'M 
— {EAria - <^£ACa)[Uix] - {EACa + <i’l^A-na,)[Wj^x] 


(5.10) 
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Mr, = J J {CEe^^)dr]dC 

= EA(^a[ux + 2^1 + 2^®^ 

+^EAC2[<i)l] - EAD2 [<Px.] - ToEAD'M 

(j)EIfiTj)[Vxx] i^Iirr] “I” (5.11) 

M^ = J j {CEe^x)dr]d( 

= EACa[Ux + -vl + -tul] 

—EAD2[<i)x:^ ~ '^oEAD'2[4>x] 

-(WdM-iM^^nwxx] ( 5 . 12 ) 

Me; = I J {-7)Eexx)d7idC 

= -EAr]a[ux + ^vl + ^wl] 

-^EACi[<l>l] + EADi[(I>xx] + ToEAD'M 

+(ET^C - 07rv)M + iWv + 0kc)M ( 5 - 13 ) 

M'^ = J J {-r)E€xx)dr]dC 

= -EAria[ux + ^vl + ^wl] 

-{-EADi[(f)xx] E ToEAD'i[(f)x] 

+ + (ET(^^)[Wxx] (5.14) 

Px = j J {'i>Eexx)di]dC 

= EADo[ux + ^vl + ^wl] 

+^EAD^[<j>l] - EAD3 [(I>xx] - roEAD5[(l>x] 
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-{EADi - 4>EAD2)M - {EAD 2 + <f>EAD,)[w^^] (5.15) 

K = j /[(C^, - 

= EAD'J^U:^ + -vl + -Wj] 

+]^EAD'Ml] - EAD,[<j^,,\ - ToEAD'M 

-(bad;' - cfEAD2)[v,,] - (Em' + cfEADl'm,] (5.16) 

[( 77 ^ + (^)Ef.^,:r]d7]dC 
= EACo[u^e]^vI + ]^wI] 

+]^EAC,[cl>l] - EADm] - ToEAD'm 

-(EAC( - (fEAC2m] - (EAC 2 + (^EI^')K^] (5.17) 

— y J [(^ r^€xrj ^ (^^x^)E!]dTjd( 

= iGJi-GJ2)[(l>x]+GM<l>o] (5.18) 

Sx = y J C^xri)G]d7]d(^ 

= {GJo-GJ,)[<Px]+GJo[<Po] (5.19) 

Sx' = y y[(77exc - Cexr,)G]dridC 

= {GJo-GJ,)[(l>x] (5.20) 

In the above expressions, the moment resultants and Sj have the same 

definitions as Mrj, and Sx respectively. But higher order terms are neglected in 
defining Mrj' , M(( and Sj ■, since these expressions are coupled with terms of order 
e^. Whereas and Sx are coupled with terms of order e. The cross-sectional 



integrals associated with the strain energy variation are defined as: 

EA = I J[E]d7]dC 

EAria = j j [E'!i\dr)d(, 

EACa = I JmdvdC 

Elr,, = j J[Ee]dvdC 

EI^^ = J l[Er,^]drjd(: 

EI,,^ = I l[ErjC]dr]dC 
EACo = I JiEiv^ + OldvdC 
EACr = I l[Ev{v'^ + e)]dvdC 
EAC2 = I l[Earf+e)]dvdC 
EACz = j l[E{v'^ + Cy]dridC 
EADo = 11 [E^dr]dC 
EADi = j I [Eri'^]dridC 
EAD 2 = j JlEC^ldrjdC 
EADs = I j[E'^^]d7]dC 
EAD, = J j[E{ij^ + e)^]drjd(: 

EAD, = I |[ii;v]/(Cd/,-r/'I/c)]d7?dC 
EAD'„ = I lm%-rj<if^)]dvdC 
EAD[ = I l[Ev{C%-V^d]dvdC 
EAD '2 = j j[EC{C<i!^-'n'^,:)]drjdC 
EAD'^ = I J[E{C%-v'^^f]dr]dC 
EAD\ = j j[E{rf + (:^){(:^r,-ri^d]dyidQ 
GJo = I j[G{v^+e)]dndC 
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GJ, = J l[GK^„-r,<Slt)]d,,di 
G.h = J /[G($; + 4<?))c!r,rfC 

Certain additional terms are defined as: 


EArja 

zzz 

EArja cos 9g — EAQa sin 6g 

EACa 

__ 

EArja sin 6g + cos 6 g 

EAC, 

m 

EACi cos 9g — EAC 2 sin 9g 

EAC 2 

= 

EACi sin 9g + EAC 2 cos 9g 

EADi 


EADi cos 9g — EAD 2 sin 9g 

EAD 2 

zzz 

EADi sin 9g + EAD 2 cos 9g 

EADi 

= 

EAD\ cos 9g — EAD 2 sin 9g 

EAD 2 

nz 

EAD'i sin 9g + EAD'^ cos 9g 

EItj^ 

= 

ElgQ cos 9g — Elgg Sm 9g 

EIr,n 

= 

Ek^ sin 9g + Elrjn cos 9g 

Eke 

= 

Eke ^0® ^9 ~ Ekv ^9 

Ekv 

= 

Eke ^9 + Ekv ^9 

GJ 


GJo - 2GJi + GJ 2 



Chapter 6 


FORMULATION OF ELEMENT 
MATRICES ASSOCIATED 
WITH KINETIC AND STRAIN 
ENERGY VARIATION 

6.1 FINITE ELEMENT DISCRETIZATION 

Tho variational expressions associated with the kinetic and potential energy of the 
rotor blades are nonlinear. The unknowns are the deformation functions u^v,w, 
and (j). These are dependent on both space and time. The spatial dependence is 
eliminated using a Rayleigh Ritz finite element formulation. The blade is divided 
into subregions(finite elements) as shown in Fig. 9, and total dynamic potential is 
calculated for each subregion. By applying Hamilton’s principle to each subregion, a 
discretized form of the equations of motion can be obtained. In this developement, 
each subregion is modelled by a straight beam type finite element. These beam 
elements are located along the elastic axis (line of shear centers) of the blade. 


31 



32 


The (iiscretized form of Hamilton’s principle is written as: 

rt 2 j!L 

^{SUi-STi- 5 W,i)dt = 0 (6.1) 

i=i 

where: 

N = Total inimber of finite elements in the model. 

5 Ui = Variation of the strain energy in the element. 

6 Ti = Variation of the kinetic energy in the element. 

S Wei = Virtual work of the external loads in element. 

Assume that the four unknown displacement functions of the beam finite element 
are expressed in the following form. 


?; 
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0 
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. {U} . 


where , {$,4} , {$„} are space dependent interpolation functions; and 

{V} , {W} , {$} , {t/} are time dependent nodal values of v, w, (j), u, respectively. 
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The nodal coordinates are shown in Fig.lO. 

Tli(! variation of the displacement functions for tlie beam can bo written as: 


5 V 


I 

0 

0 

0 

j 


■ {fy} ■ 

5 w 


0 

0 

0 


{ 6 W) 

64 > 


0 0 0 



5 u 


. 0 0 0 _ 


, _ 


( 6 . 3 ) 


In this development, Hermite interpolation polynomials are used for the space 
dependent interpolation functions. A cubic Hermite interpolation polynomial, {$c}) 
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is used for the bending deflections (v,w) and a quadratic Hermite interpolation 
polynomial,! $ 5 }, is used for the torsional rotation (0), and the axial deflection (u). 
The mathematical expressions for these polynomials are given as: 


{d>.} = {^K} = 


1 _ 3^2 + 2^3 

ie{^-2e+e) 
ze - 2^^ 
ic{-e+e) 


{'I’c} 


(G.4) 


where: 


{$4 = {^u] = 


l-Z^ + 2e 

4^ - 4^2 

-e+2e" 


{<!>,} 


(6.5) 


f = - 

X = spanwisc (axial) coordinate of the beam element. 

Ig = length of beam element. 

For bending deformations, the nodal parameters are the displacements and slopes 
at both ends of the beam element. Therefore, the resulting elements will have 
interelement continuity for both displacements and slopes. In addition, because 
of the cubic interpolation polynomial, the bending strains vary linearly over the 
clement length. 

The quadratic Hermite interpolation polynomial is used for the torsional rotation 
{(f)) and the axial deflection This polynomial has the capability of modelling a 
linear variation of strains along (,he element length and therid'ore provides ( he sainf! 
level of accuracy as the beam bending element. The nodal parameters for these 
elements are chosen as the values of the displacement function at the two end nodes 
and at the mid-point of the element. 

The resulting beam element has 14 degrees of freedom: 4 in-plane (lag) bending 
degrees of freedom, 4 out-plane (flap) degrees of freedom, and 3 degrees of freedom 
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each for torsion {(j)), and axial deflection (u). The nodal degrees of freedom are 
shown in Fig. 10. 


6.2 ELEMENT MATRICES ASSOCIATED WITH 
KINETIC ENERGY VARIATION 

The beam element matrices associated with the kinetic energy variation are derived 
by substituting the assumed expressions for the displacement functions (Equation 
6.2) in the kinetic energy variation 6T (Equation 5.5) and carrying out the inte- 
gration over the length of the beam element. The resulting variation of the kinetic 
energy can be written in the form; 




14x1 



f *• 

Rx 


^ . N 

Rx 


H 

V 


f . N 

+ \m^] 

L Jl4x3 

iky 

L Jl4x3 

Ry 

> + < 
L Jl4x3 


> + [m^I 

L J 14x3 

dy 


Rz 

x J 


Rz 

\ J 


s, 

V J 


0z 

v. y 




14x1 


( 6 . 0 ) 


Where {q} represents the vector of unknown nodal degrees of freedom 
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[V] 

m 

{cl>} 


> 


i {t/} 


(6.7) 


Detailed expressions for the various matrices defined in Equation 6.6 are given as 
follows: 



3G 


6.2.1 Mass Matrix 

[Mil] = / H {^c} + {JiTT'CC cos^ <0G + <l>k> +Imr,r, sin^ <0a + <pk> 

J 0 

-2Imr,^ sin(0G + 4 >k) cos(0g + <f>k)} {$1} dx 

'[M12] = sin(6'G + <i>k) cos(6'g + ^k) + /m^^(cos^ < 0 G + (i>k> 

-s-m^ <0a + <l^k >)} {<} {<I>',}'^] dx 

[Mja] = ^ + (/>fc) - ?>iCm i^oii{ 0 a + </»*:)) (I’c) ] dx 

[Mu] = [{-mr)m cos( 0 G + <i>k) + mCm sm{ 0 n + ^>k)] {I’l} {^7}^] dx 

/’^e 

[M 21 ] = / - Im,jn) sin(<?G + </>*=) cos{0g 4- (l>k) 

J 0 

+/m^(;(cos^ < 0G + (j)k > - sin^ < 0g + (f>k >)} {^1}^] dx 

[M22] = [ \m {<&c} {^c}^ + {dm^c sin^ < 0 G + (pk> +I'rnr,r, cos^ < 0 g + (l>k > 

Jo L 

+ 2 Imr,i sin(6>G + (l>k) cos{0g + (l>k)} {<I^} {C}^] dx 
[M23] = [{mr)m cos{ 0 G 4 - <j)k) - m(rn sin( 0 G 4 - <j>k)} {^c} dx 

[M24] = sm(t>G 4 - (j>k) 4 - rnCm cos{0g 4 - (l^k)} ] dx 

[M 3 ,] = [M, 3 ]"’ 


[M 32 ] = [M23]'^ 
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[Mas] 

-i:\ 

, {$,} {# dx 

[Ma,] 

= [ 0 ] 


[M,,] 

= [Mu 


[M42] 

= [M24 


[M43] 

= [ 0 ] 


[M44] 

= /' 
Jo ' 

[7/1 {< 1 ),,} {<!>,}''] dx 
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6.2.2 Matrix 

[Mfi] = ^ cos(0 G + (t>k) + mCm sin(^G + </>*)} {^>'c} {^c}^] dx 
W 12 ] = / [{m(/3p + /?dCos07 + /3,sin6>i)}{$c}{'J'c}^] cia; 

[M 13 ] = / [-27m^^ sin^ < 6 'g + (/>* > {$ 9 }^ 

J 0 

H-2(int^^ + <;os Oi sm{Oc + tj)k) co3{0o + ipk) {'t,} ’ ] dx 

IMfJ = [01 

fh 

[M|i] = / [{m77^ cos <60 + (j>k> —'mCm sin < 0 a + <i>k> -rn-Vm sin < 0 a ^ 4>k> 

^ 0 

rjp 

-rnCmCOS <0a + (j’k >}(/?,, cos <?f - Pd^mOj 4- 1) {(T>^} |<T>'J 
-{imjrn sin(( 9(7 + (pk) + viC„iCoa{0a + (pk)} {‘I’l} 

+{mr]mSin{eo + (pk) - mCm cos( 0 g + <Pk)} sin0/ {'i’c} {^ 1 }^] dx 

[M 22 ] = / [{mr]m sin(0G + <l^k) + mCm oos{9g + (t>k)}{Ps cos 0/ 

Jo 

= -/?dSin Oi + 1) sin Oj {$c} dx 

[MI 3 ] = jj [{-2mrijnSin{6G + (pk){Pp + PdCOsOi + PsSinO])} {^c} {^qV 

-{mrim sin( 0 G + <Pk) + m-Cm cos( 0 g + (pk)} {'^'c} {^<7}^ 

-a{rm]m cos(/?g + (pk) - rnCmSin(6»G + (pk)} {^^1} 

-2/m^^sin^ < Oq + (pk > {^’c} 

sinOi - 2 /m^^}sin( 0 G + (pk) cos{9g + <pk) {^c} {^ 9 }^] 


[Ml,] = [0] 
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[Mil] = cc + sm{Oa + <f>k) cos(0G 4- (pk) - ‘2Imr,,; cos^(0g + </»fc) 

+2(/m^^ + /m,^) sin(0G + (/)fc) cos(0g + ({>k) sin 0/ 

-{hn(^(^ - i7n^,;) sin(OG + ^!»fc) cos((?g + ^k) 

-2/m,,^cos^ <0g + <j>k >}cos0/ {$,} 

{MI2] = f cos^((?G + - Irur^ + Irrir^ sin(6>G + ^fc) cos(6'g + <?i)k) 

J 0 

-/m^^(2sin < /?g + </>fc > cos < Oq + (j)k >) 
cos^ < 0G + ^k > 

+2(Im(i; + Imrjrj) sin(^?G + M cos(<?g + (/^k) sin Oj 
+mr]mSm(dG + (t>k){0p + Pd cos 9 1 + PgSinOi)} {$,} dx 

[Mlal = / [[(mr?, „ cos(0G + Pk) - mCmSin(6>G + (j^k)} {Pd cos 0, + p, sin 0j){< n - 1 > h ■+ Xk) 

./o 

- siii((?G + ^k) - 'rnCm cos(6>g + (pk)} 

{a - {Ps cos 9i - pd sin ^/ + 1)(< n - 1 > Ze + ^k) 

-(ei + 62)} - (/m^< + Irurirj) sin(<?G + (pk) cos( 9 g + MiPp + PdCOsOj + /?« sin0;) 
-{Im(^(^ + /m^^)sin(0G + (jf»fc)cos(^G + (i)k){Pp + Pd cos Or + Ps sin Oj) 

-]m,,(^{3siiP{0c + (t>k) - 5cos\0a + (l>k))iPp + PdCosO/ + sin Oj) {<!>,} {<!>,}'] rii- 

Hr 

l-Ai/IJ = / [- 7 /;, cos 0, + { 7 / 77 /,,,, cosi^^G -I- - HiCni sm{()o + (l>k) sin 0i 

Jo 

-{mrjjn sin(0G + (j)k) - mC,m cos(0g + Pk)} cos 9 r {^9}^] dx 
Wti] — [tticos^/ {$ c}^] dx 
[^ 42 ] = ^ [7nsin6»/{$J{$c}^] dx 
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[M^j] = ^ [{-2m77„ sin(^G + ^k)} {$ J dx 

[M^^] ^ [{m, 7 ]m sm( 0 G + <Afc) - ”^Cm cos(6>g + <^<;)} cos Or {$ J {$7}^ 

-{{m'nmCOs{ 9 G + </>fc) - mCmCOs(6>G + ^6fc)}sin^/ {#^} {$7}^] dx 
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6.2.3 Matrix 

[/CiY] = ^ [-mcos<?/ {$J {<I)J^ 

+{/m^^ cos^ < 9G + ^k> +Imrp, sm"^ < Oq + (j)k >} cosOj 
-2/m^^ sin(0G + (f>k) cos{Og + (f)k) cos6i dx 

[d<n] = J ^‘{[- + (l>k) COs(/?f7 + .J&fc) 

-Jm^^{cos^(6>G + (i>k) - sin2(6>G + <i>k)} 

+rn?/„i cos{Og + (i>k) - mCn sm{OG + d>k)] cos 0, {4>c} 
+msin 6i cosdj {<I»c} dx 

[^fl] = M 

[Ktl] = lo] 

[/f|f] = j l{{I”’-a-^’'‘’m)’>"'(0o + <l>k)'xs{0a + <l>h) 

+Iui,,^{sh{^{0a + (l>k) - cos^(0g + <l)k)]]cosOi |<t^} {'^1}^ 
+m sin di {$c} {*^’c}^] dx 

[k^I] = f [-{mr]m cos{9g + (f>k) - mCm sm{9G + (j>k)} 

(p, cos 9i - Pd sin 9i 4- 1) {‘5c} 

sin^(0G + pk) + drv.r,r, cos^((?g + pk) 

+2/m,,f sm({?G + Pk) cos(9g + {<5'^} {^5^}^ 

4-m.sin0/ {‘^’c} {^c}^] dx 
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[K2i] = [0] 

[K^2(\ = [0] 

[Kfi] = ^ ' [[{mrjm sin(6'G + - mCm cos(6>g 4- + {/rn^c sin(6'G + 0fc) 

-Jm^^cos((9G + <?i>fc)}(; 0 p + PdCosOr +/?,sin0/)]cos(9/ {$,} {$c}^ 

+ ImTjr,){Pp + Pd cos df + Ps sin 9i) 

-{mrjm sin(0G + (pk) + ^Cm cos((9g 4- ^fc)} 

(f!, -I- C 2 -\- < n - 1 > /« + -Xk) {4>„} d:r. 

[kII] = ° - Im(^d sin(6'G 4- ^k) cos{9c 4- Pk){Pp + Pd cos Oj + ps sin 9r) 

+{mr]mCos{eG + <f>k) - mCmSin(^G 4- <?^fc)}(ei + 62 + < n - 1 > IgXk) 
+2/m^^{sin^ < 6g + (l>k > — cos^ < 9g + Pk >} 
iPp + Pd cos 9i + Ps sin 9i] cos 0/ {$,} dx 

Kli\ = [ 0 ] 

K1 = [0] 

[kII] = jj j-{m77„iCos(0G + Pk) - mCmSin(6>G 4- <?!>fc)}cos6>/ {$,} {^>c}^] dx 

\Kli\ = ^ ' \rn{Pp + Pd cos 9i + Ps sin 9i ) cos 9i {<!>,} { ^ j dx 

[KaI] = [0] 

[-7m(cos/?/ 4- sin^ (?/ {4>,} {<[> 5 }'] fix 
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6.2.4 Vector 

["^u] = ^ ' [-{^Vm sin(0G + (f>k) + m(m COS(0G + 0fc)} cos^ Oj {$c} 

{ma{2 cos 6 i — sin 6 i) cos 9 i — m(ei + e2)(/?j cos 6 j — Pd sin 6 i) cos^ 9 i} {$c} 
-{mr 7 „jSin( 0 G + h) + mCmCos( 6 »G + <i>k)} {"^c} 

+{mr]mCos{9a -f- pk) - mCmSin((?G + '?^fc)}(< n - 1 > le + Xk) {^c}] 

[^2^] = ^ [{-n7-(ei + 62+ < n - I > le + Xk){Pp + Pd cos 9 i + p, sin Or)} {‘I’c} 

{ma(2 cos 9 f — sin 9 i) sin 9 t { 4 >c} 

+{mr7mSin(^G + h) + mCmCOs(^G + (l)k)]{< n - 1 > k + Xk) {^c}] dx 

[Vg^] = ^ sin( 0 G + (f>k) cos((?G + Pk)} cos^ 01 

-/m,^{sin^(^G + pk) - cos^(^G + </'jfc)} cos^ 9j 
+{Im^^ sin{ 9 a + pk) + /m,,,cos(0G + 0Jfc)}sin07] {<&,}] c/rc 



fU 

/ [-m{ei + 62+ <n-l>le + Xk + {PdCos9i + /?jSin0/)(< n - 1 > k + Xk) 
J 0 

—a{Ps cos 9 r — Pd sin 9 r)} cos 9 i {$5}] dx 


Using small angle approximation for <f>k, the vector can be written in a modified 
form, which is given in Appendix B. 
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6.2.5 Matrix [M^] 


14x3 


K1 = /, i{ — m sin ipk} cos 6i dx 

12] = j [{mcosipk} cosOj {^c}]dx 


Kl = h —msmOi dx 
K] = rif — m cos tj)k) {i’c}] dx 
[Mj‘2] = jj l{-msmi:i,}{$^}]dx 
[Mj's] = lm{ic]]dx 

K‘i] = [01 


K2] = [0] 

Ks] = [01 

[iW4i] = f [{— mcosV'fc} rfa; 

[ML] = / [{msini/jfc} {$,}] dx 

K1 = i [m{l3p + I3d cos Oj + (3s sin 9i) {$,}] dx 
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6.2.6 Matrix [M^' 


14x3 


Kl = f [— m{2 cos 'ipk + {Ps cos 9r — Pd sic 9i) sin cos di dx 
KI = t [— m{2 sin + {—Ps cos 6i — Pd sin di) cos Tpk} cos 9i {$c}] dx 

l^iz] = ‘ [^i {^c}] dx 

[mIi = J [m{sin tpk (Pp + Pd cos 9i + PsSm.9j ~ 9i) - 9r cos i)k} {^>c}] dx 


Ky - / [ m{cos 'ipk{Pp + pd cos 9i + Ps sin 9i - 9i) — 9 t sin 'tpk} {^c}] dx 


[ 


Ml 


[01 


K] = [01 


K] = [0] 

Ks] = [0! 

[<] = i [m cos ipk{Ps cos 9r — pd sin 9i) cos 9i {$,}] dx 

Ky = r [m{2 cos 'ipk + {Ps cos 9i — Pd sin Oj) sin ^fc} cos 9i {$,}] dx 


Kl = [01 
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6.2.7 


Matrix [M^' 


[ 0 ] 


14x3 


Kj = [01 


K] = [0] 

KI = t [m{a cos tpk + (ei + 62) sin tpk} cos Oj {$c}] 

= r [m{a sin tpk — (ei + 62) cos ipk} cos 9i {$c}] dx 

= [ 0 ] 


ll = 


= [ 0 ] 

K] = [01 

[Mil] = [0] 

[a4] = [0] 

[.V/^j = j [— nwcos0;{$g}](ix 
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6.2.8 Matrix [M^' . 


I *<6 

/ [[—■'^{'2 COS dj — {Ps COS di — Pd sin 0/ + 1)(< n — 1 > Zg + 2 :*) 

Jo 

— (ci + 62 )} cost/)A; + m{(< n — 1 > Zg + ^k) sin^t}] {$c}] dx 


KV 


/ [[— cos — {Ps cos 61 — PdSm9r + l){<n — l > Zg' + Xk) 

Jo 

— (ei + 62 )} sin'^fc — m{(< n — 1 > Ig + Xk) cos^^}] {$c}] dx 


fh 

/ [ma cos 6i{2 cos 9r — sin 9i) {$c} {*^c}] dx 
Jo 


—m{ 2 a sin ipk cos 9i + {< n — 1 > Ig + Xk) cos tpk} {$c}] dx 


Ki = r [m{ 2 a cos ■ 0 A: cos — (< n — 1 > Zg + Xk) sin tpk} {'^c}] dx 

= f ‘ n - 1 > Zg 4- Xfc) {Pp + Pd cos 9i + ps sin 9i - 9i)} {$c}] dx 


[^^"32] 

ij' 4 i 

Ni] 


/ [— m{a cos 9j - {ps cos 9i — pdSiaOr + l){<n - 1 > lg + Xk) 
Jo 

-(ei + e 2 )}smtpk {$,}] dx 


Kl = f [m{a cos 9 1 - (i, cos i 9 / - sin ^/ + 1) (< n - 1 > Zg + Xk) 
— (d + 62)} cos tpk {$,}] dx 

[-''«] = r [7n{(ei + 62) -h 2(< n - 1 > Zg + Xk)} cos dr {-I*,}] dx 
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6.2.9 Vector 

[^u] = [-m[{3a- {ei+ 62 ) - {Ps cos 9i - PdSin9i -rl){< n - 1 > le + Xk) 

—{Pd < cos 9 1 + sin 9 1 > +Ps < sm 9r — cos 9r >)}9i 
-{{Pd cos 9i + Ps sin - 1)(< n - 1 > Ze + ^k)}9i 
4-{a + {Pd sin 9i + ps cos 9r){< n - 1 > 1 ^ + Xk)}9j] {$c}] dx 

[^ 21 ] = [-m[{a + {Pd cos 9 1 + /?* sin 9i){<n-l>le + Xk)}9i 

—{a cos 9i - {Ps cos 9i — PdSm9i + l){<n-l> le + Xk) - (ci 4- 62 )}^/ 
—{Ps sin 9r + Pd cos 9r){< n — l>le + Xk)9]] {$c}] dx 

K] = M 

= J ' [^-m[{2(< n-l>le + Xk){Pd cos 9i + Ps sin (9/)}^/ + a9i 
-{{ei + 62 + <n-l> le+Xk){l + Pd cos 9i + ps sin 9i) 

— (ei + 62 — CL)}9j + a sin 9i cos 9i9i] {$g}] dx 
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6.2.10 Nonlinear Vector 

K^] = /' [-m[2ei cos Oi {$c} {^c}^ {^} + {&] cos 0, + 9:6, - Bj) {<5,} {W} 

+9: J {^,f {U} + 29: {%f {^7}]] dx 

+$}{i,}{i,f{W}]]dx 

K'"] = [0] 

[vil^] = _('‘[ml(-«7){4,}Rf {v} -(«/ + «?) cos 9, {«,}{*„}’■ {W} 

-«/ {«,} {*,}"■ {f^})] dx 
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6.3 ELEMENT MATRICES ASSOCIATED WITH 
STRAIN ENERGY VARIATION 

The element matrices associated with the strain energy variation are derived by 
substituting the assumed expressions for the displacement functions (Equation 6.2), 
in the strain energy variation, 5U (Equation 5.9), and carrying out the integration 
over the length of the beam element. The resulting variation of the strain energy 
can be written in the form: 

6U = Ej {i 5 f { [if®] {q} + {F®}} (6.8) 

Where 

{q} represents the vector of unknown degree of freedom. 

is the linear stiffness matrix. is the nonlinear stiffness vector. 

The detailed expressions for the and are as follows: 

6.3.1 Linear Stiffness Matrix [K^] 

The linear stiffness matrix [K^] is given by the following sub-matrices; 

[irg] = 

[Kl] = KT ■ 
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KJ = /'■ ™ «») {*;} {«:}’■] ‘fa 


K] = l^[{EAD:,)m{i“f+[T,EAD,'){$:}{i',}^dx 


[Kg] = jf'’ [{-BXC4 «} {«;f ] * 


= 'K- 


Ef 

13 




^'•((£, 41 ) 3 ) { 4 ;} {*;'}"■ + (T,EADM{i'i} {*;}"■ + {■»;} { 4 "}^) 

+{r^,BAD',+GJ){i',}{i',f]dx 



/'• \i-EAD,) {4;} { 4 "f + (-t,£AD;) {4;} {4;}’'] dx 

«/ 0 L 


[irfi] = KJ'’ 

[Kg] = [Kgf 

[ifl] = [ifS]"' 


[jf£] = [(s--!) {*'.} Hr] fa 


-tiN. ‘-I3RA81I 

i 1. .-NPUR 

4lil ■ ji ^ ^-Or^ii 
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6.3.2 Simplification of Linear Stiffness Matrix [K^] 

A simplified form of the lineair stiffiiess matrix can be obtained by making the 
following assumptions: 

1. The elastic axis is coincident with the modulus weighted centroid. 

2. The e,j, Cq are chosen parallel to the principal axis. 

3. Warping function is assumed to be zero. 

With these assumptions EA,EIr^,EI/^^,EACo are the only integrals that remain 
non-zero. The linear stiffness matrix can be written as: 



where 


{EIq^ COS^ ^5)[^] 
{-EI^ sin dg cos 05 )[- 4 ] 

[O] 

[O] 


{-EIjjj, sin 9g cos 9g) [A] 
iEImCos^9g)[A] 

[ 0 ] 

[ 0 ] 


[0] 

P] 

p] 

PI 

{GJ)[B] 

P] 

P] 

{EA)[B] 


(6.9) 


[^] 



12 6le -12 6le 

-12 -61e 12 -6k 

6k 2ll -6k 412 



7-8 1 

-8 16 -8 
1 -8 7 
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6.3.3 Nonlinear Stiffness Vector 

The nonlinear stiffness vector [F^] is given by the following sub-vectors: 

^ { [(14) K) - [Sj Sm9g){Va:x COs9g -I- Wa^ajSin^j)] {$(.} 

+[-{ElQg cos 9 g - EIj^sm9g){(j>){vxx) 4- {EIqi; cos 9 g -WT^sm9g){<f)){wxx) 
—~{EAT]a) (vl + wl) + {Sj cos 9g) (-Vx sin 9g + Wx cos 9g) 

Ld 

+{M,'cose, - MiSme,)(4.) - i(S3CD(4)l 

{ [(14) {wx) + {sj cos 9g) {Vxx cos 6g + Wxx sip dg)] {$[.} 

+[- sin 9g -1- Elgr^ cos 9g) {(j)) (vxx) + {EI(^q sin 9g + EIr,<^ cos 9g) ((^) {wxx) 
~^{EAQ) {vl + wl) -1- {Sj sin 9g) {-Vx sin 9g -1- Wx cos 9g) 

+{Mr,' sin 9g + M^' cos 9g) {(j)) - {(pl)] {^'l}}dx 

■[F 3 ®| = J * {[(ilV ^9 ~ ^k' ^a) (^ 2 : 1 ) + i^k' sill ^a + cos 9g) (ti;xi)]{$ J 

-[(7-o)[(£lAi^.- - EAD{wxx){<f>) + IeAD'M + wl) -h ^EAD'Sl)] 

-{GJo - GJr)icl>o) - {fx){<l>x)] m 

-[^E^ADoivl + wl) + lEAD^{<j>l) + (EAD^Vxx - mD;wxx)(l>W^}}dx 

{if} = J^{l(E^.v„-^.v,„M + lEA(vl + v.i) + k-iCM)]W,}}<^ 

The underlined terms in -[Ff },{i4®} and are associated with the axial strain 
at the elastic axis (ux + |u^ + 5 «;;) of the blade. These nonlinear terms are modified 
to a set of linear terms by substituting the axial stress by the axial inertia force. The 
discription of the substitution approach [23,26] is given in next section for clarity. 
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6.3.4 Treatment of Nonlinear Terms Associated with Axicd 
Strain at Elastic Axis 

The equation of motion corresponding to the axial degree of freedom can be written 
symbolically as [23,26] : 

- -Z. = 0 (6.10) 

The axial stress resultant (Equation 5.10) can be rewritten as; 

K = EA[u, + + F (6.11) 

Where F contains all the additional higher order terms, which can be neglected. 
The simplified expression for Vx can be written as 

%^EA[u^+'^vlF-wl\ ( 6 . 12 ) 

Neglecting aU the higher order terms and also the time derivative terms, the dis- 
tributed inertia force Zu (given in Appendix A) can be written as: 

= Ci[{n - l)k + x] + C 2 (6.13) 

where 


Cl = 1- {Pd cos di + Ps sin 9i) 

C 2 = (ei + 62 ) - a{-Pd sin 9i + Ps cos 0 /) 


Integrating Equation 6.10, 




vnFF 


=/ 


Z„dx 


Vx can be written as 


K = /(?)=<iie" + 0,f+a5 



(6.14) 


. (6.15) 


where 
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ai, 02, 03 are given as: 

= ^1-5' 


— [-Ci(X„ - i,) - Cali. 

Jlfo 
02 1 

~ ^e)} "*■ ‘^2][-^V + (-^n — ^e)] 


where 


XN = Y.ile)j = Nh 

j=l 

Xn = Y,{Qi = nU 
i=i 

Combining Equation 6.15 with 6.12, the axial strain at the elastic axis can be written 


in terms of axial inertia force as 


r 1 2 1 21 /(6 


(6.16) 


Using Equation 6.16, the underlined terms in and are modified 

to a set of linear stffiiess matrices, which are given as: 

[' Kt.. {$;} dx=kf^ m {$y di{v} = [iff;] {v} 


I' {«;} m {«'J am = [ifl] {wy 


j‘- {$;} ir = [^] k H m {«;} {«;}^ «*} = K] {*} 

Combining these submatrices, a linear stif&iess matrix is obtained, which can 
be written as 

fK] PI PI P!l 


fR-ffl ^ PI K] PI PI 

^ ^ P! P) [Ki] PI 

PI PI PI PI 


(6.17) 
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where 


['K’n ~ ~ (^1 + ^2 + ~ (2ai + a2)[^2] + (2ai)[vl3] 


[^zsl — (°i + ®2 + <23) [-61] — (2ai + 02) [-62] + (2ai)[S3] 


where 




36 

3/e 

-36 

3/e ‘ 

[^l] = 

1 

3/e 

4/e 

-3/e 

-/f 

30^e 

-36 

-3/e 

36 

-3/e 



. 3/e 


-3/e 

4/2 



36 

0 

-36 

6 /e ' 


1 

0 

6 /e 

0 

< 

[^ 2 ] = 

■ 60/, 

-36 

0 

36 

- 6 /e 





— 6 /e 

2/2 _ 



24 

- 2 /e 

-24 

O^g 

[A 3 ] = 

1 

140/e 

- 2 /e 

-24 

6/2 

2 /e 

2 /e 

24 

-5/e 



. 5/e 


-5/e 

4/2 

3 e 


(6.18) 


(6.19) 


( 6 . 20 ) 


m= 


EAC„ 

ZEAU 


7-8 1 

-8 16 -8 
1-8 7 


( 6 . 21 ) 
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m = 


EACo 

QEAle 


11 -12 1 
-12 16 -4 

1-4 3 


( 6 . 22 ) 




EACo 

30EAle 


23 

-26 

3 

-26 32 

-6 

3 

-6 

3 


( 6 . 23 ) 



Chapter 7 


MODEL VALIDATION 


The first step in any aeroelastic response and stability analysis is the e^’aluation 
of natural frequencies of the rotor blade. Using the inertial and structural model 
developed in this study, a structural djuiamic analysis was performed. It may be 
noted that the inertial and the structural operators given in Equation 6.6 and 
Equation 6.8 respectively, are nonlinear. Since the structural dynamic analysis 
requires only linear terms, all the nonlinear terms are neglected. The corresponding 
linear equation for one beam finite element can be written as: 

Mi{9}i + [jrli{9}i = 0 (T.1) 

where [M\i represents the mass matrix of element, (given in Equation 6.6. and in 
Sec. 6.2.1). The stiflBaess matrix [iTjj consists of three components. They are 
(given in Sec. 6.2.3) , [K^] (given in Sec. 6.3.1) and [K^] (given in Equation 6.17). 
The element stiffiiess matrix is given as: 

[K]i = -f- [K^] + [K^'] (7.2) 

The element matrices are assembled to form the global finite element model for the 
rotor blade. Imposing the root boundaxi' conditions, the corresponding rows and 
columns from the global model are elinnnated. The resulting matrix equation can 
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be written as: 

[M]{q} + [K]{q} = 0 (7.3) 

Performing an eigenanalysis, the natural frequencies of an undamped rotating blade 
in vacuum can be evaluated. 


7.1 RESULTS 

A structural dynamic analysis was carried out for a uniform, untwisted blade, for 
which the data and the uncoupled natural frequencies are available in Ref. [27] . The 
blade data corresponding to two rotor blade configurations are given in Table 1. 
The terminology Soft-in-plane blade configuration indicates that the first nondimen- 
sional lag frequency is less than 1. In Stiff-in-plane blade configuration, the first 
nondimensional lag frequency is greater than 1. 

Using the data given in Table 1, the natural frequencies of a rotating blade are 
calculated. The results corresponding to soft-in-plane and stiff-in-plane blades are 
shown in Tables 2 and 3 respectively. It can be seen from these results that in the 
first flap, lag and torsional modes there is very a good convergence with respect to 
number of elements. The converged results, for first flap, first lag and first torsion 
axe in excellent agreement with the results obtained in Ref. [27]. 



60 


Soft-in-plane blade data 

Stifif-in-plane blade data 

= 0.0004 

= 0.0004 

iTnrjrj/mt^ = 0.0 

IrrinTjlml'^ = 0.0 

9g = 0.0 

do = 0.0 

m = 1 

m = 1 

O 

o 

II 

II 

P 

O 

o 

o 

II 

II 

o 

o 

o 

o 

II 

p 

o 

II 

o 

o 

II 

o 

o 

II 

GJlTr£lH^ = 0.001473 

GJ/mQH^ = 0.001473 

EA/w£lH^ = 20.0 

EA/m^H^ = 20.0 

Cl = 0.0 

Cl = 0.0 

62 = 0.0 

62 = 0.0 

a — 0.0 

a = 0.0 

EACo/EA = 0.00021036 

EAColEA = 0.0008166 

= 0.0301 

EI(^JmaH‘^ = 0.1474 

EI^r,/rnElH'^ = 0.0106 

= 0.0106 


Ofisets of mass centre and tension centre from elastic centre are zero. 


TABLE 1. Input data for soft-m-platne and stiff-in-plane blade 
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Uncoupled Natural Frequencies 


No. of Element 

Ref. [27] 

Mode 

N=1 

N=2 

N=3 

N=4 

N=5 


1st Lag 

.7605 

-7363 

.7331 

.7326 

.7326 

0.732 

2nd Lag 

6.4123 

4.4947 

4.4662 

4.4587 

4.4563 

- 

1st Flap 

1.1507 

1.1301 

1.1261 

1.1250 

1.1247 

1.125 

2nd Flap 

4.4102 

3.4418 

3.4165 

3.4109 

3.4089 

- 

3rd Flap 

- 

8.8800 

7.7140 

7.6605 

7.6376 


1st Torsion 

3.2753 

3.2640 

3.2634 

3.2632 

3.2632 

3.263 

1st Axial 

6.9798 

6.9551 

6.9538 

6.9533 

6.9533 

- 


TABLE 2. Results for Soft-in-plane blade 



Uncoupled NatureJ Frequencies 


No. of Element 

Ref. [27] 

Mode 

N=1 

N=2 

N=3 

N=4 

X=5 


1st Lag 

1.4302 

1.4174 

1.4165 

1.4171 

1.4166 

1.417 

2nd Lag 

13.414 

8.7951 

8.7478 

8.7308 

8.7258 

- 

1st Flap 

1.1507 

1.1301 

1.1261 

1.1250 

1.1247 

1.125 

2nd Flap 

4.4102 

3.4418 

3.4165 

3.4109 

3-4089 


3rd Flap 


8.8800 

7.7140 

7.6605 

7.6376 


1st Torsion 

3.5167 

3.5023 

3.5015 

3.5014 : 

3-5013 

3.501 

1st Axial 

6.9798 

6.9551 

6.9538 

6.9533 

6.9533 



TABLE 3. Results for Stiff-in-plane blade 
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Chapter 8 

CONCLUDING REMARKS 


A general structural dynamic model of a helicopter rotor blade has been formulated. 
This model consists of various geometrical complexities of the rotor system and also 
the rigid body motion of the hub. The rotor blade is modelled using beam finite 
elements having 14 degrees of freedom. The equations of the motion are derived 
using the Hamilton’s principle. Various linear and nonliaear terms associated with 
inertial and structural operators have been derived in a matrix form which is suit- 
able for application to further studies on rotor blade aeroelasticit}' and helicopter 
vibration. 


The analytical model developed ia this study was used to calculate the uncoupled 
natural frequencies of an isotropic, untwisted and hingeless rotor blade. The results 
of the analysis are found to be in excellent agreement with the results available in 
the literature. 
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8.1 SCOPE FOR FURTHER WORK 

Rotor blade aerodynamic model has to be formulated and integrated with the present 
structural dynamic model. The resulting analytical equations can be used for the 
study of aeroelastic response and stability analysis of a rotor blade, coupled rotor- 
fuselage stability analysis and vibration analysis of helicoptors. 
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Fig. 2. Inertial System - R 



Fig. 3. Body Fixed Hub System - H 
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Fig. 5. Rotating Hub System - 2k 








Appendix A 


TERMS ASSOCIATED WITH 
KINETIC ENERGY 
VARIATION 

— 

[7717], n sm{ea + (l)k) + m(m cos{0o + hW-h + w'k) 

{^z - OliPdCOsOj +/?,,sm6*/) - (cos/?; +<?;)} - (cos^; + 9i){6jWk{PdCOs6i 
+(3sSm9r)} + 9iiv'f. cos0;{^x siri'^fc - 9y cos-ipk 
+/?; {Pd sin 9 1 - /?, cos /?;)}] 

+[m7],n sm{9o + (pk) - mCm cos{9o + (f>k)][{^i{Pd < cos 9i + sin 9i > 

H- /?5 < sin 9j — cos 9j > +1) + {Pp -t- Pd cos 9i + Ps sin ^/) } 

{9z — 9r{pd cos 9i + Ps sin 9i) — (cos 9i + ^;)} — {9i{Pd cos 9j + Ps sin 9i) 

+(cos 9j + 9j)} cos 9j{9x cos ipk + sin Tpk) — (cos 9i + 9i)cos 9i9,.w'f.] 

+{7n7]m, cos{9g + 4>k) + wCm sin (^(7 + <pk)][{-^k + ^liPd < cos9j + sin 9i > 

+Ps < sin Oi - cos 9r > +1) + {Pp ^ Pd cos9i + Ps sin 9i)}{9x sin Tpk - 9y cos Tpk 
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+(}j {fii sin Of - I 5 s cos Oj - 1) + 0^) sin 6]) 

+{01 {fh sin Qj - P, cos - 1 + sin dj) cos + Oy sin V’fc)}] 

+[rn?;„, cos( 0 G + (pk) - sin( 0 G + cos 6j - 1 ){^^ - cos Oj 

+/?., sin Bj) - (cos/?; + Bi)] + {cosB, + ( 9 /){(?^cos ^774 - B cos Bj + /?,, sin ^ 7 ) 

+BiVkvQ + Brv'kiBa., sin -ipk - By cost/^^ + (?7(/?rfSin0j - iSsCosBj)}] 

—m[BiiOk + (ci + 62+ < n — 1 > Is + 3:fc)(l — ^7)] 

[B, - B, {Pd cos Bi + /?,, sin Bi) - (cos Bi + 61)] 

+ 7 n[cos Bi (i?,^ sin V’fc - Kj cos - (( 9 ^ + ftd cos 6>7 + A sin ^7) (< n - 1 > Ze + a;*)] 
[Bi{pd cos ^7 + Ps sin Bi) + (cos Bi + ^7)] 

-7/i[(/?,.t: cos V-’fc + Ry sin 'ipk)iPs cos Bi - Pd sin Bi) - 

BiiPsSinBi + PdCos 0i){ei +62) - rukiB^, cos '(j;k + Oy sin -ipk) 

+B^c.osBi{e] +62) - Bi{2{PdCosBi + Ps sin Bi){ei + 62 ) 

+a{Ps cos 61 - Pd sin 61) - Uk} - Wk{Pp + Pd cos Bi + Ps sin 9 i) 

+Uk 4 - n{ps cos Bi — Pd sin Bi) cos 07|[cos Bi + ^7] 

+777,[{7i;fc + a cos Bi + Vk - {Ps cos Bi - Pd sin ^7 + 1) 

(< ?),-!> Is + Xk) - (ci + C2)} 

{4 sin ijjk - By cos ipk + Bi{pd sin Bi - Ps cos ^7)} + {Rx 
+{Bs: sin V-'fc - By cos i)k){< n-l> le + x-jt)} 

[Oi {Pd cos Bi - ps sin l 9 / - 1) sin Z?7}] 

+ 7 ? 2 [/?a.{cos'!/’fc — sin i/)fc(— COS +/?rfsin6>7)} + Rx{—sin'tpk 

- cos ipk{~Ps COS Bi + Pd sin Bi) - sin i)k 9 i{Ps sin Bi + Pd cos ^7)}] 

— m[Ry {sin -ipk + cos'0fc(— cos07 + PdsinBi)} + Ry{cosipk 

- sin 'tpk{-P.s cos Bi + Pd sin Bi) + cos tl}k 9 r{P.i sin Bi + Pd cos ^7)}] 

—7n[R:s {Pp + Pd cos Bi 4 - Ps sin Bj) + RzBi {-pd sin 61 4 - Ps cos 61)] 
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-mlaOj (/?, sin Oj + pa cos (9/) + (/?, cos ^7 - sin Oj) + Uk] 

-rm],n[-Vk c.os{Ba + pk) + v^pk sin(0f7 + <Afc) + v'^^l cos{0a + pk) 

+2Vk<pk sm{9c + (t>k) - ^'k sin( 0 G + Pk) - 2w'kh cos( 0 g + <j)k) 

-tu'kPk cos(^G + (f>k) + ^'k^l sin((9G + (?^fc)}] 

-m(,m[v'k sin(0G + (i>k) + v'^pk cos(0g + h) - sin(0G + h) 

+2v'k^k cos( 0 G + 4>k) - w'k cos(^G + 4) + 2w'k^k sin( 0 G + h) 

+w'k(j>k sin V'fc + w'kPl cos(^G + (Pk)}] 

+ cosf?;[{m 7 ;„, sin(^?G + (yijk) 4- ?n,C,„. cos(^f; + </; 7 ,)} 

{0x + Oy) sin ij)k + 0r. - Oy) cos V->/c} 

-pkiit^Vm cos(0a + (t)k) + niCm sin(6'G 4- pk)}{0x sinT/’fc - By cos ijjk)] 

4-77?. cos Bj [iOjtl (^),r 4- By) sin 'i/jk 4- {Bx - 4) i^k] 

+Wk { B,r. sin iJk - By cos ij)k}] 

[Bx{m.{a cos 6j 4- Vk) 4- mrjm cos(^g + <Pk) + n^Cm sin(^G + h)} 

+B:,{mvk - miim^k sin(0G + A) + n^Cm^fc cos(6»g + A)}] 

4 -^';[{ 7 r?,(ei 4- 02 - wjt) 4- 77777 , ^ sin(0G + A) - ^Cm cos( 6 »g 4- A)} 

{/i, (;os Bi - Pd sin Bj } 4- {7777/, „ c()s(f 9 G 4- A) + "^Cm sin((?G + <j>k)} 

{1- {Pd cos Br + Pa sin 0/)} - mrjm sin( 0 G + Pk) + m^m cos(0g + A) 
+m{acosBi 4- Vk + Wk 4- (ci 4 - 62 — acosBi — Vk){PdCos9i 4-/3sSin0/)}] 

14?/ (!0S^/[{~-7777/)A; 4- 77/.7/,„(/>a.C0s((?g' I' <l>k) 4' 77?4rt<^fc sin(<?G 4* Pk)] [Ps COS <?; 

-pd sin 61} — 0/{77i(ei 4 - 62 — Wk) + rnrjm sin(^G + Pk) — i^Cm cos( 0 g 4 - (j)k)} 
{Ps sin 9 i — Pd cos Bj} + {mvk — mrjmA sin( 0 G + A) + '>^CmA cos( 0 g 4 - (pk)} 
{{ 1 - Pd cos Bi - Ps sin Bj) + Bi{m{ei +67- a- Wk)- mrjm cos(^g + A) 
—m(,m. sin( 0 G + A)} + ^n.7i;jk — mrjmA cos(/?g 4 - A) ~ 'i^'CmA sin( 0 G + A)] 

4- cos Bj[rn{ (ci + 62)^7 {Pa sin Bj + Pd cos ^ 7 ) - aBi sin Bj) - Vk} 
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+)n‘)hn<Pk sin(f^G + 0Jfc) - cos(6»g + (l)k)] 

+ sin 0j[7nwk - rnr]„i^k cos{Oa + (Afc) - m(m^k sin(6>G + <l)k)] 

+{- 0 i + sin 9j) sin 0/{— «4{m77^sin(6»G + (t>k) - TnCmCos(0G + <Pk)} 

cos{6g + h) - 7nC,mSm{dG + h)} 

+m{e] + 62 + < 71— I > le + Xk+ 71, k} - 771(10 j COS0/] 
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[{vnhn s\n{Oa + M cos{Og + (t)h)}w'k 

+{nnhn cos{Oa + ^k) - m(,n sin((?G + </>fc)}v'fc 

+{rmhn sm{Oa + (l>k) - ?».Cm cos{Oa + 0jt)}{0/(/?, cos 0/ - 0^ sin (?/ - 1)} 
+{nn]ni cos{9o + (pk) + n?.C„ sin(6>G + ^k)} 

{(0/(1 - 0d cos Of - /?., sin 0/) + cos 0/} 

+?n{(0/ + 1)(74 + a < 2cos0/ - sin0/ >) + 0 /t(;jt }][02 
+0/(1 + /?rfcos0/ +/5,, sin0/) +COS0/] 

+[{?n,?/„i COS(0(; + (pk) + n<m sin(0(; + (l)k)}{-dz) 

+{rM7/„, sin ( 0/7 + (pf..) - ??<,„ cos(0f7 + </+)}{ 4 sin '0* - 0y cos-0fe) cos0/ 
+?n.(/f,r (X)s ■!/>/. + sin V’fc)][^^/ {(Pi 0/ + /if., sin 0/) + (cos 9j + 0/)] 
+[{n?.7/„i sin(0/7 + (pk) + mCm cos(0(7 + cpk) ]'o'kh 
COs(0c; + (pk) - m(m sin(0G + pk)}Wki>k 
+m,{{-Pg cos 0/ + Pa sin 0/ ) (- i?-,r sin + Ry cos ipk) 

+Rz{(P, + PdCosOj + //,sin0;) + a0/(/?.,sin0; + /?dcos0/) + ujt 
-7/)fc(0;^ sin 7/ifc ~ 0j, cos Ipk) - 0.{n + 7+) } 

-niOi{{-Wk + G\ -t- C2)(/^.,cos0/ - PdShiOi) 

+(ci + C 2 - a cos 0/ — Vk){Pd cos 0/ + /?., sin 0/)} 

- 777 (ei + 62 ) (/3s cos 9 1 — pd sin 0/) cos 9j 

-{77177;^ — 77?,7;,n sin(0G + 0k) + COS(0g: + 0fc)} sin 0/][cOS 0/ + 0/] 
+[{w;^cos(0G! + 0k) + 7?7CmSin(0G + 

+9 1 {Pd < cos 9 1 + sin 0/ > +/?« < sin 0/ - cos 0/ > +1) 

+ (/3p + /3d cos 0/ + sin 0/ ) } 
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+ {mi}rn COs{6g + (f>k) “ mCm sill(^c; + 

+ {rr?.7/,„ sin(^(7 4- (pk) + cos(0c + 

+m{wk + a cos 9] + Vk- (A, cos 6>/ - sin 0/ + 1) 

(< n — 1 > /e + Xk) — (Ci + 62 )}] 

[9^ cos V^A: + sin 'tjjk + 9r{pd < cos 9i + sin 9r > +/?, < sin 9j - cos 9 i > +1)] 
+[{m7],„ cos{9a + + mC„,,sin(0f7 + h)} {cos 9 1 {9^ cos ipk + ijSmi)k)} 

+7n{Rz cos 9 1 + 9x sin -ijjk ~ 9y cos 7/4-) (< n - 1 > Zg + a^fe)] 

[^7 (A/ < cos ^7 + sin > +/?,, < sin (?; - cos (9; > +1)] 

+[{7m/,„ cos(^G + (f)k) - rnCm sin(0f? + <l)k)}{-Vk{9x sin 7/)jfc - 9y cos 
+7;[.(/9, cos 9i — Pd sin (9/)(9/} 

4-{m7;„, sin(0(7 + Pk) + w-Cm oqs{9g + sin 'tpk - &y cos ^jb) 

4-7f4(/9, cos <?7 — sin ^7)^7} 

-?7^{(/i’,^cos'l/'^• + fly sin 'il!k){Pp + pdCosO] + fi, sin 9i) +9j{pdSin9i 

-Ps cos 9i){e\ + 62 ) + 9x cos i/jk + 9y sin 'ipk){a cos 9i + Vk) + 9^ sin ipk - 9y cos -i/ifc) 

(ei + 62) + a(Pp 4 - Pd cos 9i + Ps sin 9i)} + m9i{{pd < cos 9i + sin 9i > 

4/9, < sin (97 — cos (97 > 4l)(a 4 77 a;) 4 77,^ 4 a(/9, cos^j - /?dSin07)} 

4771^7 sin 9i{ei + 624- <n-l>le + Xk + Uk[[9j] 

-7717/, „ [sin 4- Pk){~<^>k + (-4 + ^y) sin V'fc + (9x 4- 9y) cos'tpk 

~9zw[ cos 9f — 9zw'i^ cos 9j 4- 9j{Pd < cos 9i 4- sin 9i > +P 3 < sin 9j — cos 9i >) 

-{-9]{Pd < — sin 9j 4 - cos 9i > +Ps < sin 9i 4 cos 0 / >) + 9j 

+(9,xii^9iui){ej < -/9 rfSin^7 +/9 .,cos6'7 > +1) + 0j (-/3d sin 07 4/3, cos 0/) 
+w'i^ cos 9j — <j)k{—9zv'k cos 0j 4 v[9[ < Pd cos 0/ 4 Ps sin 0/ > 40/77^. 

+9iw'f.v'f^ 4 v'P) cos 0/} 4 cos( 0(7 4 Pk){^k{—^k + 0 * costpk + Oy Sinipk 
~9zw'k cos 9 1 4 9jPd < cos 0/ 4 sin0/ > 40/ A < sin 0/ - cos0j > 40/ 
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■^Ojw'f, < [id cos 0 , + [is sin Oj > + 0 iw'^. + + / 5 ^ cos Oj 

+/?,, sin Bi + w[ cos 9r) — Ozv'k cos 9i 

-9zV'kCos9i + {(hv'k + 9iv'k){9j < -fldSinBi + /?scos0j > 4-1) 

+ ‘9jw[v[. -f 9i{w[,v[, + ii![.v'k) + v'kCOfiB,}] 

-niCm[^in{9a + (/>k){9z^!k 4- 9zV'k - {Biv^ + 9iv'k){PdCos9r + [3, sin 9 j 4- 1) 

-9jv'kw[ - Biv'kw'k - Bjv'kw'k - v'k cos Bj + 4(4 + 4 cos ipk + 4 sin V’fc 
-Bzw'k cos 9 ! - Bipd < cosBj -f sine*; > -4/?s < sin6>/ - cosBi > -9i 
+9iw[. < ftdCosBi 4-/1, sin (?/ > 4-4 '"4 
-< ftp + ftd cos Or + fts sin 9r > +w[ cos } 

4 cos{9(; -f <j>k){(h{(Kvk cos 9, ~ 9,v[. < ftd cos 9, 4 ft, sin Oj > 

-44 “ 4 4 "4 - 

-4- + (4 - By) sin -t/ik - (4 4 4) i’k + 4^Cfc ncs Bj 4 4^4 cos (9/ 

-9] {ftd < — sin Bj 4 cos Bj > 4/1, < cos 9r 4 sin 9j >) — 9r{ftd < cos Bj 4 sin Bj > 
+fts < sin 9] - cosBj > -Bj 4 Brw'kiftdCosBi 4 ftsSinBi) 

■3-9jw'f.{ftd cos Bj 4 ft, sin Bi) 4 &]u)'f.{ftd cos Bj 4 ft, sin Bj) 

+(hic'f, I- 4’'4; — 4(— /lf/‘Sin Bi 4 ft, c.osBi) 4 w4 cos(?;}] 

-m[(/4 + 4/) { - sin 4/c cos Br 4 cos i^kifts cos Bj - ftd sin Bj)} 

+(— 4t + Ry){cos ijik cos Br 4 sin fj^kifts cos Bj — ftd sin 0/)} 

-4 (4 cos ijjk + 4 sin t/fc) (/Is oos Bj 4 ftd cos Br) 

—iiz sin Bj) — RzBj cos Bj — aBj sin Bj) — aB'j cos Bj 

-Br {ft, cos Bj + ftd cos Bj ) (e-i 4 62 ) -f iifc - a4 cos Bj 

-B^j {ft, cos 9 j - ftdSin 6 j){ei+e 2 ) - Wk{^xCOS^k + By sin i^k) 

-Wjfc{(-4 + By) sin -ipk + (4 + 4) cos V’fc} - (4tnfc + BjWk){ftd{cos Bj 4 sin Bj) 
+fts{sin Bj - cos (9;) 4 1} - Bjwk{ftd {- sin Bj 4 cos Bj) 4 ^^(cos Bj 4 sin Bj)} 


- 9 i{{f^d cos 0 i + p, sin 0,){2 < ci + e 2 > + < 7 i - 1 > le + Xk +Uk) 
Jra.{fj,v.osOj - flisinO,) 

-(e, + 62) - (< 71 - 1 > le) - Xk - Uk} 

—di{ 9 i{~~Pd siii + P/j cos ^/)(2 < Cj + 62 > 

+ < n — 1 >!,(,-{- Xk + 77;^.} + [Pd^^OnO] + /?,, Rin(?;)7ijt 
-nOPPs sin Qj 4 - pd cosO,) - 

-{wkiPp + Pd cos Or + P, sin (?/) + WkSp-Pd^^axOi + ps cos Or) 

+a cos Oi)dt {p^ sin Or + pd cos Bj) -|- n^, cos Oj }| 
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[{mr/ni sin(^?c7 + <f>k) + mC„i cos{9a + <l>k)}w'k 
+{vnhncos{0a + (f>k) - M!.CmRin(^r7 + h)}v'k 

-f (j>k) - ni(,n ('.()${()(; + 4) }{(?;(/?, cos (?/ - PhSuiOi - 1 )} 

+{?»,?/„, i-os{0(; + <l>k) + + f^fc)}{(^7(l - PdCos9r - ^sin^/) + cos0/} 

+m,{{9i + + a < 2 cos 9/ — sin 9/ >) + 9/Wk}][9xSini/>kCos9/ 

-9y cos i’k <'OB 9 1 + sin 0i - 9j {p, cos 9i - 4 sin 9i) + 9i + sin 9;] 

+[{rn.?/,„ cos{9a 4- 7nC,nf^in{9(; + h)}{-0z) 

+{?;),?/,„ sin(C4; 4- <f>k) - 7n.C,„ cos(^’4 4- 4)} {4 sin ipk - 9y cos tpk) cos 9i 
+m{Rr, cos V’A- + 4 / sin V’a-)] [^7 iftn cos (9/ - 4 sin ^ 7 ) -9r - sin ^ 7 ] 

4[{nM/,„ sin(r4; 4- cos{9(; 4- (l>k)]v'kpk 

cos(^G 4- <i)k) - mCrn sin(0G + (l>k)Wkh 
+Tn{{~Pg cos 9 1 4- pd sin 9i){-R^, sin V’fc 4- R-y cosV’ik) 

+Rz{Pp + Pd <'OS 9j + /?., sin 9i) + a9i{Pf! sin 9[ + pd cos 9i) 4- iik 
-Wit (4 sin i/)fc - 9y cos V^fc) - 9, {a + n/t)} 

-7n9j{{~Wk + C] 4- C2){Ps('Os 9, - Pdsm9i) 

4(ei 4- 62 - a cos 9j - Vk){Pd cos 9i 4- A, sin 6 * 7 )} 

-?/a(ci (■ C 2 )(/i, c()s/9/ - pdsin 9/) cosOi 

-{rnwjfc — 7TM/„,, sin(0G + <^A:) + ^-Cm cos { 9 g + <^*;)}sin^ 7 ][~^/ ~ sin ^ 7 ] 

+[{nw/rji Rin(^9f? 4- pk) + ”'-Cm cos((?f; + (l>k)}{-h + w'^) 

+{inr},nSh\{9G 4- (j)k) - m(rnC'Os{9G + (l>k)}{^i{Pd < cos 07 4- sin 07 > 

+P, < sin 07 - cos 07 > 4-1) + (/?„ 4- Pd cos 0/ 4- Ps sin 07 )} 

4-{rn.T/,n cos(0c; 4- <j>k) ~ ni4i sin( 0(7 4- <^a;) } Wk0i cos 07 — 1)} 
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-m{0iWk + (ei + 62 + < n - 1 > le + Xk)}{l - Oj cos 6 r)][-9x cos -ilJk - OySinipk 

-01 {Pd < COS0J + sin 01 > +Ps < sinOj - cos0i > +1) 

-{Pp + Pd cos 0} + Ps sin 0;)] 

+[{m?7m sin(6’G + Pk) + mCm cos{0g + <i>k)}{9x cos •ipk + 6y sin tpk) 

-rn{ ilx cos 0i sin V-’fc + I^y cos 0i cos 'ijik 

+{0^ cos 01 + Pd cos 01 + PsSm0j){< n - 1 > le + Xk)}] 

[OiiPd < cos 01 + sin 01 > -\-Ps < sin 0i - cos 0i > +1) - (/?p + p^ cos 0i + P., sin 0i)] 
-[{im],ncos{0a + (f>k) - 7n.CmSin((?r; + (l)k)}{-0zVk cos 0i + 0i{pdcos9i + psSin0i)Vk 
-0iv'i,w'k} + {mr]m sm{0G + pk) - mCm cos(0g + Pk)}{-0zw'k) 

+{rmu sm{0a + pk) + cos{0g + pk)] [0iw'k{pd cos 9i + ps sin (9/)} 

+m{{Rx cos 'll} k + RySm'tpk){PsCos0i - PdSinOi) - 0i{PsSm0i + Pdcos0i){ei + 62 ) 

-xok{0x cos V+ + 0y sin 'pk) + 0z cos 0i{e\ + 62 ) - {Pp + Pd cos 0i + Ps sin 0i)'Wk 
+a cos 0i{Ps cos 01 — pdsin0i)} + 7n0i [a, cos 6 i{ps cos 0i — PdS\n0i) — Ufe}][ 07 ] 

~m[{- iix cos ■ij}k + /?,csin V+ - RySin t))k - R'yCOStpk){Pp + Pdcos0i + /JaSin^j) 
—0i{Rx cos'pk + i?.j/Sint/)fc)(— /?(jrsin0/ + PsCos0i) + R^ cos0i 

1- C2)(/9,/sin - pscos0i) + 0^{c\ + e2){PdCOS0i + psSm0i) + n)*] 

+m[Rx sin sin 0i) — R^ cos pk sin 0i) + Rx0i sin 'ipk cos 9i + Ry cos xl^k sin 0i) 

~Ry sin ipk sin 0i — Ry0i cos ipk cos 0i - a0i cos 0i + a0j sin 0i] 

cos{0(; + f/n,) I- ?nCm sin(^9<; I Pk)]{0x cosV'fc - dx.s\nil)k + Oy sin tpk 
+0yCosrpk + 0i{Pd < cos 01 + sin0i > +^s < sin0/ — cos^j > +1) 

+0]{Pd < — sin 01 + cos 01 > +Ps < cos 0i + sin 0i >) + 0i{—pd sin 0i + Ps cos 0i)}] 
-{{-mr]mSin{0G + <f>k) + mCmCOs{0G + <f>k)}][Pk{0xGOSTpk + 0ySini/'A 
'i'^iiPd < COS 01 + sin 01 > +Ps < sin 0i — cos0i > +1) + /9p + Pd cos 9i + Ps sin^/)}] 
COs(^G + <J>k) - niCm sin(^c? + <f)k)}{pk - {0x sin tpk + 0x cos ‘Ipk 
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-9y COS Tpk + dy sin -ipkWk - i^x sin ^k - Oy cos ilJk)v'k - div'^iPs sin 9i + pd cos dj) 
+v'kiPs cos 9r - 0dSin9i + 1 ) - w',.{9x sin ?/>fc - 9y cos -ipk 

- < ps cos 9 1 - pd sin ( 9 / + 1 >) }] 

-[{‘’TT'Vrn sin (^(7 + 0fc) + mCm COs{9c + (j>k)}{ — ^k + (^x sin'0fc - 9y COSIpk 

- < Ps cos 9i — Pd sin 9i + 1 >)vk - {9x sin tpk + 9x cos ipk - 9y cos V'fc + 9y sin i^k)Wk 
-{9^ sin Ipk - cos 'tpk)w'k - 9j (p^ sin 9j + Pd cos 9i)w'f, 

+{p, cos 9t - Pd sin 9i + l)?^}] 

-m[0x{{<''cos9i + Vk) cofiil)k + (ci + c^) sin i/^*: cos/9/} + /?i{— (acos/9/ + sin 
+{vk + Cl + ca) cos Ipk cos 9i ) + 0 y{ {a + Vk) sin tpk - (ei + 62 ) cos Tpk cos 9i} 

+9y{{vk + C| + (' 2 ) sin V’A; + (« + Vk) cos i/ifc cos (?/ } + 9j{{Pd < -sm9] + cos/?/ > 

4-/9, < cos 9j + sin 9j > + 1 ) (a cos /?/ + Vk) — (ci + eg) — {Ps cos /?/ — Pd sin ^/ + 1) 

(< ?i - 1 > /f, + Xk 4- v.k)} + 9 1 {Of [Pd < - sin 9 1 + cos /?/ > +Ps < cos 9j 4- sin 9i >) 
{a 4 Vk) + VkiPd < cos 9[ + sin 9j > +ps < sin 9i - cos 9i > -fl) 

4^/(^sSin0/ 4 - PdCOs9j){< n - 1 > le + Xk + Uk) ~ Uk{PsCOs9i - Pdsm9j 4- 1)} 
-aBj {Ps sin 9j 4- A< cos 9j) 4- a9}{~Pd sin 9j + Ps cos 9i)] 

cos{0a 4- (f>k) - mC,n sh\{0a + <f>k)}{v'k9i cos 9} + v'^ sin 9}}] 
sin(^G + <Pk) + mCn cos{9o + (pk)}{Wk9i cos9i + w'kSm9i}] 

4?/i/?/ cos /?/[(’] 4" 024" <! ?/. — 1 > /g 4" Xk 4" v,k (i9j{Ps cos 9i PdSioBi) 

+a9j sin 9j){Ps cos 9j — y?dsin ^/)] 
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sin((?G + <t>k) + v[ cos( 6 >g + h)}] 

[0;{cos(0G + <pk) - sin(0G + <pk)} + cos{9a + <l>k)] 

[-<j)k{l7n,,^ + /rn^^){?/4 + h ) '+ ’4 cos(0g + <^fc)}] 

\(),{mn{On + + cos{Oc + 0fc)} + Rin(^G + <?!'fc)] 

-Ivi^c^[{s\n{Oa + (f>k){'f>'k + 0^ + 0, < I ~ PdCosO] - ;0sSin0/ > +cos0/) 

+ (•()s(/?f; + <l>k){-'>i>[.. 4- Ox Kill V^fc (■«>« Oi - Op cos 7j>k cos Oj 
-0i < ps cos 9j - pd sill Oj -\-l>) 

+<?2 + sin Oj } {sin((?G + (f>k){'^i>'k + 0r < Ps cos Oj - pd sin - 1 >) 

+ cos(6'g + <Pk){v'k + 9i <1- Pd cos Oj - Ps sin Oj > + cos Oj)}] 

[{siu(0G + ^k){wk + OxS'mi^k - Oycosijjk - Oj < PsCosOj - PdSmOj + 1 >) 
+ COs( 6 »g + (f>k){v'k " Oi - Oj < I- Pd^'OsOi - PsSUlOi > - cosOj)} 

{sin(6>G 4- (f>k){-v’k 

+dT < I - Pd cos Oj - Ps sin Oj > + cos Oj) + COs(0G + 

-Oj < Ps cos Oj - Pd Sill Oj - 1 >)}] 

— [?n,,i;[(sin(f?G + 4'k){i’'k + Oz + 0] < 1 — pdCosOj — PsSinOj > +cos9j) 

+ cos{9g + (l7k){-w'k + Ox sin -ijjk cos 9j - Oy cos t/>fc cos Oj 
-Oi < Ps cos 0, - /?,/sin<?; -I- 1 >) 

+0z + sin 0 , } {sin{/?f; + (l>k){-v'k + 0, < 1 - Pa cos 9j - ps sin 9j > + cos 9j) 

+ cos{0g + <pk){Wk -Oj <Ps cos9j - Pdsin9j - 1 >)} 

+{sin(0G + (l>k){wi + 0j <Ps cosOj - PdsinOj - 1 >) 

+ COs(0G + (l>k){'f^'k + Oj < 1 — PdCOs9j — PsSin9j >)} 

{sin(0G + M {w'k + Ox sin ipk - Oy cos tlJk-0i< Ps cos Oj - Pd sm Oi + l>) 
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+ cos(^Jg + (l>k){v'k - (K - Or < 1 - Pd cos Of - PgSinOr > - cos^;)}] 

— ?w;,n[{sin((?G + <i>k){Vk + + 0/ < 1 — PdCosOj — PgSinOr > +cos6j) 

+ cos{0g + (f)k){-w'k + 0^ sin tpk cos Or - 0y cos ipk cos Or 
-Or < Ps cos Oj - PdSm Or + 1 >) 

+0^ sin Or) + sin 0t}{{0j + l){vk + a) +0^?x;fc}] 

-mCn[{sm{0G + <l>k){wk + 0^ sin ipk - Oy cos ^k-0i< ps cos Or - Pd sin + 1 >) 

+ cos(^^G + (t>k){v'k - 0:r - Or <l - Pd cos Or - PsSinOi > — cos Or)} 

{{Or + '^){vk + c) + OiWk}] 

+lm.^c[{sm{0G + ^^k){v'k + 0, < I- pdCosO, - psSinOj > +1) 

+ ('Os(f?f,' + f/>A-)( - f'4 ~ fh < P.iCosOi - pdsinO, + 1 >)} 

{sin( 6'(7 + (l)k) {O-r. sin xj^k - Oy cos x))k) cos Or - cos{0g + (j)k)0z}] 

+/r7!,„„[{sin((?G + - 0, < P^oosOr - PdSmOj + 1 >) 

+ cos((9g + 0A:) (v'k - Or < 1 - Pd COS Or - Ps sin Or > -1)} 

{- sin((?G + (i)k)0z - cos{0G + 4>k){0x sin i/’fc - Oy cos xpk) cos 0/}] 

+/?/),, ,J{sin(6'f; + ^k){v'k + 0i <1- pdCosOr - PsSinOr > +1) 

+ cos(^?G + (l)k){-xb'k - Or < Ps cos Or - pdsin 0r + l >)}{- sin( 6 >(j + (l>k)0z 

- cos{()g + (f)k){0si sin xjik - Oy cos V’fc)} 

^■{sin(^'y(; + (l>k){w'k - (h < psCosOi - pdSmOi + 1 >) +cos{0g + (l>k){'^'k 
-Or < 1- pdCosOj - psSmOi > -l)}{sm{0G + MiOx sin xpk - Oy cos xpk) cos Or 

- cos{0g + (pk)0z}] 

+?7).7?,n[{sin(^o + (j)k)iv'k + 01 <1- PdCosOj - PsSinOr > +1) + cos{0g + Mi-^'k 
~0i < Ps cos Or — Pd sin Or + 1 >)}{Rx cos V’fc + sin xpk}] 

+m,Cn[{sm{0G + (Pk){w'k -Or<Ps cos Or - Pd sin 0/ + 1 >) + cos(0g + <l>k){v'k 
-Or < 1 - pdCosOr - PsSinOr > -l)}{i?xCosi/)fc + ^sin^^}] 
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+77^^.,;(;[{sin(^^G + 0fc)(l + Oj) — cos{Oa + <i)k)di} 

{sin(0f; 4- (i)k)v'k^k — cos((?f7 + (pk)w'k^k + sin Oj}] 

-Im^rj[{sm{ec + (i)k)Oi + cos(0G + <f>k){l + 0;)} 

{sin(^G + ^kWk^k + 008(00 + (l)k)vk'i)k + sin^/}] 

+/m^J{sin(0G + + di) - cos(0G + 4)0/} 

{Rin(0G 4- 4 )'»’a: 4 + ('08(00 4- 4)44 - sin0;} 

-{sin(0o 4- 4)4 4- cos(0o 4- 4)(1 + 4)} 

{siu(0o + 4)'(44 ~ (‘os(0o 4- 4)4:4 + sin 0/}] 

47?7.7;m[{sin(0o 4- 4)(1 + 4) - cos(0o 4 4)4}{-4cos0/ 4- 4sin0/ 4- 4 
-Wk{0x sin V’A: - 4 ('OS V’fc) - 4(n + Vk) - (o, + e2)(4 cos 0/ - 4 sin 0/) 

44(-'(nA: 4- ei 4- e2)(4 008 0/ - 4sin07) 

40/ (ei 4- 02 - a - Vk){Pd<~'os6i 4- 4 sin 0/) cos 0/ - Wfe}] 

-mCm[{sin(0o 4- 4)4 4- cos(0g 4- 4)(1 4- 4)}{-4cos0/ + 4sin0/ 4- 4 
-Wk{Ox sin 4 - 0j; cos ilk) - 4(a 4- t/fc) - (ej 4- e2)(4 cos 0/ - 4 sin 0/) 

40/(-?njfc 4- Oi 4- 62) (4 ^^(^s 0/ - /^d sin 0/) 

44(ci + (’2 - a - Vk)(/ii cos ff I 4- /i, sin Oj) cos - Wk}] 

+I7n^(;[{sin{0a 4- <f>k)iOzVk + 4^} < Pdcosdi +4sin0/ > +9jv'kWk) 

-cos(0G 4- 4)(4^«fc 4- Ojxu'k < Pdcosdi + 4sin0/ >)}{4sin(0G + 4)}] 
+/rn^J{sin(0o + 4)(-4t«fc + 4K < 4cos0/ + 4 sin 4 >) 4- cos(0g 4- 4)(4ufc 
44n(. < PdCOsOi 4- 4 sin 0/ > 4-4^^fe'!i'lt)}{“4cos(0G 4- 4)}] 

+/m^J{sin(0G 4- 4)(4wfc 4- 47/fc < 4 cos 0/ + 4 sin 0/ > 4-4<^fe) 
-cos(0g 4-4)(444-44 < 4 cos 0/4- 4 sin 0/ >)}{-4 cos(0g - f 4)} 
4-{sin(0o + 4)(-4Wfc + 4w4 < oos 0/ -F 4 sin 0/ >) -1- cos(0g 4- 4)(4'Ua: 

+4^^fc < 4 cos 0/ 4- 4 sin 0/ > +Oiv'kW'k)}{ei sin(0G 4- 4)}] 
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+?n7/m[{si»(^^C7 + (l>k){Oz'l!k + Ojv’^ < ^dCOSOi + PsSmdi > +0iv[w'^) 

-cos{9a + (f>k){Ozw'k + divt'k < PdCos9i + p,smdi >)}{(< n - 1 > le + Xk){l + 9j)}] 
+mCm[{sin(0G + (f>k){-Ozw',^ + 9iw'k < l3dCos9i + PsSm9i >) + cos{9g + 

+9iv'k < Pd cos 9i + p, sill 9, > +9,v'kWk)}{{< - 1 > /« + Xk){l + 0/)}] 

+Im^d- sin(0G + PkWk + cos{9g + (l>k){-k + {dxcosxpk + 9ySmil)k) 

-9,{pd < cos 9, + sin 9, > +/3, < sin 9] - cos 9j > +1) + {pp + pd cos 9i + Ps sin 9j) + w’k}] 
[sin(6'G + Pk){{-Pk + w'f.) + 9j {Pd < cos 9j + sin 9j > +p, < sin 9r - cos 9i > +1) 

+{Pp + Pd cos 0, + p,, sin 0 , ) } -1- (.■.oi^{()G + h)Wk{0i cos /?/ - 1) }] 

+Iin,,,,[^u\{9(; + pk) { </>k + 0.r cos VU: + 9y sin + 9j {Pd < cos 9, + sin 0/ > 

I /^, < sin 9i vohOi > 1 I ) (/i,, | pdV(M()i I- p„H\n9i) — - cos{0g 4- pk)'ij'k] 

[sin(^/G + Pk){-v'k{9, - 1)} + cos{9o + pk){-Pk + w'k - ^i{Pd < cos 6»/ + sin 6>/ > 

■^ps < sin 9 1 — cos 9/ > +1) + {pp 4- Pd c, os 9} +/33sin0/)}] 

+/?n^(;[- sin((?G + (l)k)Vk + cos{9g + (j}k){-h + {Ox cos tpk + OySinipk) 

-Oi{Pd < cos 9r + sin 9j > +ps < sin 9] - cos9j > +1) 

+{Pp + Pd cos 6, + p, sin 9}) + w'k}] 

[sin((?f? 4- Pk){~v'k{dj -- 1)} -t- cos{0a 4- Pk){-Pk 4- w'k 

-9r{Pd < cos 9i + sin (?/ > 4-/3,, < sin 9j — cos 0/ > 4-1) 4- {Py 4- Pd cos 9j 4- P, sin 0;)}] 
4-/?n,,(;[sin(f'4; 4- <l>k){Pk 4- Ox cos V'fc + Oy sin V4 4- +0i{Pd < cos Oj 4- sin 9j > 

+Ps < sin 9j - cos > -f 1) - {Pp + Pd cos 9j 4- Ps sin 9j) - w'k} - cos{9g 4- (f>k)v'k] 

[sin(f/f; + (l)k){{-k 4- O 4- 0,{Pd < cos 0i + sin 9i > +Ps < sin9j - cos 9r > 4-1) 

+{Pp 4- pdCosOr + Ps sinOi)} 4- cos(0g 4- Pk)WkiOi - 1)}] 

+mr}„i[— sin(<?G 4- pkWk + cos(6>g + Pk){~Pk 4- {Ox costpk + OySinipk) 

-Oj {Pd < cos Oj 4- sin Oj > +Ps < sin Oj - cos Or > 4-1) 

4-(/3p -f Pd cos 9j + Ps sin 9j) 4- nifc}] 
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+ (f! + <' 2 + <n~ l>l,+ Xk){l - e, COS0/)}] 

+mCm[sin(0f; + cos + By sin xjjk + 9i{0ct < cos 6j + sin dr > 

< sin Oj - cos 0/ > +1) - (/?p + cos Or + sin dj) - - cos(0g + h^k] 

[-{9]Wk + (Ci + 62 + < n — 1 > + Xfc)(l - 9i COS^j)}] 

+/m<J-?4 sinida + (j^k) - cos(0g + (l)k){k + 0/(/?d < cosdi + sin dj > 

+/?, < sin di - cos ^?/ > +1) - (/3p + cos 9: + sin dj) - njjt}] 

[sin(6'G + (i>k){0x cos ijik + sin 

\ f VI,,, i\s\n(0(; I (j)k){<h- I dii/li < cokOi I sill > -f-/3, < sin — cosOj > +1) 

-{0p + Pti cos Or + A, sin 9j) - xv'k} - cos((?g + h^k] 

[- COs(/9(; + (4 cos 0jt ()y sin x/>k)\ 

+Im,^(;[-Vk mi{da + (j>k) - cos(0g + (i>k){^k + 9r{Pd < cosdj + sin 9j > 

< Hindi - v.ondi > +1) - {Pj, + PdC-OHOj + PsSindi) - u^}] 

[- cos(0G + 4>k){dx cos V-’jt + Oy sin xpk)] 

+/mpJsin(6>G + Pk){Pk + 9i{Pd < cosdr + sindj > +Ps < sindi - cosdj > +1) 

-{Pp + Pd cos dr + Ps sin dj) - tcjj.} - cos(6'g + <?ifc)nfc][sin(^G + Midx cosxpk + 9y sin xpk)] 
+fnx}m[-Vk sin {da + <pk) - cos{da + <i>k){Pk +9i{Pd < cosdr + sm9j > 

+P, < sindj - cosdr > +1) - (/3p + /?rfcos0/ + Asin^/) - 'U)jfe}][-^a:SinV’)fc cos^z 
+11,1 (X)s xpk cos di + {d^ + cos 9} + ps sin /?/)(< n — l>lr. + -t*;)] 

+7nC,„[sin(<';G + (l>k){i>k + 0i{pd < c.oh9, +sin<?z > +/?, < sin0; - cosdj > +1) 

-(/ip + 0d cos d] + /i., sin dj) - xu'k) - cos{do + PkWkl 

[-/?x sin xjik cos di + lly cos xjjk cos dj 

+{dz + Pdcosdj + Ps sin 9 j){< n-l > le + ^A;)] 

+Im^c^[—di cos{da + (^A:)][sin(0G + Pk){^z'^k cosdj + djWkiPdCosdj + psSindj)} 

+ cos(0G + (l>k){-0zv'k cos dj + djv'kiPdCOsdj + Ps smdj) - djv'kW'k)] 
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+hn,,,j\(h sin(fV; + (l>kWn\{(^a + cos0; - (/J^cosl?; + /?,sin 0;) 

+(),v[.w[.} + coH{0a + <h-){0,vi cohO, + 0,m[{p,icos9, + p.sinO,)}] 

cos{0a + (l}k)]\sin{da + <j)k){0zv'k cosOj - eiv[{PdCos9i + /JjSin^;) 

+0iiiw'k} + co^On -f- (f>k){9zvicof^9, + 9,uMicos9, + p,sm9i)}] 

+Im,r,c[di sm{9c + (?ifc)][sin(0G + (f>k){0zw'kCos9i + 9iw'k{pdCOs9i + PsSm9i)} 

+ cos((?(7 + <l^k) (h + 9jv'^.{P,i COS 9j + Ps sin (9/) — 

+m.r]m[-0r cos(^fi' + cos'0a: + Ry sin tpk){Ps cos 9 r - PdSin 9i) 

-9i{ps sin 9j + pd cos 9i){e\ + (‘2 + 2asin 9j) — Wk{9x cos^|Jk + 9y sin'ipkWz cos 0/(ei + 62 ) 
+ 0 , {a cos 0 , < 1 % (H)s 0 , ~ pd sin 0 i > -Uk) - {pj, + pd cos 9i + p^ sin 9i)wk 
+a cos 9 1 (/9., cos 9/ — p,i sin 9i)] 

+w4,[/')/sin(^9f; + (/)k)][{Rr('o^i’k + Ry sm i!k){Ps cos 9r - Pd sin 9 1 ) 

-9i{p., sin 9; 4- pd cos 9i){c] 4- (’2 + 2a sin 9i) - Wk{9x costjjk + 9y sin tj)k)9z cos 0/(ei + 62 ) 
+9i{a cos 0i < p,, cos (9/ - pd sin 6 ] > -Uk) - {Pp + Pd cos 9] + Ps sin 9i)wk 
+a cos 9j{Ps cos $j - Pd sin 9j)] 

+Ini(^(;[sin{9f; + ^^fc){?4(cos^?/ + sin/? 7 )((')r sini^fc - 9y cos 'ij^k 
+9] < pi,cos9r -- /9,;sin^9/ >)} 

+ c.os(f ?(7 f (f>k){ -w[.{9,: sin i/’fc - 9,, cos i/’fc + ^/ < P.<, cos 9, - Pd sin 9j >)}] 

[f); cos(r9(,- + </>k)] 

-\ Im„„\{w[.s\ni9a + ch) f v[. (M9a \ Pk)]icos9, + sin 9,) [9x sin i^k - 9y cos i^k 
+9 1 (/I, cos 9] ~ pd sin 9i ) }] sin(6i(; + (j)k)] 

+/m„Jsin(6>G + (j}k){vk{cos9i + sin6'/)(4sinV;fc - 9yCosxl>k 
+9j < p_^cosO] - Pd sin 9 1 >)} 

+ cos(0G + (l>k){-w'k{9xSin'ipk - 9y cos ipk + 9i < PsCOs9i - PdSin9r >)}][0jsin(0G + 4 
M(^g + 4>k) + UfcCOsC^G + Pk)}{cos9i + sin0j){4sinV’fc - ^yCos^Pk 
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+0i{fi, cos 0i - 0^ sin 0i)}][0j cos( 0 c + ^k)] 

+ 777 , 7 /, „ [sin + (f)k){cosO, + sin 0i){v'^{0^, sin ^Jk - OyCosipk + Bt < /J^cos^j - ; 5 dsin 0 / >)} 
+ cos{9a + (lyk){cose, + sin 0j){-w'^{9r.simpk - 0yCosipk + 9i < P^cosBr - PdSinBi >)}] 
[< n - 1 > 4- .r*] 

+rnCm[{t(4sin(<?G + (l>k) + v’^cos{Ba + h)} {cos 6i + smBi){B^simPk - ^yCosipk 

+Bl{0s cosfl, - /?,,sin ^ 7 )}][< 77. - 1 > lr + Xfc] 

+Im(^(_[sin{0a + (i>k){i>k — Ox cos i/jk - Oy siiiipk — 0i{0d < cosBj + sin 0 / > 

+/?., < sin 0; — cos /?/>+!+ 77 ^) — {pp + fl,i cos 6r + Ps sin ^ 7 )} — Biw'^ cos(^g + ^k)] 

[sin((?(? + (j)k)(h'^0k + cos{0a + i>k){-0k + 0i{Pd < cosBr + sinBi > 

+/?,, < sin 9r - cos 0[ > + 1 ) + (/?p + A/ f^os Bf + A, sin Br) + OiVk}] 

+lmy,,[0rw{.sm{0(:! + (fik) 4- cos{0a 4 <j)k){-h 4- BxCosipk + Oy sin ipk 
+Bi{pd < cos Oj 4- sin 0/ > +0^ < sin Br - cos Bj > +1) 4 0!v[ 

+{Pp + PdCOsOj + 0, sin 07 )}][sin(<?G 4 (f)k){-h 4- 6i{Pd < cosOi + sin0/ > 

4At < sin Or - cos Of > +1) ■+ {pp 4 Pd cos 0i 4 A sin Or) - OiVk} 4- Bjw'k cos(Bg + <i>k)] 
4/7npAsin((?G + ^k){^k- Ox cos T/'jt - 0y sin ijjk - 0j{Pd < cos Bi 4 sin Bj > 

4 A^ < sin Or - cos Or > I- 1 -f- 74 ) -- {0p I- Pd cos Bj 4 P, sin Bj)} - Bjw'f. cos{9a + </»fc)] 
[sin{0a 4- M{~^k + OiiPd < cosOj +smBr > +Ps < sinBj - cos0/ > +1) 

4(A, 4 0d cos Or 4 0., sin Bj) - 4 Biw[ cos{Bg 4 <f>k)] 

4/7n,,,Asin(A; + <l>k)(hfi>'k 4 cos(0c 4 0k){~^k + OiiPd < cos Br 4 sin Br > 

4 A < sin Or - cos Or >41)4 (A> + 0<i cos A 4 A sin Oj) 4 0rv',,}][Brw'k sin(0G 4- (f>k) 

4 cos (^(7 + <f)k){—^k + Ox COS V^fc 4- By sin ipk + 0i (A < cos Br 4 sin Br > 

4A < sin 0; - cos > 4 1) 4- Orv'k + {Pp + A cos Br 4 A sin 9r ) }] 

4m7/„,,[sin(0G 4- (j>k){h - Ox cos V4 - A sin V’fc " A(A < cos Br 4 sin Br > 

4A < sin Or - cos > 4l 4 v’k) “ (A + A cos Or + A sin Or)} - Orw'^ cos(0g + (f>k)] 
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[w/fc + acof^O, + - {/i, cosOi - PdSin 0, + 1){< n - 1 > l^ + Xk) - (ei + ej)] 

+mCm[^/«4 sin(0o + ^k) + cos(6Ig- + (l)k){-i>k + O^cosipk + OySinipk 
+di{lh < cos 9 1 + sin 9j > +/5,, < sin Oj - cos 9} > +1) + Orv^ 

+(/^p + III cos 9 1 -f 0, sin 9 1 ) }][7/;fc + a cos Oj + Vk- {0^ cos Oj - 0^ sin 9j + l) 

(< 7), - 1 > Ir + Xk) - (ei + 62)] 

+/7n'CcN**(^’^f7 + 0k){<f>k ~ 9i{0d < coFi9j + sin^; > +0g < sm9f — cos9i > +1) 
+()iv[. - ( 0 ,, \ 0,1 cos 9, +0,Hm9,)} - co^{ 0 a + 0k)9}w'0 
[(•os((/<; I- (l>k) cos 0, {9, cos t/'k + 9y sin 7 /)^)] 

+/ 777 ,,,,[sin( 6 ’r 7 4- + cm{9a 4- 0k){-i>k + Oi{0d < cos9r + sm9i > 

+II, < sin 9t cos 9 1 > •+ 1 ) -4- 9iv^ + {0^, -f 0^ cos 9i + /?,, sin 9 } ) }] 

[sin(/?^; + (f)k) cos 9; {0^ cos t/jjfc 4- 9y sin 7 /)^)] 

\Iv)y0s\u{9(; I <l>k){>hk - (h{0,l < COsOi -l-sillfl; > +0., < sill — 008(9; > +1) 
+(9/74 - iftv + f^d cos 9i + 0, sin 9j)} - (‘Os{9a + 0k)Oiw'k] 

[Hin{f9f; + (I>k)i9:r COH 'tj^k + Sill '07- )] 

+ /77?.,,c;[sin{(9(-; + (/>k)9[Wk + C0S((9f; + 0fc){-0fc + 9i{0d < COS 9i + sin 9i > 

+0S < sin 9; - cos 9/ > +1) + (9/74 4" (/^p + f^d cos 9j + 0s sin 0/)}] 

[C0S{(?(; + (f)k){9r. cos V’7 + 9y Sill ^fc)] 

4-//).7/„,[sin(r9(/ I (/)/.) {07 - 9,{0d < cos (9, +sm9j > +0s < sin (9/ - cosOj > +1) 
+(9/74 "■ W]> 4- 0d cos (9; + 0s sin (9/)} - cos((9(; + 0fc)(9/u4] 

[/?2 cos 9 1 + {9s-. sin '07 — 9y cos V-’fc)(< n — 1 > h + ^k)] 

+777,0,„[sin((?(7 + (j>k)9i'i0k 4“ cos(( 9(7 + 0fc){-0fc 4- 9i{0d < cos^j + sin^/ > 

+0S < sin (9/ — cos 0/ > +1) + 9jVk + {0p + 0d cos 9i + 0s sin 9j)}] 

[9^2 cos 9j + ((?j; sin t/j^ — 9y cosipk)i< n-l> le + ^fc)] 

+Imc^(^[~-9j sm{9a + 0fc)][sin(0G 4" <l>k){-'^'k{^xS^^'^k - ^yCOS'ipk) 
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-}-w'f,{fi,vi)sO, - fi, I (),)() j} 

+ cosiOa + <Pk){-v’k0ar.sm ipk - 4^:08 -i/jfc) + w;^(/3,cos0; -^rfSin0j)0j}] 

+lrn,0Oi cos{fla + 0 *-)][.sin(^?G + <f)k){vk{0T.sm'tpk - OyCostpk) 

!)(),} 

+ vM0g + </>A:){-<(^)Tsin - ()yCosi)k) + w'^iPsCosOi - I3dsmei)ej}] 

sin(6^f? + (j)k)][sm{dc + (l>k){v'k{OxSinil)k - dycostpk) 

-v'kiftsCOi^Oi - fl,is\nO,)(),} 

+ C()s(r/f; + (f>k){-1l>[.{(}rf^mij'k ~ COR ■(/'<.■) + ’»fc (A cos 6', - /5rf Rin }] 

+lnh,^[d, cos{()a + (j>kWM(h + (l>k){-^40x sin tjjk - Oycos'ipk) 

+w[.{li/,vosOi 0])(}i} 

+ coH{()a + (f>k){-A0xSin i^k ” OyC-osijik) + VkiPsCosOr - Pasm9j)9j}] 

+[-7mj„0i sm{Oa + </>k) + niCnOj coh{9g + ^k)][-{Rx cos Tjjk + RySinipk) 

0p + ftd cos Oi -f /?., sill Oj) + 9i{ft,i cos 9i - ps sin0j)(ei + e^) 

+4 COR V’A; 4- Oy sin 7j’k)0<'Ofi 9] + 71k) + 0x sinTpk - 9y cos i/'fc)(ei + C 2 ) + Uk 
+a{pf, cos 0, ~ fii sin ^?/ ) } + a(/5p + Ai cos 9j + A sin 9j)] 

+}v^(^[{{-v'k - w'k(h) sii^A; 4- (j)k) + {-v'kh + w'k) cos{9g + A)} 

4- (h-)iw[. 4- Asin ipk c()s(?/ - 9y cos ijik 00 s 9 j 
-Oi < psCosOi ~ A/ sin/?/ — 1 > 4- sin/?/) 

4-cos((?(; 4- Mi-v'k -9z-0i <1- l3dCOs9j + Asin^/ > - cos A}] 

+Iniyy[[{7)'^.(^f. — n’4) 810(0/7 4" “ 0'k + '^'kM COS(0G 4- <f>k)} 

{sin(0G 4- (i)k){v'k - A - A < 1 - A cos A - Asin0/ > -cos9i) 

4- cos(0G 4- 4>k){-‘0i 4- /?x sin 7j)k - 9y cos V'fc + A < A cos 9i - A sin A + 1 > - sin A)}] 
+Im,,^[{{-v'k - w'M sin(0G + h) 4- Hkh + 4) cos(0g + h)} 

{sin(0G 4- <l)k){v'k - A - A < 1 - A cos A - A sin A > - cos 9i) 
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+ cos{ea + (f>k){-K + k sin - Qy cos-^fc + < /3, cos 0/ - sin + 1 > - sin dj)}] 

+{{v'kk - uQ sin(0G + <l>k) - {v'k + w'^^k) cos{9g + (i)k)} 

{siu(/?f; + 4>k){w'k + Ox sin V'fc cos Or - 6y cos ipk cos 0/ 

-9j < (5s cos 9j - sin (9/ - 1 > 4- sin 9j) 

+ cos(0G + 4>k)[-v'k -9^-9i <1- cos 9r + /?* sin 9i>- cos 0/}] 

+nnhn[{{-v'k - K^k) sin(6'G + h) + {-v'kh + ^'k) cos(9g + <^fc)} 

{R^ cos ij^k + Ry sin ijjk - cos 0;(acos 9T + Vk){l + 9i) - sin 9TWk]] 

+n'C>i[{ (njt<^A: - w() sin(0G + (}6fc) - {v'k + w'^^k) cos(9o + (('k)} 

{R^ cos ij)k + tiy sin V^fc - cos 9 ] (a cos 9i 4- t;fc) (1 4- 9i) - sin 9jWk]] 

+Im^c_[{vk sin(<?G + ^k) + u/^ cosV-'fc}{sin(0G + <(k){wk+ < 9a-. 4- 9y > sin ipk cos Oj 
+ < 9x — 9y > cos V-'fc cos 9 1 — 9 1 < l3s cos 9i — Pd sin 9i — 1 > 

-h9j < ps sin 9i + Pd cos 9i > 

+^k < v'k + 9z-\- cos 9 1 > +9i^k <^- Pi cos 9i 4- /?« sin 9i >) 

4- cos(<?G + (j>k){^k <Wk + 9x sin 'ipk - 9y cos 'ipk > -9i^k < Ps cos 9r - Pd sin 9j - l> 
-tifc - dx- Bj <1- pd cos 9i 4- Ps sin 9r > -9] < Pd sin 9i 4- Ps cos 9] >)}] 

-lin„„[{v'k cos{()a + Pk) + w'k sin V^jt} {sin(^>G + h){v'k - 9z 
-9j <l - Pd cos 9 1 + Ps sin Or > -O] < Pd sin - Ps cos 9i > 

-<j)k < -w'k - sin Oj) 4- OxSinijjk - Oy cos ipk > -OiPk < Ps ^^osOj - /3d sin + 1 >) 
4 - cos ( 6 »g + Pk){^k < v'k -Oz- cosOj > -Ojpk < 1 - Pd^osOi - PsSinOj > -w'k 
■{■ < 9x 5- 9y > sin -ipk cos 0i+ < Ox — Oy > cos V'fc cos Or 
+9[ < Ps cos 9[ — Pd sin 0/ 4- 1 > 

4-^; < -psSinOr - PdCOsOj >)}] 

4-/m,,^[{t;^ sin(0G + Pk) + w'k cos ipk} {^'^^{O g + Pk)ivk "" Oz 
-Or <1- Pd cos 9j + Ps sin Oj > -9] < Pd sin 9i - ps cos Oj > 
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< -w'k - sin Bi) + 4 sin t/'fc - By cos > -Bi^k < Ps cos 0/ - sin 0/ + 1 >) 

+ cos(^f7 + (l>k){<f>k< v'k -Bz- cos Bj > -Br^k <l-Pd cos Br - sin Bi > -w'^ 

+ < (K. + By > sin V’fc cos 6j+ <6^-ey> cos cos 0; + 6i < p, cos Bj - Pd sin 0/ + 1 > 
-\-B] < -/?.,sin(9/ - PdCOsBi >)}] 

—l7nij^[{vi. cos{6c + <pk) + 'Wfc sin V'/:}{sin(^G + <f>k){w'i^+ < Bx + By > sin'ipk cos 6j 
+ < Ox- By > cos Tpk cos - §1 < p, cos Bj - Pd sin 6»/ - 1 > 

< /?,, sin 0; + Pd cos 9f > 

+Pk < v'k + + cos (9/ > < 1 - pdCosB, + /?«sin0; >) 

4- cos(0G + Pk) i<j>k < w'k + Bx sin Tph - By cos "tpk > —Bjpk < Ps cos Bj — pd sin Bi — 1 > 
“’4 “ - h <l - Pd cos Bj + Ps sin Bj > -Bj < Pd sin Bj + Ps cos Bj >)}] 

+7m]m[{v'k sm{Ba + pk) + w[ cos^fc}{(-i?,^ + R,,) sini/’fc + (Rx + Ry) cos^jt 
-nA;(cos Oj +B[ - (a + Vk)Bj - Wk sin 0/}] 

-m,Cr7i[{Vk COS{Bg + (f>k) + W'k Sm7pk}{{-Rx + SinV’Jfc + {Rx + Ry) COSTpk 
-Vk{cos Bj + Bj - (a + 7Jk)Bj - Wk sin Bj}] 

+I77i^(^[{(^k cos{Bg + <j)k)}{sm{BG + (l>k){-h 
+01 c.os 'ij)k + By sin '0fc “ Bz7u'i^ cos Bj + Bjpd cos Bj 

+ sin Bj > -\-Bjps < sin Bj — cos Bj > +Bj + Bj7ju[ < Pd cos 6j + Ps sin Bj > +6jw'j^ 

+ < Pp + Pd cos Bj + Ps sin Bj > +70^ cos Bj) + COS(0G + ^k){—Bz'^'k ^OS Bj 
+ v'^.Bj < Pd cos Bj + Ps sin Bj > +Bjv',^ + Bjiv'^v^ + njj. cos 0;)}] 

-Imppli^k sin(0G + h)}{sHBG + Pk){BzVk cos Bj - 9jv[ < Pd cos Bj + Ps sin Bj > 

-jj'Jj _ - t 4 cos Bj) + cos(0f; + (l>k){-k " COS 7j)k - By sin 7j7k + 0.w]t COS 6j 

—Bjpd < cos Bj + sin0/ > —BjPs < sinBj — cosBj > —Bj 
+0/n/fc < Pd cos Bj + Ps sin 0/ > 

+BjWk- < Pp + PdCOsBj + /?sSin0j > +WfcCOS0/)}] 
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+/777,„c:[{0fcCos(6»G + (?!»fc)}{sm(6'c; + <l>k){OzVkCosei - Biv’k < PdCosdi + Asin6'/ > 

-v'kOj - 6 1 Vf.Wk - Vk cos di) + cos{0o + h) {-h - 9x cos t/'jfc - 6 y sin tpk + Ozw'k cos Oj 
-Brpd < cos 6 , + sin 6 ; > < sin Bj - cos Bj > -Bj + B:w', < /?, cos Bj + /?, sin Br > 

+Biw[- < flj, + 0d cos Br + /3, sin Bj > +w[ cos ^j)}] 
sin(0G + (j)k)}{sm{BG + 

+<?.p (X)s V'fc + By sin V'fc - cos Br + Brfta cos Br 

+ sin Br > +Brfh < sin Br - cos Br > +Br + Brw[ < Pd cos Bj + /?., sin Bj > +6rw',^ 

+ < Pp + Pd cos Br + Ps sin Br > +?i;[ cos Br) + cos(0g + 4>k){—Bzv'r^ cos Br 
+v'j^Br < Pd cos Br 4- Ps sin Br > -^Brv'^ + Bjw'j.v'r. ‘^k cos ^j)}] 

+7/77/„,[<;^fc cos(^?fj + ^ 7 fc)][/?.x{-sinV'fcCos 6 >/ + cos ipkiPs cos Br — PdSinBr)} 

+7?.„{c()s ijjk cos Or + sin i>k{Ps cos Br - pdSm Br)} - sin (9/ 

-(h { [Pn cos Or + Pd cos Br ) (ci + Ca) } + h - aBr sin Br 

-Wk{Bx cos Tpk + By sin i)k) - BjWk{Pd{<^os Bj + sin Bj) + ps{smBj - cos Bj) + 1} 
-Br{{PdOQsBr + /?., sin Br){2 < Ci + 62 > + < n - 1 > ig + 

+a(/?.,cos<?/ - Pd sin Br) 

-(ei + C 2 ) - (< 71 - 1 > k) -Xk-Uk} - {{Pp + PdCOsBr + PsSinBi)wk 
-a cos Or {Ps cos Br - Pd sin Br) + cos Br{ei + 62+ <n — l>le + Xk + «*)}] 

-viCmlRx { - sin V’fc cos Or + cos 4>k{Ps cos Or - Pd sin Br ) } 

-I /?,/{cos ij)k cos Or I sin il>kiPn cos Or A; sin (?;)} - 7 ?,, sin Or 
~Br { {Ps cos Or + pd cos Or){ei + 62) } + Vk — O’Br sin Or 

—Wk{0x cos t/’fc + Oy sin i/’jfc) — BrWk{Pd{^<^^ Br + sin Br) + /?s(sin Or — cos Br) + 1} 

—Br{{Pd cos Or + Ps sin Br){2 < ei + 62 > < n — 1 > le + Xk + Uk) 

+a{Ps cos Br - PdSinBr) 

-(ei + 62 ) - (< n - 1 > le) - Xk - n,} - {{Pp + PdCOsBr + PsSmBr)wk 
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-a(X)s^^;(/i, cos(?/ /idsinf^/) +cos 6 >/(ei + 62 + < n- 1 > Ig + Xk + Uk)}] 
+Irna{^m{Oo + </>k)][sm{(h + h){~k + (-0. + ^„)sinV>jfc + + 9,) cos^ 

-§,w'k cos 9, - 9,w'k cos 9i + 9i{pd < cos 9i + sin 9j > +/ 3 , < sin 9i - cos 9i >) 

-V9]{(5d < - sin 9j + cos 9r > +/?, < sin 9i + cos 9j >) + 9j 

-\-(9]XUj^ + 6jxhi^)[9i < Pd sin 9j + Pg cos 9j > -t-1) + — Pd sin 9i + Pg cos ^/) 

+u4 cos Or - M-Gzv'k cos 61 + v'^qOi < pd cos Oj + pg sin 9j > 

+div'^ 4- (pxv'kv'k + v'k cos 0/) } 

+ cos(9(; + cos x/>k + Oy siti x/)k - 9^xVk cos9r + 9ipd < cos 9r + sin 9i > 

+01 Ps < sin 9 1 — cos 6j > +9i + 9}u/k < Pd cos di + Pg sin 9i > 

+0] xo[, + Pp + Pd <-os 9 1 + pg sin 9, + xu'^ cos 9j) - 9^Vk cos 9/ - cos 9/ 

+(0rx/k + (9r < -pd sin 9j + /?,, cos 9i > +1) 

+ '6,wWk + 01 + iPk^h) + v'k cos 9j}] 

+Impp[c.os{9G + (/>fc)][sin(^c? + Pk){0zVk cos (9; + OgV'k - (Ojv'k + Ojv'k) 

{Pd cos 0! + Pg sin (9/ + 1 ) - 9iv'kXv'k - Oiv'kxn'k - 9iv'kxh'k 

-v'k cos 9 1 + pkiPk + 0T cos xj)k + Oy sin xl>k - Oz^i'k cos 6] - 9iPd < cos 9j + sin 9j > 
—Oipg < sin 0] — cos 0i > —Oj -t- (hxo'k < PdCosOj + PgShi9i > -\-9jXJUk 
- < pp + Pd cos 9 1 + pg sin 9j > +'w4 cos 9,)} + cos{9g + pk){M9zVk cos 9i 
-9iv'k < Pd cos (9/ -I- Pg sin 9i > -v'k9] - 9iv'kXii'k - v'k cos 9i) - (pk 
+ {9.r. — 9y) sin xpk — {9.,r -I- 9y) cosi/’a: !■ 9zV>'kV-^>^9i + OgXj/kCosOi 
-9]{Pd < - sin 9 1 + COS0/ > 

+pg < cos 6 >/ + sin <9/ >) - 9j{Pd < cos 9j + sin 9i > +pg < sin 9i - cos 9r > 

-9 1 4 9iw'k{Pd cos 9i 4 Ps sin 9i) 4 9ixh'k{Pd cos 9j 4 Ps sin 9j) 

+9'iXu'k{Pd cos 9i 4 Ps sin 9i) 4 9ixv'k + 9ixh'k 
-9i{-pd sin 9i 4 ps cos 9j) 4 xh'^ cos 9i}] 
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+/7n„(:[sin((?G + </>/=)] [sin {6 »g + 4){0.WfcCos0; + O.v'k - {6jv'^ + Oiv'k) 

(/^rfCOs^?/ + fl, sin + 1) - (hv'kW'k - (hv'kw[ - djv'kw'k 

—v'k C'OS 6j + 0fc(0fc + cos '0fc + Oy sin ipk — cos — OjP^ < cos 0/ + sin 9i > 
-OjP, < sin 0, - cos Oj > -Or + em/^ < PdCosOi + ^,sin0; > +9iw'k 
— < Pp -h Pd COS 9; + Pg sin 9j > +w'k cos 0/)} + cos(^g + </’ik){<^fc(^zn[ cos 9j 
-9jv'k < Pd cos 9 1 + Pg sin 9i > -v'kOi - 9iv'^w'k - v'k cos 9j) - ipk 
+(^x ~ ^4) sin V>A- ~ {9g., + 9y) cos tjjk + cos 9j + 9zw'k cos 0; 

-(PiPd < - sin 0i 4- cos9i > 

+pg < (;os 9i + sin <?/>) — 9j (/j,^ < cos 9j + sin 9j > -yPg < sin 9j — cos 9j > 

~(ji (),w'^{p,icoi^()i pgs\n9,) + 0iw'k{Pd<^oa9i + pgsm9j) 

+9]v{.{pii cos 0, + pg sin (?/) + 

-9r{-pd sin 9j + Ps cos (?/) + 1/4 cos 9 ] }] 

+/rn^^ [cos (^?G + (l)k)][^m{9c + pk){-pk + {-0x + Oy) sin + {6^., + 9y) cosipk 
-9zxu'k cos 9j — cos 9j + OpPa < cos 9j + sin 9i > +Ps < sin 9j — cos 9i >) 

+9]{p,i < - sin Oj + cos 9i > +Pg < sin Oj + cos 9j >) + 9j 

+(9jw'k + Oiv{P0i < -pa?-\n9) +/i,cos(?; > +1) + ^;(-/5dSin(?/ 4- cose*/) 

4(c^, cos 9i — cos 9] + v'l^Oj < p,i cos 9j 4- Pg sin 9j > 

+9jx{. 4- Ofw'kV'k 4- r>'kcos0,)} 

+ cos(6>g + <Pk){<Pk{-k + O^gcospk + OySin^Jk - 0,w'kcos9i 4 Ojpd < cos9i + sm9i > 
+0}pg < sin 9j — cos 9j > +9i + 9jw'k < Pd cos 9i + Pg sin 9i > 

+9jw'k + py + pd cos 9 1 4- Ps sin 9j + w'k cos 9j) - 9,v'k cos 9r - 9,v'k cos 9i 
+ 0iv'k){9i < -/?dsin^/ + (5,,cos(9/ > 4-1) 

+9iw'kv'k + Opw'kv'k + w'kV'k) + v'k cos 9i}] 

+mVm[sm{9G + <l>k)][{Rx + Ry){- sin V’it cos 9r + cos MPs cos 9^ - Pd sin 9i)} 
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+(-7?,. + ny){cosxl>kCosfl, +smi}k{P,v.osO, - /?,isin0;)} - R,sm9r) - RJicosOj 
- 0j{Rx cos V'A: + Ry sin i)k){P, cos Oj + cos Oi) 

-aO, sill 0 i ) - aOi cos Oj - 0 , { 0 , cos 0 , + 0 ^ cos 6 j){ei + 62 ) + h 
-O'j) 0 ,cosO, - Pdsme,){ey+e 2 )-wk%cos'ipk + 0 ysmi}k) 

-Wk{{-dx + Oy) sin ipk + (Ox + Oy) cos 'ipk} 

-{OiWk + OiWk) {/?rf(cos 9j + sin 9r) + 0s{sm 9j - cos 9i) + 1 } 
sin 0/ + cos/?/) + 00cos9[ 4-siii0/)} 

-0, { { 0,1 cos 0, t 0, sin (),){2 < C) + c^ > I- <7i-l>l^ + Xk + Uk) 

+a(/?,, cos 9 1 - 0,1 sin <?/) - (ct + (•■ 2 ) -{<n-l>le)~ Xk - Uk] 

- 0 i {Oi{- 0 ,i sin Oj + /?., cos 9j){2< c, + c-i > + < n - 1 > + Xk + Uk} 

+{ 0,1 cos 0, + 0., sin ()i)uk - aOi{lh sin 9, + 0,icos 6 ,) - m*} 

- {Wfc { 0 p + 0 d cos 9 1 + 0 s sin 9 1 ) + WkOi {- 0 d sin 9i + 0 s cos 90 
+a9i cos 61 { 0 s sin 61 + 0dCos9j) + iik cos 9;}] 

+mCm[(-^.a: + 7?,/) { - sin l/^A: COS 0/ + COS l/.)k{0s COS 9 [ - 0dSm9j)} 

+{—Rx + ij,/){cos^fc cos9} +smiJk{f^sCos9r — 0dSm9i)} — RzSm9i) — Rz9j cos9j 

-9/ (/?.,. cos V’A- + l\, sin V’fc) (/?s cos 9i + 0 d cos 90 

—o.(0 sm90 — a(?f cos/?/ — /?/ (/7,, cos 0/ + 0dcos90{ei + 62 ) +Vk 

—900s cos 9] — 0,1 sin 90{e\ + C 2 ) — cos V-'jfc + 9y sin V-’fc) 

-Wkii-Ox + 9y) sin 0 k + (Ox + Oy) cos 0 k} 

-ijhwk + 9,Wk){0d{cos9, + sin/?/) + 0s{sin9i - cos90 + 1} 

-9jwk{0d{— sin 9i + cos90 + 0s{cos9i + sin 0 /)} 

-9i{{0d cos 9j + 0s sin 90{2 < ei + e 2 > + < n-\ > + Uk) 

+a{ 0 s cos 9 1 — 0 d sin 00 — (ei + 62 ) - (< n - 1 > le) —^k — 

-O!{9,{-0d sin 9j + 0s cos 90{2 < ei + 62 > + < n - 1 > le + ^k + Uk} 
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+(/5rf cos Bj + p, sin Bj)uk - adi% sin Bj + cos Bj) - u^} 

- { iPp + Pd cos B, + p, sin Bj ) + WkS, (-/?d sin 6r + /?., cos Bj) 

+aBi cos Bj (/i, sin Bj + p^ cos Bj) + cos0;}] 

+/m<c[0fcSin(/?G + (A^-)]Nn(0c? + '?!»ifc){-4wj;sin4cos0j + ^,t/;;cos^fccos0/ 

+Biw',,{Ps cos Bj - pdSinBj + 1 ) - w[ sin Bj} 

+ cos{Bg + (l>k.){+Pk + 0x(cost/)fc - KfcSinV’fcCos^/) + 0j,(sin'0jfc 4- v^cos ipj, cosBj) 

+h{Pd < cos Bj + sin Bj > +Ps < sin Bj - cos Bj > +1 

+v[. < p, cos Bj - p,i sin f?, -f 1 >) -t- p,, + p,, cos Oj + /?., sin Bj - v',, sin Bj}] 

+l7n,jrj[<Pk: COS ( 0(7 + </’fc)][sin( 0(7 + ^)k){-^k + 0i(cOS Ipk + v[ sin Pk cos Bj) 

+0,Xsin V'fc - cos V’fc COS Bj) + Bj{Pd < cos Bj + sin Bj > -\-p, < sin Bj - cos Bj > +1 
< /?, cos Bj — Pii sin 0; + 1 >) + /3p + /3d cos 0f + Ps sin Bj + v[ sin 0;} 

+ cos(0(7 + (/>*•){ -OxV’k 7j)k + ByW'^ cos '0A; + BjWkiPs cos Bj ~ Pd sin 0; + 1) — w'jc sin 0/ }] 
+lm„(^[^k sin ( 0(7 + (/>fc)][sin( 0(7 + (l)k){-<i>k + 0i(cos V’fc + v',,smpk cos Bj) 

+0,/(sin 'ipk — v'j^ cos V'fc cos Bj) + Bj{Pd < cos 0/ + sin Bj > +Ps < sin Bj - cos 0j > 4-1 
-Vk < Ps cos 0 / - Pd sin 0 / 4 - 1 >) 4- /3p 4 /3d cos 0/ 4- Ps sin 0 / 4 v'^ sin 0 ;} 

4 cos(0(7 4 0x^4 V-’fc + 07/^4 cos'ipk 4 Bjw'f.{Ps cosBj — pdSin 0/4 1) — sin 0/}] 

4/m,,^[0fc cos(0(7 4 ^fc)][sin(0(7 4 pk){-^xWk sin pk cos Bj 4 Byw',. cos pk cos 0/ 

40/?i4(/3s cos Bj — Pd sin 0/ 4 1) - ^ sin 0/} 

4 cos(0G 4 Pk){+Pk + 4(cos 'pk - v'k sin pk cos Bj) 4 0y(sin pk + v[ cos pk cos 0/) 

+0r(/3d < cos Bj 4 sin Bj > 4/3., < sin Bj — cos Bj > 41 

+V',, < P, cos Bj - Pd Sin Bj 4 1 >) 4 ^p 4 /3d cos 0/ 4 Ps sin 0/ - v[ sin 0/}] 

+7mu[k sin(0(? 4 </>ifc)][-i?x{cos pkiPp + Pd cos Bj 4 Ps sin Bj) 4 sin Pk sin Bj} 
-RyisinPkiPp 4 PdCOsBj 4 PsSinBj) - cosPkCosBj} 4 R^cosBj 
40/ {Pd sin Bj - ps cos 0/) (ci 4 62 ) 
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—aOi cos 0 i + ii>k + Ox { <‘os •i/.’k {a cos 0 i + v/.) 

+ sin ijjk 0 [(ei + 62 + <n — l>l^-\-Xk + Uk)} 

+^„{sin cos 0, + Vk) - cos xl^k cos 0,{ei +62+ <n~l> l, + Xk + «fc)} 

+0i{Pd < cos 0, + sin Oj > +/?, < sin 9, - cos Oj > +l)(a + Vk) 

+{ei + 62) (A, cos 0/ ~ PdSin 6 [) 

-(/?., cos Oj - fii sin + l)(ei + 62 + < n - 1 > + a;fc + Uk) 

+a{(3, cos 0, - f3a sin 9,)} + a{pp + cos 9i + p, sin 9j) 

+ sin 0i){('\ t- C2 f- < V ~ 1 > I.,. + Xk 4- 11 .^,)] 

+mCm [<k cos 0A-][-- { cos x/)k{p„ + Pd cos 9, + P, sin 9{) + sin -ipk sin 9i} 

-ny{s\n%lik{p,, + PdCos9i + P,sin0,) - cosi/>fc cos0;} 

+/?;; (X)s Oj + 0i {P 4 sin 0] - /?., cos 9i){ey + 62) - ak cos 9j 

+ Wk 4 4:{cOS |/'A-(o COS/?/ + Vk) + sitl V’A COS (?/ (cj + 62 + < H - 1 > Zg + Xfc + Wjfe)} 

+9y{sin xj}k{a cos (?; + Vk) - cos ifjk cos 9j{ey + 62 + < n - 1 > k + Xk + Uk)} 

+k{Pd < cos 9j + sin (?/ > +/?., < sin 9j - cos 9i > +l)(o + Vk) 

+(ci + 62 ) (/i, cos 9j - p,i sin 9i) - (/?,, cos 9} - Pd sin 9i + 1 ) 

(ci + f’24- < 11 — 1 > Ip + Xk + Uk) 

+a(/^, cos 9] — pd sin ) } + o.{Pp + Pd cos + P,, sin 9i) 

+ sin Oi){^Ci + C 2 + < 71 — 1 > Ip + Xk + Uk)] 

cos{9c + 0A:)][sin(<?(7 + </>*;) {cos 0/ (sin 7/jfc < —k - k ^ 

+ cos ij^k 

- sin 9jWk + {-k. sin ipk + k '^k)w'k cos 9i + {9jWk + kw'k) 

{Ps cos Or - Pd sin + 1 ) 

—9jiu'k{Ps sin Or + Pd cos Or) — ^kik + < cos ipk — Uk siniik cos Or > 

+k < sin iPk + v’k cos V+ cos Or >) - kWd < cos Or + sin Or > +Ps < sin 0 / - cos Or > 
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+1 + cos 6 1 /?rfsin + 1 >) — sin di) — ^k(^/3p + (3d cos 9i ^3- Ps sin ^/)} 

+ cos(6'o' + ^k){^k{-0xw'k sinV'Acos^/ + 6yw[cospkCos6i 
+0iiv'i. < /?, cos Or - Pd sin 6'/ + 1 >) 

- sin Oiw'k + ^k + {Ox + Oy){cos pk - sin pk cos Oj) 

+(-0x + ^,/)(sin ipk + v'k cos pk) 

+{-9x sin pk + Gy cos pk) cos Ojv'k + Or{Pd < cos Oj + sin Or > 

+Ps < sin Or - cos Or > 

+1 4- v[. < Ps cos Or - Pd sin + 1 >) + 0'j{Pd < - sin Or + cos Or > 

-\-ps < cos Or + sin Or > -v'^ < /?., sin Or + pd cos Or >) 

+0rv'k < (i,cos0, - PdsmOt + 1 > 

+0t{-pd sin Or + /?., cos Or) + ?)[}] 

+Im,,,j[sm{0a + (l)k)][sin{0G + (j)k){-h + {Ox + Oy){cospk + sin ptc co^ Or) 

+{-0x + Oy) + (sin - v'k cos pk cos Or) 

+(<'4 sin pk - Oy cos pk) cos Orv'k + 0i{Pd < cos Or + sin Or > 

+/?,, < sin Or — cos Or > +1 - v'^ < P,, cos Oj - Pd sin Or + 1 >) 

+0‘j{pd < — sin Or + cos Or > 3-ps < cos Or + sin Or > +Vk < Ps sin Or + cos Or >) 
-(h'K'kiP.iCOsOr — /ii,||Sini?/ + 1) + (-/?,/ sin 0/ + P,,cosOr) 4 -nji;Sin^/ 

-^ki-Oxw'k sin i^k + Oyw'k cos V-'jt 

+0rv/k < Ps cos Or - Pd sin ( 9 / + 1 >)} + cos(^>g + <l>k){hi-^k 
+0X < cos ij)k + v'k sin pk cos Or > +0y < sin pk — v'k cos pk > cos Or) 

+0ipk{Pd < cosOr + sinOr > +Ps < sinOr - cosOr > 

+1 - v'k < Ps cos Or - Pd sin + 1 >) + ^k{Pp + Pd cos Or + Ps sin Oj) 

+(- < Ox + Oy > sin 'tpk+ < -Ox + Oy > cos pk)wk 

+(~^3:Sin'0A: + COS ) ic jj. + {Orw'k + 0rWk){Ps cosOr — PdSinOr) 
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-w',s\n0,) - (yfw'MsinO, + (j.cosO,)}] 

+(-(^x + dy) + (sin 'tpk - v'^ cos cos dj) 

+(4 sill V^s - cos *) cos (?,»[ + |j, (/J^ < c„s ^ ^ 

+A < sin e, - cos «, > +1 _ sin 9, + 1 >) 

+9?W, < -sin,?, +cos<i, > +A < e„s«, +sin«, > +„; < A sin 9, + ft cos 9, » 
-/?/7A.(/1,cos 6;; - /Irfsinri; + 1) + ^M-/lrfSin0, +^,cos(?;) +^;sin^; 

sin V’A: + COS V^a, 

+(?An.-' < cos 0 , - fi , sin <?; + !>)}+ cos(^yc + <t>k){M-k 
MK < <-os V'A- + 7-(. sin V'A- cos r;, > +4, < sin V?, - n' cos > cos 0;) 

+^iMlk < cos 0A + sin Oj > +/5,, < sin Oj - cos 0/ > 


+1 - r'(. < /I, cos (h - fh sin 0^ + 1 >) + ^os + /?, sin Oj) 

+(- < 0x + Oy > sin i)k+ < -k ^Qy> cos 'tPk)wk 
+(-4 sin ijk + Oy cos i;k)w'k + (Oiw', + e^w '^) (/?, cos Oi - sin Or) 
-w'k sin 0 , ) - O'fn/^ {0, sin Oj + pd cos }] 


+lni,,(‘[mi{0a + </>a:)] [ sin (/?(-? + (/AA){cos6i;(sin V^*. < -0^ -0y> 

+ cos (/'A- < -r),, -f 0y >)w[ - sin 0,Wk 

+{-0r sin V’A- + 0y cos i>k)w'^ cos Oj + {Oiw^ + 0iWk){p, cos Oj - Pd sin 0i + 1) 

-0iWk{P.., sin 0} + pd cos 0j) — <j)k{^k + 0X < cos i/jk - v'^. sin -ipk cos 0j > 

+0y < sin V'A + v(. cos V’A cos 0J >) - (PkO, {Pd < cos 01 + sin 0, > +0, < sin 0r - cos 0j > 


+ 1 + Vj. < /I, cos 0j — Pd sin + 1 >) — (Va? 4 sin 0i) — ^k{Pp + Pd cos 0j + pg sin 0j)} 
+ COS{0G + 0a) {^k{—0xWk sin V’A cos 01 + 0yW'k cos V’A cos 01 
+0jw'k < /?_, cos 01 - Pd sin + 1 >) 

- sin 0jw'^ 4- 0A + {0x + 4) (cos 0^ - v'^ sin V’A cos 0i) + {-Ox + i9j,)(sin 0 a + 4 cos 0 a) 
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+( K. sill V’fc + COS V-'fc) COS Qjv'k + 0 i [fid < cos 6 i + sin 6 j > +/5j < sin — cos di > 
+1 + Wfc < Pa cos 6j — pd sin (9; + 1 >) -f 0‘^[p^ < _ gin 0j 4 . cos 0/ > 

+Pa < COS Or + sin 9i > -v'^ < p^ sin 6 t + Pd cos 9, >) 

+9iv'k < Ps cos 9 1 - Pd sin (9/ + 1 > 

+0i{-pd sin Of + Pa cos 9i) + i;^}] 

+nnu[- cos[9c + <f>k)][{[Rx - Ry) sin pk - [Rx + Ry) cos'ipk} 

[Pp + Pd cos 9j + Pa sin 9i) 

—9i[n,rV.osij'k t- sin '0fc)(— /9,/sin(9; Pa cos 9 j) 4-iizCos^; + 9r[PdSin9i — PgCos9i) 
(ci + C 2 ) + 0/ [Pd cos 9j + Ps sin Oj) + Wk - aOj cos 9j H- [6^, + ^y){cos ipk[o, cos 9j + Vk) 

- sin V’Ar cos 9 1 (f'l -4- ('2+ <)).-!> + Xk + TAjfc) } + [-9^ + ^j/) {sin pk (a cos 9i + Vk) 

- cos ijik cos 9i{ei + 62+ < n - 1 > /^ + .-rjt + Ufc)} + 9o,{vk cos ipk + Uk sin ipk) 

+9y{vk sin V’fc - ilk cos i>k) + 0i { [Pd < cos 9j 4 - sin 9i > +Ps < sin 9i - cos 9i > + 1 ) 

(a cos 9i +Vk) + (ci + 62 ) [Ps cos 6 j - Pd sin (?;) - [Ps cos 9i - Pd sin Oj 4- 1) 

(ei + 62+ <n-l>le + Xk + Uk) + a{Ps cos 9i - PdSin 9i)} 

+9'^j{{Pd < - sin 9i + cos 9i > +Ps < cos 9i + sin 0/ >)(a + Vk) 

- (ci + f’2) [pa sin 0i + pd cos 9 1 ) + (/{, sin 9, + Pd cos 9]) 

(Ci + 62+ < n - 1 > + Xk + Uk) 

-a[Pa sin 9i + pd cos 9/)} + 9i{{Pd < cos (9/ -f- sin 9i > 

+pa < sin 9i — cos <9/ > +l)?)fc — (/?,, cos0/ — PdSm9i + l)wfc} 

+a9j{—PdSm 9i + psCOs9j) + UkSin9!) + 9jUkCos9j] 

+7nC,^[sin(^G + </’A:)][{(^x — Ry) siuifk — {Rx + ^) cos'ipk}iPp + Pd cos 9 1 + Ps sin 9i) 
— 9 i[Rx cos ifik + -/?.y sin't/>fe)(— /9dSin0/ + PsCOs9j) + RzCOs9i + 9i{PdSin9j — PsCos9i) 
(oi + 62) + 9] {Pd cos 9 r + PsSin9j) + Wk - a9iCOs9i + {9x + 9y){cosipk{(icos9i + Vk) 

- sin iJk cos 9j{ei + 62+ < n- 1> le + Xk + Uk)} + (-4 + ^y){sini/>fc(ocos^r + Vk) 
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- cos V’fc cos 0i (ei + 62 + < n - 1 > /g + Sjt + Wfc) } + 4(^jt cos i/'fc + sin i/>fc) 


+^j/('!’'a- sin V’fc — n/; cos V-'t) + { (Ai < cos + sin B] > +A, < sin {?/ — cos Bj > + 1 ) 


(acos(?; +?;/;) + (ei + e2)(Acos0; - /^rfsin^/) — {Ps cos Bj - /3d sin 0/ + 1) 
(ci + 62 + < 71 - 1 > k + Xk + Uk) + a(A9Cos0; — P^sinBr)} 

+Bi{{P,i < cohB] +sin0/ > 

+P,, < sin Bj — cosBj > +l)iA — {Ps cos 0; — ^rfSin 0/ + l)ufc} 

+ 00 / {—Pd sin 0; + P, cos 0/) + wjt sin 0;) + 0/iijt cos 0/] 



sin(^r; + + ^k){ii>'k + 4 sin V^jt cos 9, - 9y cos ipk cos 6r 

-9j < Ps cos 6i — (3d sin - 1 > 4- sin 9j) 

+ cQs{9a + (l)k){-v'k - 9^ - 9] < I - Pd cos Oj + Ps sin 9 1 > - cos0/}] 

+Im.,pj[{dkCOs{9a + (^fc)}{sin(^G + ((k)d''k — 9^ 

-6] <1 — Pd cos 9] — /?., sin 9i > — cos 9i) 

+ cos(^?c7 + (j)k){-w'i. + dxsui'ipk - Oycostj)k + 9] < Ps cos9j - Pd^in9j + 1 > — sin^/)}] 
+Imy(\{^k sin(6>G + 0fc)}{sin(6'G + ((kdv'k - 9z 
-9 1 <1 — Pd cos 9 1 — Ps sin 9j > — cos (?/) 

+ cos(0G + (f>k){-w'k + 9xSini/)k - ^yCostpk + 9^ < PsCOs9i - Pdsm9i + 1 > — sin^/)}] 
+{^k cos(0G + (fjk) } (sin {9g + ?!>fc) {w'k + 4 sin ipk cos 9i - 9y cos i>k cos 9i 
-Oj < Ps cos 9] — Pdsin 9i — 1 > +sin<';^/) 

+ cos(0G + dk){-v'k - 9z - 9j < 1 - pd cos 9j + Ps sin 9i > - cos 9j]] 

+7m]m[{^kSm{9G + (f>k)}{Rx cos ipk + Ry sin Tpk - cos9j{acos9i + nfc)(l 4- 9t) - Wksin9i}] 
+m((m[{i>k cos{9g 4- ((>k)}{Rx cos Tpk + ^j,sin'0fc - cos9i{acos9i + Vk){l + 9r) - Wk sin0/}] 
4- /»?(;(; [cos (^4; + <Pk)[—(h + (hiPdCos9, + /?,, siti 0/ 4" 4- 1) 4- COS<?/}] 

[sin(0G 4- <l)k){-<(>k + Wfc) 4- 9j{Pd < cos9i + sin0/ > 4-/3* < sin^/ - cos0j > 4-1) 

+{Pp 4 - Pd cos 9r + Ps sin 0/)} -f- cos(^g 4 - </>*) {9r cos 9j — 1)] 

[sin(0G 4- 4>k){9z - 0/(/?dCosi9/ 4 - PsSin9i -w'f. + l)- cos0/}] 

[sin(6'G 4- (i>k){Vki^ - 9} cos 6>;)} 4- cos(6'g 4- (()k){{-^k 4- w'k) 

~^iiPd < cos 9i + sin 9j > +Ps < sin 9i — cos 9j > -f 1) — {Pp 4- Pd cos 9i 4- p, sin 9j)}] 
4-/m^Jcos(0G + <j7k){—9z 4- 9j{PdCos9i 4- pssin9i 4- 4- 1) 4- cos0/}] 

[sin(0G -i- <f>k){Vk{l -(?/)} 4- cos{9g 4- (l>k){{-h + w'f.) 
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—(h {lii < cos 0i + sin 6] > +fl, < sin Oj — cos Oj > +1) — (/3p + 0^ cos 6j + 0s sin Gj)}] 
+hn,,i^{sm{0G + (l>k){Gz - Oj{{0dCosOi + 0sS\n6] - w[ + l) - cos (9/}] 

[sin((9c; + <Pk){—^k + w'/.) + 6}{0fi < cos(?/ + sin Gj > +0s < sin Gi — cosGj > +1) 

+(/5p + Pd cos Gj + 0s sin 0/)} + v'k cos{Gg 4- (l)k){Gr - 1)] 

+771?7,„[cos(^?g + (pk){-0z + Gi{0,icosO, + j0,sin Gj+w[ + l) + 1}] 

[—GiVJk — {ci + 624- < Ti — 1 > Zg 4- .x^){l — Gj cos &[)] 

4-mCm [sin (6 'g' 4- (l>k){Gz - Gr{0dCOsGj 4- 0sSmGr - wj|. 4- 1) - 1}] 

[ — OlWk — (ci 4- 024- < 7?. — 1 > I,, 4- .rA:)(l - 0i cos Oj)] 

+lvi(^(^[cos{0a + (f>k) [G I {0d cos 0 1 -t- /?., sin/9/ 4- 4- 1) -f cos Oj] 

{sin(/9fj 4- (t>k){Gs. cos 4- Oy sin V’/t) }] 

4-/mpp[sin(0(; 4- <j)k){-Gi{0dCos6] 4-/3., sin 0/ -w[ + l) - cosGj} 

{- cos{Gg + (l>k){Gx cos V-'fc + Gy sin -ipk)}] 

+Irnr,<^{cos{GG 4- (f)k) {Gi{0d cos Gj 4- 0sSm Gj + w'^ 4- 1) 4- cos i9;} 

{ - cos (/?(7 4- (f)k){GT. cos i’k 4- Oy sin tj)k ) } 

4-sin(/?f7 4- (f)k){-Gi{0d cos 0[ -f 0s sin Oj - w[ 4- 1) - cos0/} 

{s\\\{Og -t- (l)k){Gx cos V’A- 4- Oy sin ■(/'fc)}j 

4-7717;, „ [cOs(/9f7 -f (l)k){Gj {0d cos Gj -i- 0s sin 0/ 4- 7/7j[. -f 1) 4- 1} 

[—Rx sin tjjk cos 0] '+ Ry cos cos Oj 

+{G^cosGj 4- 0dCOsOj 4- 0sSmO]){< n-l>le 4-xa)}] 

4-777Cm[sin(0c' 4- (l)k){-Gi{0d cos Oj 4- 0s sin ^/ - 4- 1) - 1} 

{ — i?.,; sin t/’a cos 4- cos cos 0/ 

4-(02 cos 4- /3d cos 0/ 4- /3s sin 0/)(< 77 - 1 > /g 4 - xa)}] 

4-/777(;^[cos(0g 4- (f)k)0- + ^7)]lsin(9G 4- (t>k){-GzWk CosGi 4- 6 j w’,,{0d cos Gj 4-/3,sin 0/)} 
4- cos{Og + 4>k){-Gzv'k cos Gj + Giv'f.{0dCQS 6j -1- 0s sin Gj) - Gjw'kv'^}] 



109 


sin(6^G + (j}k)[l + f?;)][sin(0G + (t>k){0::v'^ cos0/ - 9iv'f.{ftd cos Gj + A sin Gj) 
■^Gjw'f.v'k] + cos (<?(7 + ^k){Gzw'k cos G, + Gjn/^ilS^, cos Gj + /?., sin 6’/)}] 

+hnnc_\^os{Ga + h){^ + ^/)][sin((?G + (pk){Gzv'k<^o?,Gi - Giv'^i^dCosGi + PgSinGr) 
+9,w'^v'f,} + cos((9g + (l)k){Gzw'k cos 0, + Giw'^if^d cos 6j + sin 0/)}] 

+/77?.^J-sm(6'G + ^k){l + ^/)][sin(6»G + 4>k){-G^w'kCosGi + 9jw'k{0dcos Gj + 15^ sin6>/)} 
+ cos(0G + <i>k){-0z'^'k cos^/ + Giv'kiPdCosGj + PsSinG,) - Gjw'kv[}] 

+Tn77m[cos(0G + <?^fc)(l + ^/)] [(i?i COS + i?,y sin COS 0 / - PdsinGr) 

-GiiP^sinGj + PdCos 6 ,){ei +62 + 2asin Gj) - WkiG^cosil^k + GySmipk) 

+ 0 zCOS 0 /(ei + 62 ) 

+Gi{a COS G, < /?5 cos (?/ - Pd sin Gi > -Uk) - {Pp + Pd cos Gj + /?,, sin Gi)wk 
+a COS Grips cos 9r - Pd sin Gj)] 

+mCm[-sin(0G 4- 0fc)(l + (?;)][(7iiCosV^fc + Ry sin 'tpk){Ps cos Gj - PdSinGr) 

-Gj(Ps sin Gr + /?rfCos^/)(ei 4- 62 + 2 a sin ^;) - Wk{Gx cos i>k + GySinipk) 

+GzCosOr{ei 4- 62 ) 

+Gi{acosG] < PsCosGr - PdSinGr > -Uk) - {Pp+^dCOsGi + /?., sin 
4« cos Gj [ps cos Gj — pd sin Oj)] 

+Im(^(;[cos{Gc 4- (?ifc) (cos Gj + sin Gr){ —Gy, sin ijkk + Gy cos V-’Jt + Gj {Ps cos Gj — Pd sin 0/) }] 
[f)fH4sin(/?G + (^k) + cos(('1g + (l^k){—4>k 4- GiiPd < cosOj +sin 0 } > 

+Ps < sin 61 - cosGj > 4-1) + Pp + PdCOsGj 4 Ps sinGr + Gjv'k)] 

+lmjjr,[sin{Ga + (l)k){cosGj 4- sin0/){^i sin^fe — Gy cos xpk — Gj{Ps cos Gj — pdSinGi)}] 
[sm(0G 4 - <i>k){-^k + hiPd < cos Gr 4- sin Gj > +ps < sin Gj - cos Gr > +1) 

+/^p + Pd cos Gr + Ps sin Gr — 0/?4} + cos(0g + <^fc)] 

+Imyc^[cos{Ga -f <j>k){cos9t 4- sin^ 7 ){- 0 isin tpk + Gy cos ipk + GiiPsCosGr - PdSinGr)}] 
[sin((9G 4- + GtiPd < cosGj + sin^r > 4^0, < sin^j - cosGi > 4-1) 



+/Jp + fidi^os 0! + (3s sin 0, - 0iv'f.] + 0/w4 cos[0o + 0*)] 

+ hn,,clsmi0a + (f>k){co^0i + sin ) { 4 sin ?/>fc - 0yCosi)k - 0 i{PsCQb 6} - 0dsm9j)}] 
ie,iv'k fim{0a + (f)k) + cos ( 6 'g + (t)k){-^k + 0i{Pd < cosOj + sin0/ > 

+/?, < sin 0] - cos f ?7 > +1) + ftp + Pd cos 0j + /5,, sin 0j -f 0iVk}] 

+mr],n[cos{0a + (pk) (cos 0] + sin di){-0x sin + Gy cos ij)k + Oi iPs cos 0/ - pd sin 0i)} 

{icjt- + acoR0i + Vk - iPs COS0J - Pdsm 0, + 1)(< n - 1 > Ig + Xk) - {ei + e^)}] 

+mCn [sin {0G + (pk) (cos 0j + sin 9 ] ) { 4 sin rpk - 0y cos rpk - h [Ps cos di - Pd sin 0/ ) } 

+ acos0i + Vk - {Ps COS0J - /^rfsin Oj + 1 )(< n - 1 > + Xk) - (ej + 62 )}] 

+Im(^/^[{cos{9G + ^k){cos0[ + sm0j){0j < PsCosOj — PdSmOj >)} 

{cos(i?G + <Pk ) cos 01 {0x cos ?/ife + 0y sin xf)k ) }] 

+/m,^({sin(6>G + (pk){cos9[ + sin6'/)(-(9/ < Ps cosOj - PdSinGi >)} 

|sin(<?f; + (l)k) cos 0, {0r cos i’k + 0y sin Ipk) }1 

+/m,,(;[{cos((9G + </>fc)(cos0/ + sin0/)(^7 < PsCosOj - PdSinOi >)} 

{sin(^G + <pk) cos 0r(0s: cos ipk + 0y sin V-'jt) } 

+{mi{0G + <pk)icos6j +sm0i){-Oi < pgCosOj - Pdsm0r >)} 

|cOH{r;f ; + (pk) cos 0, {0s: COS V’fc + Oy sill (/’fc) }] 

+mT/m[{cos(^G + <Pk){cos9i + sm0i){0i < PsCOsOi — pdSinOj >)} 

{/?; coi^Oi + Oj sin t/jfc — 0y cos V^fc)(< — 1 > /e + 3:^)}] 

+rnCm[{sm(0G + (pk){cos0i +sm0,){-0i < PsCOs0r - PdsmOj >)} 

{Rs cos 0j + 0x sin V’fc - cos ‘tpk){< n - 1> le + iCfc)}] 

-■ ! n)(_(_[(j)k sm{0(; + (pk){fim{0a + <pk){‘>’k^k - n'4 4- Ox sin i}>k cos 0i — 0y cos tpk cos 0j 
-(L sin 0i) - 9j < pg cos 0f — pd sin 0i > +9i + sin 9i) 

+ fos(^G + (pk){-v'k - w'k^k - Oz COS dj - 0! <1- Pd COS 6i - Ps sin 0i> - cos 0/)}] 

-InirjylPk cos(0G + !;&fc){sin(0G + (pk){v'k + w'k^k — ^zcos - 0/ < 1 - PdCOs9j[ 



Ill 


-fj^sindi > - cosOj) + cos(0(7 + <j>k){v'k^k - u>k - 9xSmtpkCos6i 4- 9y cos tpk cos 9} 
+0,sm9,) - 9] < -P,,cos9r + /3rfsin0/ + ^7 > -sin0/)}] 

-Inhf({<f>k^'M9a + (l^k){^o\{9c + (th-.){v’k + w'^.^k - 9^cos9j -9i <1- /3dCOs9i 
-f3f sin 9} > - cos 9r) + cos( 6 'g' + (pk){v'k^k - w'k - 9^ sin i)k cos 9i + 9y cos -ipk cos 9i 
+ii.sin^;) -9] < -P,,cos9j + PdSmOi +9j > -sin0/)} 

+4 cos{9g + (j)k){sin{6G + <i>k){v'k^k - w'k + 4 sin 'tpk cos 9i - 9y cos ipk cos 9i 
-9 2 sin 9i) -9i < P, cos O; - Pd sin 9i > +9[ + sin 9i) 

+ cos (^ f7 + (l)k){-v'k - y^'k'i^k - 9,008 6} - 9} < I - Pd cos 9} - P.^sm9f > — COS 0 ;)}] 
-m7]my>k sin( 6 'G + (pk) {Rx{cos ipk - sin ipk < -Ps cos 9} + Pd sin 9} >) 

+i?y(sin ijik + cos ijjk < -Ps cos 9} + pd sin 9} >) 

+R: {Py + PdCos9]+ Ps sin 9j ) + a9j {Ps sin 9} + Pd cos 9}) + Uk- 9xWk sin ipk cos 9} 
-\-9yWk <‘os V’fc cos 9} — 9z cos 9}{o. cos 9} + Vk) 

-9}(< -Wk + ei + 62 >< Ps cos 9} - Pdsm9! > 

+ < Vk + acos9}) >< 1 — PdCOs9} — PsSm9} > 4- < ei +62 >< pdCos9i 4- PsSm.9i > 
+Wk) - cos 0/ (a cos ^ 7 4- < ei 4- 62 >< PsCos9} — PdSm9r > +Vk) — sin^j) 

(< Cl 4- 62 >< Pd cos 9} 4- Ps sin 9} > 4-it;fe)}] 

-rnCmlPk cos{9g 4- ^fc) {i?i(cos ijik - sin ijik < -Ps cos 9j 4- Pd sin 9j >) 

+Hy{sm ij)k 4- cos ipk < -Ps cos 9} 4- Pd sin 9} >) 

+Rz{Pp + PdCos9j 4- /?,,sin^/) 4- a9i{PsSm9i + pdCosBi) + Uk- 9xiUkSimpkCos9j 
■¥9yWk cos i-’k cos 9} — 9z cos 9} {a cos 9] +Vk) 

-9}{< —Wk 4- Cl 4- 62 >< Ps cos 9 1 — Pd sin 9 1 > 

4- < Vk + acos9i) >< 1 - pdCos9j - PsSin9i > + < 61 + 62 >< PdC 0 s 9 i + /lasin^z > 
4 -i 7 !jt) - cos ^/(a cos 0/4- < 6 / -f 62 >< Ps cos9j - ^dsin0/ > 4-Ujfc) - sin 0/) 

(< 6 } 4- 62 >< PdCos9i 4- psSin9i > 4-t«fc)}] 


i ; 

i 
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-/?/!<(;[— cos{0a + + (l^k){f^'k^>k + 2v'k^k — w'k + 

-cos^^/(< Ot + ^1/ > sin V^it+ < Ox - > cosV^fc) - ^JzCOsOi 

+0,{< 0,cos()! - l3d sin Oj > -1) +0i{-0i < p^smOi - pdCOsO, > 

+^k < -i + Pd cos 0i + /!,, sin Oj >) - cos 0/} + cos(0g 4- <Pk){-K + Vk^l 
-2w'k<f)k - w'k^k - ^k cos 6j < 6x sin x/jk - By cos xlxk > +Bz cos 0; 

+0/ {1- < Pd cos 0/ + p, sin 0/ > ) 

+0iiPk < Ps(^os9i - Pd sin dj > +0/ < Pdsin9[ - P^cosBi > -^k) + ^fcSin0/}] 

-Imyy[sm{9a + ^fc)][sin(0G' + (^>k){i!k - + “^Wk^k + w'kpk 

-pk cos 0/ < 9a: sin ^Jk - 9y cos i)k > +02 cos 0/ + 0/(1 -< Pd cos 0; + Ps sin 0/ >) 
+9i{pk < Ps cos 0/ - pd sin 0; > +0/ < pd sin 9j - Ps cos 9j > -^k) + ^k sin 0/} 

+ cos(0G + <i>k){v'kPk + 2t;'fe4 - w'k + xu'kPl 

+ cos0/(< 0T + 9y > sin ?/>/:+ < 0x - 9y > c.QSxjjk) + ^k9z<'os9i 
■\ 9ii-PsCas9, + pds'm 0? + 1) + 0/(0/ < /i,sin0/ - pdCosOj > 

+^k < 1 - PdCosOi + PsSm9r >) + 0/cCos0/}] 

cos(0G + <^fc)][sin(0G + Pk){vk - + “^^'k^k + w'k4>k 

-^k cos 0/ < 0x sin t/^fc — 9y cos > +(L cos 0/ + 0/(1— < ,0ii cos 0/ + /?s sin 0/ >) 
+9riPk < j0,, cos0/ - PdSin9[ > +0/ < PdSmBj - Ps cos 9 j > -4>k) + <^fcsin0/} 
+ COS(0G + <Pk){VkPk + 2v'k^k - w'k + '^'k^l 
+ COS 0/(< 01 + 0t, > sin xj}k+ <9x — 9y > cos xj^k) + <^fc0z cos 9f 
+9j{-psVos9i + pdSinB] + 1) +0/(0/ < 0.,sin0/ — PdCosBj > 

+^fe < 1 - pdCosBj + pgSinB] >) + 0fcCos0/}] 

-/m,,/;[sin(0G + (j!>fc)][sin(0G + ^fc){44 + 2u'fc^fc - w* + w'k^l 
- cos 9j{< 9x + 9y> sin ?/jfc+ < 0x - 9y > cos xpk) - ^k9z cos 0/ 

■^0/(< /^scos0/ — /?dsin0/ > —1) +0/(— 0/ < PsSinB} — PdCOsBj > 
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< - 1 + /irf cos 9 1 + ps sill 9i >)- (/ifc cos 9j} + cos(6>g + { -life + Vk<i>l 

- 2 u’[ 4 h - w'f. ' 4 >k - (l>k cos < 9^ sin i/>fc - cos ■(Jjk > +9^ cos 
+///(!- < Pd<^os9i + sin0/ >) 

< 0, cos 9i - PdSm9r > +9, < /?rfsiri(?; - p,cos9j > -(^t) + ^ sin 0;}] 
-m7j„4- cos(6>f7 + ^fc)][/?.x{cost/ijfc - sin tJjki-Ps cosdj + /?dSin 0/)} 

+Rr{- sin i>k - cos i^ki-Ps cos 9, + 0^ sin 9i) 

-9 1 sin ipk {Ps sin 9i + pd cos ) } + Ry { sin -ijjk + cos -tpki-Ps cos 9j + Pd sin 0/) } 

+i?„{cos iJk - sin iik{-Ps cos 9j 4- Pd sin 9^ + 9j cos tpk{Ps sin 9j + Pd cos ^/)} 

+Rz{Pp + Pdccs 9j + P, sin 9j) + 7?. /?/(-/?,; sin 9i 4- Ps cqs9j) 

4{a6i/{/7, sin0/ 4- pdCOs9i) 4- a9'j{psCos9, - PdSin9i) + ilk} 

-Wit cos + ^i/)sin 4- {9^ - cos - Wk cos 9 1 {9^. sin ipk - 9yCos‘ipk) 
-9.vos9i{aoos9! + Vk) - Ki’k - <'7/{(-n'i(- 4- ci 4- e 2 )(A^cos<?; - pds\\\9}) 

4(t!t + a cos 9i) 4- (ci + e 2 -Vk- a cos 9j){Pd cos 9i 4- Ps sin 9j) 4- Wk} 
-9i{-Wk{PsCos9r - Pd sin 9 1 ) - 9j{-Wk 4-ei 62 ) {Ps sin 9 j - Pd cos 9 r) 

+Vk{l - PdC0s9i - Ps sin0;) 4- 9i{ei + 62 - acos0/ - Vk){-PdSin9i 4- Ps cos9t) + lOfc} 
-9i cos 9i{ei+ 0 - 2 ) {Ps sin 9] 4- Pd cos 9 ] ) 4- Vk 
-{(e, +e 2 )(-/?dsin 0 / + /?,,cos6'/) + sin0/] 

-?nCni[sin(;?G + pkMRxicosipk - sini’k{-PsC0s9j 4- ^SdSin^;)} 

+i?x{ - sin t/jfc - cos tpk{-Pii cos 0/ 4- Pd sin 9]) 

-9} simj'k{PsSin9j + Pd cos 9 j)} + R,, {sin '(pk + cos xpk{-Ps cos 9 j + Pd sin 9 j)} 

4i?j,{cos tpk — sin 'tpk{—Ps cos 9i + Pd sin 9j 4- 9i cos i)k{Ps sin 9; + Pd cos 0/)} 

+Rz{Pp + Pd cos 9i + Ps sin 9i) 4- Rz9i{-Pd sin 9i + Ps cos 9j) 

+{a9!{PsSin9i 4- PdCos9i) 4- a9]{PsCos9j — PdSin9j) 4- Uk} 

~n}kCos9i{{9s, 4- ^y)sin ipk + (4 - 9y)cos^k} - WkCOs9{{9s:Simpk - 9yCostl>k) 


I 

4. 

I 

I 

I 

if 
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-0, cosOi{n cos0j + Vk) - OzVk - (h{{-T'k + ej + e2)(/?5cos0/ — PasinOi) 

+(vk + acos0/) + (e, + ^2 - Ufc - n cos 6'r)(/?d cos 0; + ^sSini9/) + tOfc} 

cos/?/ - /Irfsin/?/) - Oi{-Wk + ej + e2)(/?5 sin - ftiCosOi) 

+!'’*.(1 - l5dCOsQi - /3s sin 0/) + /?/(ei + 62 - acos0/ - t;fc)(-/3rfsin0/ + (5, cos0/) + 10*} 
-/?/ cos 0/(6] + e2)(/3., sin 0/ + Pd cos 0/) + iifc 
-{(c] + e 2 )(-i 0 dsin 07 + /3sCos0;) + tn;fe}0/ sm0/] 



cos((9g + }{sin(6>G + (l>k){w'k + 4 sin tpk cos Oj 

- 0,1 cos fk cos Oi — 9} < /?., cos Of — sin Oj - 1 > + sin 0/) 

+ cos{do + (j>k){-v'k - - 0i < 1 - PdCosOi + PgSinOi > -cos0/}] 

+I'm.^r,[{^k sin((9G + 4>k)}{sm{^G + 4>k){v'k - Oz 

-9i <\- PdCosOj - l^gSmOi > -cosOj) 

+ cos(6^g + <t>k)[-'^'k + 0xSin'0fc - 6y cos -ipk +dj < PsCosOj - pdS'mOi + 1 > - sin^/)}] 
+/m^^[{-4 cos(0G 4- <^jfc)}{sin(6'G + 0fc)(w'fc - 9^ 

-9j < 1 — /?rf cos 9 1 — /?., sin 9j > — cos Oj) 

+ cos(0g + (t>k){-w'k + Ox sin V-’fc - Oy cos ipk + Oi < Ps cos Oj - PdsmOj + 1 > - sin 0f)}] 
+ {^fcsin(0c? + >;/>*:) }{sin(i?G + <i^>ifc)(n4 + sin cos - <?„ cos i/’jt cos 0/ 

-9t < PuCosOj — PdS'mOj - 1 > +sin0/) 

+ cos(^G + (i>k){-v'k - 9^ - Oj < I - Pd cos 9i + P, sin 9j > - cos 0/}j 
cos{6'g + 0jt)}{i?.iCOsV->fc + i?,,sin t/)jt 
-cosOpacosOj + nA;)(l 4- dj) — tCfcSin Oj]] 

+m(^m[{^k sin{9c 4- (l>k)}{Rx cos tpk 4- Ry sin 'if^k 
- cos9{{acos9i + 7Jfc)(l + 9j) — Wk sm 9j}] 

4- /??!<;(; [sin (l?G 4- (i>k){—9z + 9i{pdCos9i + /^,sin0/ - 4- 1) 4- cos0j}] 

[sin(0G 4- (i>k){-h 4- w’k) + 9i{Pd < cos 6'/ + sin 07 > 4-A < sm Oj - cos0/ > 4-1) 

4 iPp 4- Pd cos 07 4- Ps sin 07) } -f vjj. cos(0g 4- 4>k){9i cos 07 - 1)| 

+/t?i^^[cos{0g 4- <pk){9z + 9;{PdCOS0i -H PsSinOr 4- 4 4- 1) 4- cos07}] 

{sin{0G -f <^fc){nfc(l - 07 COS07)} - 1 - cos(0g 4- ^k){{-^k 4- w'k) 

~9i{Pd < COS07 + sin 07 > +ps < sin 07 — cos 07 > 4-1) — {Pp + PdCOsOi -f PgSin 07)}] 
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+ hu,,^[sm{Oa + <Pk){-(h + (hi/3dCos()i + /^,,sin /?7 - + 1) + cos0/}] 

|sin{^f7 + cosOj)} + cos{0c7 + <^k){i-^k + 'i‘’k) 

-fl/(/lrf < cos Oi + sin Oj > +/1, < sin Oj — cos Oj > + 1 ) — (/?p + p^cos 9i + /?,, sin 0;)}] 
+ 1 m + (l>k){9; + OrifiicosOi + /I, sin 9, + + 1) + cos0/}] 

[sin(^?f 7 + (f)k){-^k + ?'4) + ^liPd < «>s0/ +sin0/ > +/?* < sindj — cos9i > +1) 

+{(ip + /id cos di + ps sin 9i)} + v',, cos{9a + (j)k){9i cos 9j-l)] 

+m??m[sin{0G + <l>k){-dz + 9r{PdCos9r + ^sSin^z - 4 + 1) + cos0/}] 

(-(i/ii'A- ~ (f’l + C2+ < n. — 1 > + :rA)(l — (?/)] 

+mCm(cos(<9Gr + {Pd COS 6j + Ps sin Or + v'^ + l) + cos 9j]] 

{—OlWk ~ (ci + C2+ < 71 “ 1 > /p + •'J'a)!! ~ 9l)] 

+/7n((;[{sin(<?(; + (f)k){(h < Pd<^osOi + /i,,sin (?/ + 1 — > +cos/?/)} 

{sin((?f; + (f)k){ 03 -. cos xjjk -f 9y sin i>k)}] 

+Im„y[{cos{9o + <f>k){(h < Pdf’osOi + /i,sin(?; + l + v'f.> +cos0/)} 

{- cos (/?(7 + (?ife)(4 cos V^fc + 9y sin V^a)}] 

+Imy(;[{sm{9a + <f}k){9i < PdCosOj + PsSinOj + 1-Vk> +cos^/)} 

I - cos{9c + ^k){9x cos V>fc + 9y sin t/'a) } 

+{cos{0G + 4>k){0j < Pd cos 9} 4- /is sin ^7 + 1 + Ua > + cos 0/)}{sin(^G + (j)k) 

(9x cos V’A + 9y sin ipk)}] 

+mrj,„[{sm{9(; + (j)k)i9i < PdCOsO] + P>,sm6i + 1 - ujfe > +1)} 

{-i? 3 .sin ijfk cos 9i + Ry cos ■0* cos 0/ + (^jcos^; 4- Pdcosdj + P,sm 9i) 

{< n - 1 > If, 4 JTfc)}] 

^k){9l < Pd cos 0 1 4 Ps sill Oj 4-14- Ua> 4-1)} 

{-Rx sin 0A cos Oj + Ry cos 0^ cos Oj + {OzCOsOj + Pd cos 6j 4- Ps sin 9j) 

(< n - 1 > le + Xfc)}] 
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+/?/?« [siu(6^f? + <|!>fc)(l + + (/>fc){-02u4cos0/ + Ojv/^iPdCOS&i + 0gSmdj)} 

+ COs(/?f; + <f>k){~^z'^''k + 0iv'j^{/3d('OS0j + f3gS\n9i) — ^/w4^4}] 

+/7 H„,,[cos(/?g + <?!>fc)(l + /?/)][sin(/?(7 + (l)k) {9 cos Oj - 0iv'^{l3d cos 9 j + /?,sin0/) 
+0/u44} + cos(0G + (t>k){9zw'k cos 9j + 9iw^{/3d cos 9r + 0, sin 0/)}] 

+lm^t;[sin{0a + ?^jfc){l + ^/)][sin(^f; + (j)k) cos 9 j - 9iv[{0dcos9i + 0sSin9i) 
-¥9iw[.v'f,} + cos{9a + (l)k){9zWk cos9j + O}w'k{0dCos9j + /?, sin 9j]}] 

+/7n^[cos(6'G + (l)k)il + (9;)][sin(0G + (l)k){-92Wk cos9i + 9j'w[{0dcos9i + /?, sin^/)} 
+ ros{/?f; 4- 4ik){—9zv[.cos9i + div'f.{ftdcos0j + ft,s\\\9j) — 9l'ft^'kv'^}] 

+m?/,„[sin(0G + <j)k){l + 9 i)][{R^, cos ■)j}k + RySmipk){PsCos9] - 0dSm9r) 

-9,{li,sm9i + 0dcos9,){c] 4- C2 + ‘iasin 9i) - r/)fc(0iCOS t/’/t + 9yS\ntpk) 

+9zCOs9]{ei + 62) 

+9i{acos9j < 0.,cos9i -/Irfsini?; > -v.k) - {0p + 0dcos9j + 0,,sm9i)wk 
+a cos 9 1 {0^ cos Oj - 0d sin Oj )] 

+mC,„[cos(0G + (l>k){'i- + ^;)][(^TCOsV7jt 4- RyS\nil)k){0sCOs9i - y^dSin^j-) 

-6!{0,sin9i + 0dCos9i){ei 4- 62 4- 2nsin <?/) - Wk{9xCOS’i})k + sin'^*) 

4^jCos0/{ri 4- C 2) 

+9t[acos9i < 0scos9i - 0dS\n9i > -n.k) - {0p + 0dCos9i + 0,,sm9j)wk 
4« COS 9} {0s cos 9 1 — 0d sin 9i )] 

4/?7J^(;[sin(<?G 4 (I>k){cos0i + sm0i){-ds^,sinij>k + Oycosipk +O}{0.,cos9i — /5dsin0/)}] 
[^j; ?/4 sin {9a + <j)k) + cos(<?g 4- { -<^fc 4- 9i {0d < cos 0; 4- sin 9[ > 

+0S <sin9i - ccm9i > +1) + 0p + 0dCos9} + 0s sin 9i 4- 

4/m^;,[cos(0G + <i>k){cos9i + sin^/){-0iSin^fc + 9yCm^k + 0r(/?sCOS0/ — 0dSm9j)}] 
[siii(^G 4- 4>k){-^k 4- 9i{0d < cos9i + s\n9j > +0s < sin0/ - cos 6/ > 4-1) 

« 

+ 0d COS 9j + 0s sin 0/ — 9iv'^} + 9}w'^, cos{9a + <^ifc)] 
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+li]i,,^[sm{()a + <f>k){cosOi + ain Oi) {-Ox sin ipk + Oycostpk + Oii^sCOsO; - /?dSin0/)}] 
[sin(^?(7 + 4>k){-^k + ^liPd < cos Oi + sin ^7 > +Ps < smOi - cos 6i > +1) 

+/5p + Pd cos Of + ps sin Oj - + Ojw'i, cos{9g + <l>k)] 

+Jmrj^[cos{dG + <i>k){cos6i + sm9!){~dxSmipk +OyCos'tpk +6i{PsCos9} - PdSin 0;)}] 
[<? 7 t/ 4 sin(/?r; + (f>k) + cos{6g + <j>k){-^k + 9i{Pd < cosOi + sinOj > 

+P,, < sin Oi - cos Oj > +1) + pp + pd cos Oj 4- /?., sin Oj + ^/?4}] 

+m 7 ?„i[sin( 0<7 + (f>k){cos9i + sin 0i){-Ox sin 'ipk + 9yCostlJk + Oj{PsCos9i - pdsindj)}] 
+ (icosOf + Vk — [P/iCosOi — /Irfsiji 9i + 1)(< 7i — 1 > /^ + Xk) — (fix 4- 62 )] 
4-77).Cm[TOs(/?f7 4- (j)k){c0s9i + SlU Oj) {-Ox Simpk + 9y COS i/^k 4- OPPsCosOi - PdSlnOi)]] 
[lOfc + a cos 9j + Vk- ips cos Oj - pa sin 0/ 4- 1)(< n - 1 > t^ + Xk) - (ei 4- 62 )] 
4-/m^J{<?/ sin(0(7 4- 0fc)(cos^7 4- sin 9 }){ps cos Or - Pd sin 9 1 )} 

{(•osC/lf; 4- (l)k) con 0 1 {Ox cos ij>k 4- OySin tj^k)}] 

+ImTjn[{9j cos{9g 4- <j>k){cos9j 4- s\n 9 i){Ps cos 9 1 - PdSin9i)} 

{sin(0(7 4- <t)k) cos 9 1 {9x cos ^Jk + 9y sin T/^fc)}] 

4 -/m^[{ 07 sin( 0 G 4- <^fc)(cos0/ 4-sin0/) (/0., cos ^7 - Pd sin 9 1 )} 

(sin(£?f; 4- <pk) cos 9, {Ox cos ^Jk + Oy sin V-’t) } 

4 -{ 07 Cos( 0 g + <i>k){cos9i +sin9j){PsCos9] - PdsinOi)} 

{cos(<?f 7 4 - (kk) cos Of {Ox cos 4)k 4- Oy sin t/'fc)}] 

sin{^G 4- ^fc)(cos 9j + sin 9i){Ps cos 9i - pd sin 0/)} 

{/?, cos 9i 4- {Ox sin t/j* - Oy cos i/kk)(< ?i - I > le + a:*;)}} 

+”'Cml{^/ cos{9g + <^jfe)(cos <?7 4- siti 9i){P,cos9i - PdSin9j)} 

{/?* cos 9i + {Ox sin V’jt - 9y cos ^k)i< n-l> le + x*)}] 

~f^<d{-^kCos{9a + ^fc)}{sin(0fi' 4- <f>k){Vkh - wjb + 4 sin cos ^7 
-9y cos tpk cos 9 1 - Ox sin 9i) - 9i < ps cos 9j - pd sin Oj > 4-0/ 4- sin Oj) 



119 


+ cos(^G + (l>k){-v'k - - O^cosOi - 01 <1- pdCosOi - psSmdi > -cos0/)}] 

-Im,jr,\{^kSh\{9G + <j)k)} {siniOa + <j)k){y'k + - OzCosOj -0j <\- PaCOsOi 

> —cos/?/) + cos{0a + (f>k){yk^k — ’"jt — ^isin ^■’it<-os/?/ + 0y cos tpk cos 9 1 
+0jSin^/) — 9} < —/3s cos 0/ + f5dsm.9i +9] > — sin^/)}] 

+ /r?!,,/;[{-<^JtCOS(6/G + (?ifc)}{sin(6/G + 4>k){v'k + “ 0zCOS9j - 6} <1- l3dCOs9! 

sin 9i > — cos 9i) + cos(0g + Mi^'k^k — — 0x sin ipk cos0/ + 9y cos ipk cos 9[ 

+9z sin 9j) — 9j < —/3s cos 9] + /3rfsin 9j + 9} > - sin 0/)}] 

+{^k sin{9G + ^fc)}{sin(0G + (t>k){v'k^k - w'k + 4 sin cos 9j - 9y cos ijjk cos di 
-9z sin df) — 9i< Ps cos 9j — Pd sin 9i > +9i + sin 9r) 

+ cqs{9g + <Pk){-y'k - w'k^k - 0zcos/?/ - 9[ < 1 - /3d cos/?/ - /3s sin 0; > — cos0/)}] 
-mr],„[{-^kCOs{0G + <?^jt)}{3?3-(cosV’fc - sin < -P,o.os9j H-^^dSin/?/ >) 

I /3j,(sin j/’fc t cosi/'/t < -/3s cos/?/ | /3d sin 33/ >) 

-\-Rz{Pp + PdCosOj + pf,sm9i) + a/?/ (/3s sin 0/ + Pd cos 9 j) + Uk - 9xWk sin tpk cos 9 j 
+9yWk cos cos 9j — 9z cos 0} (a cos 9i + Vk) 

- 9 j{< -Wk + 61+62 >< PsCos 9 i - Pdsm 9 i > 

+ < Vk + a cos 9 j) >< 1 — Pd cos 9i — p,, sin 9i > 

+ < Cl + 62 >< /3d cos 0/ + Pg sin 9i > 

i Wk) - cos/?/ (a cos/?/ 1 < f'l + ('2 >< /3s<'os/?/ — /3dsin/?/ > + 1 /^) — sin/?/) 

(< C] I (’2 X /3d cos 33/ i- /3s sin 3?/ > t »»*.■)}] 

- nKml { t ^) iKin (33 G + </>*•) }{ Hi ( cosV ’* - sini/)jfc < -/3 sC <^/?/ + PdSinO} >) 

+/?„(sinV'ji- + cosV>jk < -/?,cos<?/ + /3d sin/?/ >) 

+/?z(/3p + /3d cos 0/ + /3sSin 0/) + a0/(/3sSin0/ + /3dcos0/) + «jfc - 9xV}ksintpkCOs9i 
+9yWk cos cos 0/ - 0z cos 0/ (a cos 0/ + t/jt) 

+ufk) - cos0/(acos0/+ < 61 + 62 >< PsCos Oj — Pd sin 9 1 > +wjt) 
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-shi9i){< e-i +62 >< pdCOsO] + /3,,sin0/ > +iOfc)}] 

-/7/!^ sin(6>o + ^i!>fc)][sin(6'G + A)Wk^k + 2v'k4k - Wfc + yJk<}>l 

- cos6i{< 6x + ^y> sin '4)k+ <9x — 9y > costpk) — ^k9z cos Of 

+0;(< /?., CO.S0/ — /?dsin0/ > —1) + < /?sSin0/ — (idCosOj > 

+4 < -1 + PdCosOi + PsSinOj >) - ^kCos9f} +cos{dc + ^fc){-iifc + Vk^\ 
-2w[.^k - ic'k^k - ^kCos9, < 9x sin -i/>k ~ dyoos-ipk > +0zCos9j 
+9i{l- < l^it cos 9} + fi,sm9, >) 

+9i{^k < (hcos 9j - Pdfim9, > +9] < /?rfsin0/ - P,cos9i > -^k) + i^itsin^/}] 

cos(^G + ^fc)][sm(^G + 4){iifc - Vk4>l + + w'k^k 

—^k <'os Qj < 9x sin t/’a- — 9y cos ijif. > +9^ cos 9j + 0/(1— < ftd cos 9j + /5j, sin 0/ >) 
+9ii^k < PsCos9j - ^d^m9j > +0/ < /3rfSin0/ - /3.,cos0/ > -^k) + 0ifesin0/} 

+ fos(0G + ii>k){v'k(i^k + 2v'k^k - n’A- 

+ cos0/(< 9x + 9y> sin V’’jfc+ < 0i — 9y > cos ipk) + ^k^z cos9j 
+0/(-^jCos0/ + pd^^ci9j + 1) + 0/(0; < /5.,sin0/ - /3dcos0/ > 

+^k < 1 — cos 0/ + /?s sin 0/ >) + cos 0/ }] 

sin(0G + (?!>fc)][sin(0G + (pk){Vk - + ^w'k^k + w'kh 

- 4k cos 0/ < 0s sin 4k — 9y cos 4k > +^'z cos 0/ + 0/ ( 1 — < pd cos 0/ + sin 0/ > } 
+9 1 {4k < PsCOs9i - /3dsm9} > +0/ < ;0dsin0/ - (3, cos 9} > -4k) + <;0Jt8in0/} 
+ COS(0G + 4k){Vk4k + 2v'k4k - w'k + ^k^l 

+ cM0/(< 0s + 0J, > sin V’fc+ <9x-9y> cos4k) + 4k9zCos6j 

+9i {-flg cos 0/ + (3d sin 0/ + 1 ) + 0/ (0/ < (3^ sin 0/ - (3d cos 0/ > 

+4k < 1 - (3dCosdi + PsSinOj >) + ^jfccos0/}| 

-/r?7„cf- cos(0« + <?i,fe)][sin(0G + 4k)Wk4k + ^v’k4k - Wfe + ^k^l 

-«>s0/(< 0x + 9y > sin V’jt+ < 9x — 9y> cm 4k) — 4k9z cxis 9 1 
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+(?/(< fi^cosOj — (i^sinO} > —1) + <?;(—<?; < I3,,sindi — ftdcosdi > 

4-<^s- < — 1 + fid f'os Bj + /?., sin 0j >) — 0i} + cos{6a + <t>k){ ~^k + 

-2w'f.^k - ^I’k^k - cos 6] < Or. sin t/’jt - Oy cos il>k > +0^ cos 9i 
+0i{l- < PdCosOj +■/?., sin >) 

+6i{^k < Ps^osOj - /3dsin^/ > +0j < PdSmOj - PsCosOj > -^k) + ^fesin0/}] 
-rmu\- sin(0G + ^fc)][i?i{cos V^fc - sin i)k{-0s cosOj + /?dsin 0/)} 

+Rr { - sin ii)k — cos ipk{—^.t cos 6j + 0d sin 9i) 

-9ifi\ntl’ki(is^^^^(^r 4- /Irfcosfil/)} 4- /\|sin ?/'*, + cos cos + ftdHmOj)} 
+/?y {cos V'fc - sin ipki-Ps cos 9i + pd sin Oj + Oj cos ^k{Ps sin 6[ + Pd cos di)} 
+n:iPj, + Pd <‘os Oj + /?, sin 9i) + llJi{-pd sin Or + A, cos Or) 

+{a6r{Pi,sm Or + PdC.osOj) + a(9f (/?., cos 0; - PdSmOr) + iifc} 

-Wk cos Or {{Ox -i-Oy) sin V’fc + {Ox - Oy) cos xpk} - wjk cos Or {Ox sin ipk - Oy cos 
-0xCOsOt{acos9r + Vk) - O^Vk - ^/{(-wit + e\ + e2)(/?sCos^/ - PdSin Oj) 

+{wfe + acos0/) + (e, + 62 - w* - a cos Or) {Pd cos 9 1 + PsSinOt) + wt^} 
~6r{-ibk{Ps cos Or - PdSin Or) - Oi{-Wk + Ci + e 2 ){Ps sin Or - Pd cosOj) 

+t>fr(l — pd cos Or — ps sill Or) + Oi{c\ + — a cos Oj — Vk) 

{-PdSin Or + P, cos Or) + wt} 

-Or cos Or (cj + C2) (/?5 sin Or + Pd cos 0/ ) + Ufc 
-{(cj + e2)(-PdSin9r + Ps cosOr) + Wk}0t sin^/] 

-n<m[- cos(0G + <^fc)][Hi{cos rpk - sin i^k{-Ps cos Or + Pd sin Or)} 

+i?j I - sin tj}k - COS t/)k{-Ps cos Or + Pd sin Or) 

-Or smrpk{PsSmOr + pd cos Or)} + Ry{smtpk + c(^ i>k{-Ps cos Or + Pdsm0r)} 

■» m* 

+i?„{cos - sin tpki-Ps cos Or + /?dsin Or + Or cos ^rdPa sin Or + Pd cos 9/)} 
+ih(Pp + PdCOsOt + psSinOr) + RzOr{-PdSm0i + PsCcmOr) 
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+ { aOi (/i, sin O] + fid ^os Oj ) + aO'j cos Oj - 0d sin <?/ ) + ik } 

-xDk cos 9 j {{Ox + 9y) sinipk + {Ox — Oy) cos ipk} — Wk cos Oj {Ox sin tpk — 9y cos ipk) 

-0^ cos0i{ncos9, + Vk) - Ozh - Oi{{~Wk + ei + e 2 ){ 0 sCosO[ - PdSinOj) 

+{vk + acosOj) + (cj +62 - Vk - a cos Oj){ 0 d cos Or + /?.,sin0/) + «;*} 

-01 {-tbk{ 0 s cos Or — pdS'in Or) — Or{—Wk + Ci + e2)(/0sSin Or — PdCosOr) 

+n{l — ftdCosOi — /?,,sin 0 /) + 0t{(’\ + 62 — n.cos0r — ?+){— + ftgCOsOr) + w*} 
-Or cosOr{e\ + 62 )(/?., sin <?/ + 0dCosOj) + Vk 
-{(f’l + <’2){-0,i^io 0, + ftg cosOj) + Wk}0i sin 0,] 
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The section integrals are defined as: 


m 


(A.l) 

m.7]m 

= 1 I^PVdrjd^ 

(A.2) 

rnCm 


(A.3) 


= ll^fK'^dijdC 

(A.4) 



(A.5) 

Irrtrj^ 

= I 

{A.6) 


Where m is the mass per unit length of the blade; m7]m and are the first 
moment of the mass of the cross-section about the elstic axis; /m^, and 
are the mass moments of inertia per unit length of the beam about the elastic axis. 



Appendix B 


MODIFICATION OF THE 
TERM [Vg^] ASSOCIATED 
WITH KINETIC ENERGY 

Expanding sm{0a + <Pk) cos{0a + ^k), and assuming that 4>k is small- 

sin -1- (j)k) ^ sin Oc + (pk cos 9g 
cos{9g + (f^k) ~ (X)s 9a - 4>k sin 9g 

Sulwtituting the above approximation in j , and rewriting as 

[vS] « [v^;*] + 

where 

^] = ^ [{ - /m^tj)(cos^ 9a - sin^ 9 g ) cos^ 9 j 

+4/m,^ sin 9a cos 9g cos^ 9i + (/m^^ cos 9a — sin 9a) sin 0/} {<!>,} dx 

[K/i#] = J sin9acos9Gcm^ 9i - Im^{sin^ 9a — cos^ 9a) cos^9i 

4- (/m^^ sin 9g + cos 9g) sin 9i} {#,}! dx 
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The .subinatrix [A 33 #] Ii<xs (,o be added U) (.he matrix of linear stiffness matrix, 

given in Sec. 6.2.3. 




